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Preface 

This was the fourth workshop in our series of annual "retreats" of the Sonderforschungs- 
bereich Astroteilchenphysik (Special Research Center for Astroparticle Physics), or SFB for 
short, to the Ringbcrg Castle above Lake Tegernsee in the foothills of the Alps. These 
meetings are meant to bring together the members of the SFB which arc dispersed between 
four institutions in the Munich area, the Technical University Munich (TUM), the Ludwig- 
Maximilians-University (LMU), and the Max-Planck-Institute for Physics (MPP) and that 
for Astrophysics (MPA). We always invite a number of external speakers, including visitors 
at our institutions, to complement the scientific program and to further the exchange of ideas 
with the international community. 

This year's topic was "Neutrino Astrophysics" which undoubtedly is one of the central 
pillars of astroparticle physics. We focused on the astrophysical and observational aspects 
of this field, deliberately leaving out theoretical particle physics and laboratory experiments 
from the agenda. Each day of the workshop was dedicated to a specific sub-topic, ranging 
from solar, supernova and atmospheric neutrinos over high-energy cosmic rays to the early 
universe. A session on future prospects served to conclude the workshop and provide an 
outlook on the field in the next decade and beyond. We started every topical session with one 
or two introductory talks, reviewing the status of theory and experiment and to providing 
some background for the non-experts. 

For the entire program we interpreted "neutrino astrophysics" in a broad sense, including, 
for example, the physics of 7-ray bursts or the recent observations of TeV 7-rays by the 
imaging air-Cherenkov technique. Some of the after-dinner-talks went significantly beyond a 
narrow interpretation of the field! From our perspective the profile of neutrino astrophysics 
as defined by our program worked very well, even better than we had hoped. We are proud 
that the main complaint of the participants seemed to be that they did not get enough 
mountain-hiking done because the sessions were too interesting to miss. 

Besides regular SFB resources this workshop was made possible by a direct grant from 
the Max-Planck-Society and additional funds from the Max-Planck-Institute for Astrophysics. 
Special thanks go to the SFB secretary, Maria Depner, for her smooth and skillful management 
of all practical matters related to the workshop. 

We thank the participants for their high-level contributions and for being extremely co- 
operative in submitting the "extended abstracts" of their talks on time and in a format that 
allowed us to produce these proceedings in electronic form. Anyone interested in a printed 
version should write to the SFB secretary at the address given on the previous page. We 
hope that you will find this booklet a useful and up-to-date resource for the exciting and 
fast-developing field of neutrino astrophysics. 

Michael Altmann, Wolfgang Hillebrandt, Hans-Thomas Janka and Georg Raffelt 
Munich, January 1998 
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History of Neutrino Physics: Pauli's Letters 

Rudolf L. Mdfibauer 

Physik Department El 5, Technische Universitdt Miinchen, 85747 Garching, Germany 



Editors' Note 

Professor Mofibauer gave an evening lecture on the history of neutrino physics during which 
he read PauU's famous letters on the neutrino hypothesis. Rather than writing a formal con- 
tribution to our proceedings, Professor MoBbauer suggested that we reproduce these letters. 
(They were taken from Ref. but can also be found in Ref. ||2|.) Because much of the 
humour of Pauli's writing is lost in the translation we quote the original German text. An 
English translation can be found, for example, in Ref. |^]. 




Figure 1: Wolfgang Pauh (1900-1958) in Vienna 1933. 
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Pauli's Letters 

Brief an Oskar Klein, Stockholm, vom 18. 2. 1929 

Aber ich verstehe zu wenig von Experimentalphysik um diese Ansicht beweisen zu konnen 
und so ist Bohr in der fiir ihn angenehmen Lage, unter Ausnutzung meiner allgemeinen 
Hilflosigkeit bei der Diskussion von Experimenten sich selber und mir unter Berufung auf 
Cambridger Autoritaten (iibrigens ohne Literaturangabe) da etwas beliebiges vormachen zu 
konnen. 

Brief an Oskar Klein, Stockholm, 1929 

Ich selbst bin ziemhch sicher (Heisenberg nicht so unbedingt), daB 7-Strahlen die Ursache 
des kontinuierlichen Spektrums der /?-Strahlen sein miissen und daB Bohr mit seinen dies- 
beziiglichen Betrachtungen iiber eine Verletzung des Energiesatzes auf vollkommen falscher 
Fahrte ist. Auch glaube ich, dafi die warmemessenden Experimentatoren irgendwie dabei 
mogeln und die 7-Strahlen ihnen nur infolge ihrer Ungeschicklichkeit bisher entgangen sind. 

Brief an die Gruppe der "Radioaktiven" 1930 

Physikalisches Institut 

der Eidg. Technischen Hochschule 

Ziirich Ziirich, 4. Dez. 1930 

Liebe Radioaktive Damen und Herren! 

Wie der Uberbringer dieser Zeilen, den ich huldvollst anzuhoren bitte, Ihnen des naheren aus- 
einandersetzen wird, bin ich angesichts der falschen Statistik der N- und Li 6-Kerne, sowie des 
kontinuierlichen /3-Spektrums auf einen verzweifelten Ausweg verfallen, um den Wechselsatz 
der StatistikQ und den Energiesatz zu retten. Namlich die Moglichkeit, es konnten elektrisch 
neutrale Teilchen, die ich Neutronen^ nennen will, in den Kernen existieren, welche den Spin 
1 /2 haben und das Ausschliefiungsprinzip befolgen und sich von Lichtquanten aufierdem noch 
dadurch unterscheiden, dafi sie nicht mit Lichtgeschwindigkeit laufen. — Das kontinuierliche 
/3-Spektrum ware dann verstandlich unter der Annahme, daB beim /?-ZerfaIl mit dem Elektron 
jeweils noch ein Neutron emittiert wird, derart, dafi die Summe der Energien von Neutron 
und Elektron konstant ist. 

Nun handelt es sich weiter darum, welche Krafte auf die Neutronen wirken. Das wahr- 
scheinlichste Modeh fiir das Neutron scheint mir aus wellenmechanischen Griinden dieses 
zu sein, dafi das ruhende Neutron ein magnetischer Dipol von einem gewissen Moment fi 
ist. Die Experimente verlangen wohl, dafi die ionisierende Wirkung eines solchen Neutrons 
nicht grofier sein kann als die eines 7-Strahls, und dann darf fi wohl nicht grofier sein als 
e • (10~^^ cm). Ich traue mich vorlaufig aber nicht, etwas iiber diese Idee zu publizieren, und 
wende mich erst vertrauensvoll an Euch, liebe Radioaktive, mit der Frage, wie es um den 
experimentellen Nachweis eines solchen Neutrons stande, wenn dieses ein ebensolches oder 
etwa lOmal grofieres Durchdringungsvermogen besitzen wiirde wie ein 7-Strahl. . . . 

Also, liebe Radioaktive, priifet, und richtet. — Leider kann ich nicht personlich in Tiibingen 
erscheinen, da ich infolge eines in der Nacht vom 6. zum 7. Dez. in Ziirich stattfindenden 
Balles hier unabkommlich bin. . . . Euer untertanigster Diener W. Pauli 

^Heute Pauli'sches Ausschliefiungsprinzip 
^Heute Neutrinos 
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Brief an Oskar Klein, Stockholm, vom 12. 12. 1930 

Ich kann mich vorlaufig nicht entschliefien, an ein Versagen des Energiesatzes ernstlich zu 
glauben und zwar aus folgcndcn Griindcn (von denen ich natiirlich zugebe, daB sie nicht 
absolut zwingend sind). Erstens scheint es mir dafi der Erhaltungssatz fiir Energie-Impuls 
dem fiir die Ladung doch sehr weitgehend analog ist und kann keinen theoretischen Grund 
dafiir sehen, warum letzterer noch gelten sollte (wie wir es ja empirisch iiber den /3-Zerfall 
wissen) wenn ersterer versagt. Zweitens miifite bei einer Verletzung des Energiesatzes auch 
mit dem Gewicht etwas sehr merkwiirdiges passieren. 

Telegramm von Reines und Cowan vom 14. 6. 1957 an Wolfgang Pauli 

We are happy to inform you that we have definitely detected neutrinos from fission fragments 
by observing inverse beta decay of protons. Observed cross section agrees well with expected 
6 • 10-''^ cm2. 

References 

[1] W. Pauli, Fiinf Arbeiten zum Ausschliefiungsprinzip und zum Neutrino, Texte zur 
Forschung Vol. 27 (Wissenschaftliche Buchgesellschaft Darmstadt, 1977). 

[2] W. Pauli, Wisscnschaftlichcr Bricfwechsel mit Bohr, Einstein, Hcisenberg, u.a.. Vol. II: 
1930-1939, ed. by K. v. Meyenn, (Springer- Verlag, Berlin, 1985). 

[3] W. Pauli, On the Earlier and More Recent History of the Neutrino (1957) in: Neutrino 
Physics, ed. by K. Winter (Cambridge University Press, 1991). 
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What Killed The Dinosaurs? 

Arnon Dar^'^ 

^ Max-Planck-Institut fiir Physik ( Werner- Heisenberg-Institut) 
Fdhringer Ring 6, 80805 Miinchen, Germany. 
"^Department of Physics and Space Research Institute, 
Technion, Israel Institute of Technology, Haifa 32000, Israel. 

The early history of hfe during the Precambrian until its end 570 million years (My) ago 
is poorly known. Since then the diversity of both marine and continental life has increased 
exponentially. Analysis of fossil records shows that this diversification was interrupted by five 
massive extinctions and some smaller extinction peaks [|^. The largest extinction occurred 
about 251 My ago at the end of the Permian period. The global species extinction ranged then 
between 80% to 95%, much more than, for instance, the Cretaceous-Tertiary extinction 64 My 
ago which killed the dinosaurs and claimed ~ 47% of existing genera |Q] . In spite of intensive 
studies it is still not known what caused the mass extinctions. Many extinction mechanisms 
have been proposed but no single mechanism seems to provide a satisfactory explanation of 
the complex geological records on mass extinctions These include terrestrial mechanisms 
such as intense volcanism, which coincided only with two major extinctions [Q] or drastic 
changes in sea level, climate and environment that occurred too often, and astrophysical 
mechanisms, such as a meteoritic impact that explains the iridium anomaly which was found at 
the Cretaceous/Tertiary boundary but has not been found in any of the other extinctions, 
supernova explosions and gamma ray bursts Q which do not occur close enough at a 
sufficiently high rate to explain the observed rate of mass extinctions. 

The geological records, however, seem to indicate that an accidental combination of drastic 
events |^] occurred around the times of the major extinctions. For instance, the dinosaur 
extinction coincides in time with a large meteoritic impact, with a most intensive volcanic 
eruption and with a drastic change in sea level and climate. The origin of these correlations is 
still unclear. Meteoritic impacts alone or volcano eruptions alone or sea regression alone could 
not have caused all the major mass extinctions. An impact of a 10 km wide meteorite with 
a typical velocity of 30 km s~^ was invoked Q in order to explain the Cretaceous- Tertiary 
(K/T) mass extinction 64 My ago, which killed the dinosaurs, and the iridium anomaly 
observed at the K/T boundary. But neither an iridium anomaly nor a large meteoritic 
crater have been dated back to 251 My ago, the time of the Permian/Triassic (P/T) mass 
extinction, which was the largest known extinction in the history of life where the global 
species extinction ranged between 80% to 95%. The gigantic Deccan volcanism in India that 
occurred around the K/T boundary Q and the gigantic Siberian basalts flood that occurred 
around the P/T boundary have ejected approximately 2 x 10^ km^ of lava They were more 
than a thousand times larger than any other known eruption on Earth, making it unlikely 
that the other major mass extinctions, which are of a similar magnitude, were produced 
by volcanic eruptions. Although there is no one-to-one correspondence between major mass 
extinctions, large volcanic eruptions, large meteoritic impacts, and drastic environmental 
changes, there are clear time correlations between them. We propose that near encounters 
of Earth with "visiting planets" from the outer solar system are responsible for most of the 
mass extinctions on planet Earth and can explain both the above correlations and the detailed 
geological records on mass extinctions 
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Recent observations with the Hubble Space Telescope of the Helix Nebula (Fig. 1) the 
nearest planetary nebula, have discovered that the central star is surrounded by a cir- 
cumstellar ring of about 3500 gigantic comet-like objects ("Cometary Knots") with typical 
masses about 10~^Mq, comparable to our solar system planets (MEarth ~ 3 x 1O~^M0 and 
Afjupiter ~ 9-6 X 1O~^M0). It is not clear whether they contain a solid body or uncollapsed 
gas. They are observed at distances comparable to our own Oort cloud of comets but they 
seem to be distributed in a planar ring rather than in a spherical cloud like the Oort cloud. 
It is possible that these Cometary Knots have been formed together with the central star 
since star formation commonly involves formation of a thin planar disk of material possessing 
too high an angular momentum to be drawn into the nascent star and a much thicker outer 
ring of material extending out to several hundred AU. Evidence for this material has been 
provided by infrared photometry of young stars and also by direct imaging. 

It is possible that such Cometary Knots and the recently discovered gigantic asteroids |jl^ 
in the outer solar system between the Kuiper belt and the Oort cloud are the high mass end 
of the vastly more numerous low mass comets. The massive objects are more confined to 
the ecliptic plane because of their relatively large masses, and form a circumstellar ring, 
while the very light ones are scattered by gravitational collisions into a spherical Oort cloud. 
Gravitational interactions in the ring can change their parking orbits into orbits which may 
bring them into the inner solar system. In fact, various "anomalies" in the solar planetary 
system |]Tl| could have resulted from collisions or near encounters with such visitors in the 



early solar system. These include |12] the formation of the moon, the large eccentricity and 
inclination of some planetary orbits, the retrograde orbits of 6 moons of Jupiter, Uranus 
and Neptune and the tilted spin planes of the Sun, the planets and moons relative to their 
orbiting plane. 

Strong gravitational tidal forces can cause frictional heating of planetary interiors leading 
to strong volcanic activity, as seen, for instance, on Jupiter's moon lo, the most volcanically 
active object known in the solar system. The moon, the Sun and the known planets are too 
far away to induce volcanic eruptions on Earth, but relatively recent "visits" of planet-like 
objects near Earth could have generated gigantic tidal waves, large volcanic eruptions, drastic 
changes in climate and sea level, and impacts of meteorites which were diverted into a collision 
course with Earth by the passage through the astroids and Kuiper belts |jl3|. Thus, visiting 
planets may provide a common origin for the diverse mass extinction patterns as documented 
in the geological records. 

Although exact calculations of surface tidal effects are a formidable scientific effort that 
was begun by Newton and has continued and improved since then by many of the great math- 
ematicians and physicists, an approximate estimate of the flexing (h) of Earth (radius Re) by 
a passing planet (distance d) can be easily obtained by neglecting the rotation of Earth and 
the speed of the passing planet (mass Mp) and by assuming quasi hydrostatic equilibrium: 



4Me \ d 

The maximal crustal tide due to the moon is 27 cm. However, a visiting planet with a typical 
mass like that of the Cometary Knots which passes near Earth at a distance comparable 
to the Earth-Moon distance produces gigantic oceanic and crustal tidal waves which are a 
few hundred times higher than those induced by the moon. Oceanic tidal waves, more than 
1 km high, can flood vast areas of continental land and devastate sea life and land life near 
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continental coasts. The spread of ocean waters by the giant tidal wave over vast areas of land 
and near the polar caps will enhance glaciation and sea regression. 

Flexing the Earth by /i ~ 100 m will deposit in it ~ aGM^h/R^ « lO^"^ ergs, where 
a ~ 0.1 is a geometrical factor. It is approximately the heat release within Earth during 
10^ y by radioactive decays. The flexing of Earth and the release of such a large energy in 
a very short time upon contraction might have triggered the gigantic volcanic eruptions that 
produced the Siberian basalts flood at the time of the P/T extinction and the Deccan basalts 
flood at the time of the K/T extinction. 

A reliable estimate of the masses and the flux of the visiting planets/planetesimals is not 
possible yet. However, we have fixed them from the assumption that the unaccounted energy 
source of Jupiter and its tilted spin plane relative to its orbital plane are both due to accretion 
of visiting planets/moons. Prom the 3.13° tilt of Jupiter's spin and the accretion rate, we 
inferred that Nj 16 planets of average mass Mp 0.5Me have crashed into Jupiter during 
its ~ 4.57 Gy lifetime. Similar estimates for other planets, although yielding the correct order 
of magnitude, are less reliable because the inferred number of accreted planets is too small. 
The 7° tilt of the spin of the Sun could have been produced by the impact of ~ 3 x 10^ such 
planets (ignoring possible angular momentum loss by the solar wind). This means that the 
Sun has accreted ~ 5% of its mass after its formation, at a rate of ~ 7 planets per My. In 
each capture episode, ~ 6 x 10^^ erg of gravitational energy is released in the Sun's convective 
layer. It produces optical and x-ray flashes at a rate ~ 7 x 10~^ Lg^y"^ for Sun-like stars. It 
also causes a significant luminosity rise for an extended time which may have induced climatic 
and sea level changes on Earth, and extinctions of species which could not have adapted to 
large environmental changes. The predicted rate is consistent with the observed rate of large 
changes in O^^ concentration in sea water sediments which record large changes in sea water 
level and total glacier volume. 

Using our inferred planet flux from Jupiter and its collimation by the Sun, we obtained that 
a "visiting rate" of once every = 100 My for planets with Mp ^ 0.5Me which fall towards 
the Sun implies a passing distance of approximately 170,000 km from Earth which produces 
crustal tidal waves of h > 125 m and water tidal waves of ~ 1 km height. The combination 
of tidal waves, volcanic eruptions, meteoritic impacts and environmental and climatological 
changes can explain quite naturally the biological and time patterns of mass extinctions. Por 
instance, the giant tidal waves devastate life in the upper oceans layers and on low lands 
near coastal lines. They cover large land areas with sea water, spread marine life to dry on 
land after water withdrawal, and sweep land life into the sea. They flood sweet water lakes 
and rivers with salt water and erode the continental shores where most sea bed marine life 
is concentrated. Amphibians, birds and inland species can probably survive the ocean tide. 
This may explain their survival after the K/T extinction. Survival at high altitude inland 
sites may explain the survival of some inland dinosaurs beyond the K/T border. Volcanic 
eruptions block sunlight, deplete the ozone layer, and poison the atmosphere and the sea with 
acid rain. Drastic sea level, climatic and environmental changes inflict further delayed blows 
to marine and continental life. But high-land life in fresh water rivers which are fed by springs, 
that is not so sensitive to temperature and climatic conditions, has better chances to survive 
the tidal waves, the volcano poisoning of sea water, and the drastic sea level, environmental 
and climatic changes. 

Altogether, visiting planets offer a simple and testable solution to the puzzling correlations 
between mass extinctions, meteoritic impacts, volcanic eruptions, sea regression and climatic 
changes, as documented in the geological records. Perhaps the best test will come from the 



9 



MACHO sample of ~100,000 light curves of variable stars: Planets/moons crashing onto stars 
produce a big flash of light and soft x-rays from the hot spot at the impact point on main 
sequence (rotating) stars, a nova-like thermonuclear explosion on the surface of a white dwarf 
(rare), or a strong gamma ray flash from the surface of a neutron star (very rare). 
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Figure 1: A small section of a distant ring of gigantic comet-like objects around the central 
star in the Helix nebula, the nearest planetary nebula, as seen recently by the Hubble Space 
Telescope [9]. These objects have a typical mass of SM^;. Their gaseous heads have a typical 
size of 10^^ cm. It is not clear whether they contain a solid body or uncollapsed gas. 
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Solar Models 

M. SUx 

Kiepenheuer-Institut fiir Sonnenphysik, Freiburg, Germany 



Introduction 

A number of solar model calculations has been presented in recent years, especially in view of 
the predicted flux of neutrinos from the Sun |5|, ^, |9|, |2^. In the present contribution 

I shall discuss mostly the Standard Solar Model. This model, including some refinements, 
is confirmed by seismological tests. Non-standard models, on the other hand, face severe 
difficulties. I close with a remark on the Sun's magnetic field in the neutrino context. 



The Standard Model 

The standard solar model is a gas sphere in hydrostatic equilibrium. This may seem trivial. 
However, some of the earliest non-standard models were based on the assumption of a very 
strong internal magnetic field, or of rapid rotation of the solar core; both of these assumptions 
violate the spherical symmetry as well as the hydrostatic equilibrium. 

The input to the standard model should be to the best of our present knowledge. Let me 
begin with the mass, luminosity, and radius, 

m© = (1.989±0.001)-103°kg , Lq = (3.845±0.008) • lO^^W , r© = (6.960±0.001) • lOV . 

The value given for the luminosity is from the ACRIM experiment on the Solar Maximum 
Mission, with the error as quoted in [^. The radius is for the optical depth r = 2/3 where 
the temperature is equal to the effective temperature (LQ/AirarQ)^/'^. 

The age of the Sun is found by measuring the decay of long-lived radioactive isotopes in 
meteorites, 

tQ = (4.57 ± 0.05) • lOVears . 

Wasserburg, in an appendix to gives an even smaller error; the main source of uncertainty 
is the not well-known state of the formation of the Sun at the time of the melting and 
crystallization of the meteoritic material. 



a) Chemical Composition 

The composition of the solar surface (and the entire perfectly mixed outer convection zone) is 



derived from solar spectroscopy and from meteoritic data. Recent results [16] show a slightly 
decreased content of the heavy elements in comparison to earlier results, e.g. ||l|. For oxygen, 
the most abundant element after H and He, the change is —0.06 in the usual logarithmic 
scale, or ~ —13%. Altogether, the ratio of the mass fractions Z and X of the heavy elements 
and hydrogen is 0.0245, instead of the earlier value 0.0267. 

For helium the most accurate determination is based on helioseismology. This is because 
both the sound velocity and the acoustic cut-off frequency depend on the adiabatic expo- 
nent Fi: 

2 ^Fi fdlnT' 

c = , Ll>A 



2c ' \dlnP 
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In the depth range where an abundant element is partly ionized Fi is greatly affected by the 
energy of ionization. Thus the He abundance, in particular by its effect in the zone of partial 
Hell ionization, has some influence on the eigenfrequencies of the Sun's pressure (p) modes 
of oscillation. The inversion of observed frequencies yields a mass fraction I5 = 0.242 it 0.003 
( |]2^ , other authors find similar results). It should be noted that the original He abundance, 
Yq, is adjusted so that the luminosity of the present model, at age i©, equals Lq. The result 
of this procedure is Yq = 0.27 . . . 0.28, depending on other input to the model. The difference 
to the seismically determined Ys is appropriate in view of helium settling in the radiative 
solar interior, see below. 



b) Nuclear Reactions 

Concerning the nuclear reactions I shall concentrate on the pp chains, which provide ~ 99% of 
the energy. The reaction rates, and therefore the branching between the three chains leading 
to helium, depend on temperature (see below), and on the "astrophysical S-factor" S{E). 
This factor depends weakly on the center-of-mass energy E and measures the cross section 
after separation of 1/E times the penetration probability through the Coulomb barrier. For 



the most important reactions Parker |21] reviews results of 5- factors (at zero energy): 



p(p,e+z^)d 5ii(0) = 3.89 • lO'^s MeV-b ±1% 

3He(3He,2p)4He 533(0) = 5.0 MeV-b ±6% 

3He(4He,7)7Be ^34(0) = 533 eV-b ±4% 

^Be(p,7)8B 5i7(0) = 22.2 eV-b ±14% 







It is difficult to assess the errors. Sn is so small that it can only be calculated. The other 
S'-factors are measured in the laboratory but must be extrapolated to zero energy, and a 
correction must be applied for electron screening at low energy. In particular the 14% error 



of 5i7 has been criticized, e.g. |17], as it is derived from diverse experiments with partially 
contradicting results in the range 17.9-27.7 eV-b. Perhaps better results will soon become 
available as the measurements are extended to lower energy, cf. the contribution of M. Junker 
to these proceedings. 



c) Equation of State and Opacity 

It is appropriate to discuss these two important ingredients to the solar model together, since 
both depend on the knowledge of the number densities of the diverse particles, and especially 
on the ionization equilibria and the electron density. To first order the Sun consists of a 
perfect gas; but significant corrections, at the percent level, arise e.g. from the electrostatic 
interaction of the particles (especially at the depth where abundant species are partially 
ionized) and from partial electron degeneracy (in the core), cf. [23|. Modern standard models 



are usually calculated with a tabulated equation of state and opacity. The most recent tables 
from the Lawrence Livermore Laboratory |18] still exhibit unexplained discrepancies of up to 



20% as compared to various other calculations; in the energy- generating region of the Sun 
the uncertainty probably is much less, 2.5% according to 0]. 



Generally, the recent opacities [18| are somewhat increased in comparison to earlier results, 
e.g. from Los Alamos, mainly because more elements have been included in the calculations. 
The increase renders the radiative transport of energy less effective. The solar model responds 
with a slightly raised central temperature (cf. discussion below), and a slightly lowered base of 
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the outer convection zone. As for the depth of the convection zone, there are two other input 
modifications that may lead to an increase: convective overshoot (see below) and element 
diffusion. 



d) Element Diffusion 

Driven by the gradients of pressure, temperature, and composition, helium and the heavy 
elements diffuse toward the solar center, while hydrogen diffuses upward. The process is 
slow, with a characteristic time exceeding the Sun's age by a factor 100 or more. Nevertheless 
diffusion should be included into the standard solar model since it produces significant effects, 
especially in view of the details that can be seen by helioseismology. The following table 
summarizes some results [01: 



Diffusion 


Yo 


Ys 


n/rQ 


Tc{W^) K 


snu(Cl) 


snu(Ga) 




0.268 


0.268 


0.726 


15.56 


7.0 


126 


He 


0.270 


0.239 


0.710 


15.70 


8.1 


130 


He &: heavy el. 


0.278 


0.247 


0.712 


15.84 


9.3 


137 



The obvious effect of a decreased surface mass fraction Ys of He is essentially in agreement with 
the results of other authors, and with the helioseismological result (see above); the increase of 
Ys obtained if heavy element diffusion is included is due to the larger initial helium content, 
Yq, which is required in this case for the luminosity adjustment. Also, the central temperature 
rises: Helium diffusion changes the mean molecular weight, which must be compensated by a 
higher temperature in order to maintain the hydrostatic equilibrium; heavy element diffusion 
increases the opacity, which must be compensated by a steeper temperature gradient (and 
hence larger Tc) to maintain the radiative transport of energy. Together with Tc, the neutrino 
rates predicted for the chlorine and gallium experiments are increased, as listed in columns 
6 and 7 of the table. It should be noted that the Tc effect found in other calculations is 
somewhat smaller, cf. the contribution of Schlattl and Weiss to these proceedings. Also, the 
effect of helium diffusion on the depth of the convection zone (column 4) is in contrast to the 



result of Cox et al. |12]. 



The depth of the convection zone has been determined from helioseismology by several au- 
thors. Below the convection zone the temperature gradient is determined by the requirement 
of radiative energy transport; within the convection zone the gradient is nearly adiabatic be- 
cause of the large heat capacity and hence large effectivity of the convective transport. Thus, 
a prominent transition occurs in dT/dr and, accordingly, in d(? jdr. As the sound speed c(r) 
can be obtained from the p-mode frequencies, the transition can be located. Christensen- 
Dalsgaard et al. [|l^] find rb/r© = 0.713 it 0.003, a more recent study Q yields an even smaller 
uncertainty, 0.713 it 0.001, with the claim that systematic errors are included! — The tem- 
perature at the base of the convection zone slightly depends on the helium abundance; for 
the value of Ys quoted above a good estimate [|lO| is (2.22 it 0.05) • 10^ K. 



Neutrinos 

The standard model of the Sun predicts the flux of neutrinos, as a function of neutrino energy, 
that should be observed on Earth. So far 5 detectors have been used for such observations: 
One using chlorine, two using gallium and two using water. The detectors as well as the 
observational results are reviewed elsewhere in these proceedings. Briefly the results are 
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as follows: (1) All the detectors measure a significant neutrino flux, and so confirm that 
nuclear energy generation actually takes place in the Sun. (2) The results of the two gallium 
experiments agree with each other, and the results of the two water experiments agree with 
each other. (3) All experiments measure a flux that lies significantly below the prediction of 
the standard solar model. (4) The deficit depends on neutrino energy. This last point is of 
particular interest for the solar model builder, because it is for this reason that it appears to 
be impossible to repair all the deficits simultaneously by modifications of the solar model. 

The significance of the neutrino deficit and its energy dependence has been illustrated by 
a calculation of 1000 standard models Q based on input parameters having normal distribu- 
tions with appropriate means and standard deviations. For the diverse neutrino experiments 
these 1000 models predict results that have a certain spread but are clearly in conflict with 
the actual measurements. If the predicted flux of energetic neutrinos originating from the 
beta decay of is replaced by the value obtained in the water experiments (~ 42% of the 
predicted), then the predictions for the chlorine and gallium experiments become smaller 
but are still significantly above the measurements. In other words, the existing experiments 
cannot simultaneously be reconciled with the standard solar model. 



In two recent calculations |13, 14 1 all the input parameters have been pushed to the extreme 



in order to minimize the neutrino flux; nevertheless the prediction is still significantly too high. 
Non-Standard Models 

The aim of most non-standard solar models is to lower the temperature in the energy- 
generating central region and thereby change the branching ratios of the three pp chains, 
in particular in order to suppress a part of the high-energy neutrino flux. Castellani 
et al. Q discuss in detail the influence of the diverse input parameters on the central temper- 
ature of the Sun. For example, a 45% increase of Su, or a 50% decrease of Z/X, or a 29% 
decrease of the opacity would result in a 4% smaller central temperature. For the neutrino 
fluxes resulting from ^Be and *B the temperature dependence is 0] 



$^(^Be) oc , $^(^B) oc T, 
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Hence, a 4% cooler solar core reduces the predicted ^B neutrino flux to 0.48 of the standard 
value, and the flux of ^Be neutrinos to 0.72. This may help to remove the discrepancy for the 
water experiments, and to reduce the discrepancy for the chlorine experiment, but it provides 
almost no help for the gallium experiments. Of the 137 snu (above table, last line) there 
are only 16 from the decay of ^B, but 38 from the electron capture of ''Be. Thus a large 
discrepancy remains; more specifically, the gallium experiments appear to leave no room for 
the predicted ''Be neutrinos. This is the major difficulty of the non-standard solar models. 

Perhaps the neutrino discrepancy will finally be resolved by a combination of various 
effects. A slight decrease of the central solar temperature Tc may be one of these effects, 
although it is entirely unclear at present how such a decrease of Tc could be achieved. Heavy 
element diffusion increases Tc, as we have seen. Mixing of the solar core apparently helps, 
but the mixed-core model seems to fail the seismological test [^]. Other handles, such as 
the opacity or the equation of state, permit only variations of Tc that are too small for a 
substantial effect. 

The conclusion is that most of the discrepancy should rather be resolved by non-standard 
neutrinos. The energy-dependent conversion of electron neutrinos into other neutrino flavours 



by the Mikeyev-Smirnov-Wolfenstein effect is a possibility; for a recent review see |17] 
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Convective Overshoot and Magnetism 

Non-local versions of the mixing-length theory allow for overshooting flows at the base of 



the solar convection zone |24, 25]. The overshoot layers calculated so far have a nearly adia- 
batic temperature gradient. The thickness of such a layer should not exceed a few thousand 
kilometers, according to helioseismological results [|l^]. On the other hand, some convective 
overshoot will certainly occur and must be included into a calculation that is designed to 



exactly reproduce the depth of the convection zone. The following table, adapted from |25|, 
gives a few examples of standard models, calculated without element diffusion. The input 
consists of the opacity, the type of convection formalism (MLT: mixing-length theory, CM: 
Canuto & Mazzitelh), equation of state (SS: Stix k Skaley (2|], MHD: Mihalas et al. |9|), 
initial helium mass fraction Yq, mixing- length parameter a. The results are the relative radius 
rft/r© where convection ceases, the central temperature Tc, and the predicted neutrino rates 
for the CI and Ga experiments. There is very little variation of these neutrino rates. 



Nr 


Opacity 


MLT 


EOS 




^0 




a 




Tc6 


snu (CI) 


snu (Ga) 


1 


Alamos 


local 


SS 





2733 


2 


403 


0.7361 


15.56 


7.51 


128.3 


2 


Alamos 


local 


MHD 





2731 


2 


164 


0.7325 


15.55 


7.49 


128.2 


3 


Alamos 


nonl. 


SS 





2733 


2 


409 


0.7210 


15.56 


7.51 


128.4 


4 


Alamos 


nonl. 


MHD 





2732 


2 


168 


0.7201 


15.55 


7.49 


128.2 


5 


Opal92 


local 


SS 





2769 


1 


847 


0.7231 


15.56 


7.68 


129.4 


6 


Opal92 


local 


MHD 





2768 


1 


677 


0.7194 


15.56 


7.66 


129.3 


7 


Opal92 


nonl. 


SS 





2770 


1 


850 


0.7114 


15.57 


7.68 


129.4 


8 


Opal92 


nonl. 


MHD 





2768 


1 


679 


0.7096 


15.56 


7.66 


129.3 


9 


Opal95 


local 


Opal95 





2774 


1 


660 


0.7220 


15.59 


7.81 


130.0 


10 


Opal95 


nonl. 


Opal95 





2774 


1 


661 


0.7116 


15.59 


7.81 


130.0 


11 


Opal95 


CM 


Opal95 





2775 


1 


033 


0.7218 


15.58 


7.80 


129.9 


12 


Opal92 


nonl. 


MHD 





2768 


1 


680 


0.7059 


15.56 


7.66 


129.3 



The overshoot layer appears especially attractive because it is the only place that can accom- 
modate a toroidal magnetic field of ~ 10^ G, suitable for the 11-year sunspot cycle. A similar 
field strength is inferred from the latitude of field emergence at the Sun's surface, as well as 
from the tilt angle of bipolar spot groups with respect to the east-west direction ||8|. With 
10^ G, and a layer thickness of, say, 5000 km, a neutrino should have a magnetic moment of 
~ 10"^'' //Bohr in order to suffer a noticeable spin precession. This would be in conflict with 
supernova SN1987A observations which indicate an upper limit of 10~^^//Bohr for the mag- 
netic moment of the neutrino. The reasoning may be even stronger as only a fraction of the 
overshoot layer normally will be filled with magnetic flux. Matter-enhanced (resonant) spin 
flip should also not occur since it requires an electron density exceeding that at the base of 
the convection zone. Therefore, there appears to be no reason to speculate about a solar-cycle 
dependence of the neutrino flux from the Sun. 
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The Garching solar model code is designed to calculate high precision solar models. It allows 
to control the numerical accuracy and has the best available input physics implemented [l|. 
It uses the OPAL-equation of state and for the opacities those of complemented by |0| 
in the low-temperature regions. Pre-main sequence evolution is also taken into account. The 
microscopic diffusion of hydrogen, helium, the isotopes participating in the CNO-cycle and 
some additional metals (Ne, Mg, Si) is incorporated following the description of for the 
diffusion constants. The nuclear reaction rates were taken from Q. 

In this work we want to emphasize the sensitivity of the structure of a solar model on the 
interpolation technique used for the opacity tables. As the run of opacities (k) with temper- 
ature (T) and/or density (p) may show very rapidly changing gradients, choosing a suitable 
interpolation procedure is not trivial. Our program uses two-dimensional bi-rational cubic 
splines to interpolate in the T-i?-grid of the opacity tables {R = p/T^, = T/IO^K). 
Apart from the general problem to choose suitable outer boundary conditions for the calcu- 
lation of the spline functions one can introduce an additional parameter to avoid artificial 
unphysical oscillations, as they are typical in the case of cubic splines and rapidly varying 
slopes. Increasing this damping parameter leads to almost linear interpolation just between 
two grid points and to a very rapidly changing gradient at the grid points themselves. The 
higher the damping parameter the more extended gets the linear region (a value of corre- 
sponds to standard cubic splines) . The disadvantage of this damping is that in regions where 



7' = 2.0 X 10'' K 




-2.0 -1.5 
logio ^ 



2.0 -1.5 
logio 



Figure 1: The run of k and dlog^/dlogi? with R at fixed temperature of 2 million Kelvin 
using pure cubic splines (dashed line) and damped cubic splines (full line). The arrows mark 
the value of R at this temperature in a typical solar model. 
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the opacity has a very slowly changing gradient too high values for the damping parameter 
lead to an interpolated run of opacity with T and/or R which shows step-like gradient vari- 
ations. This is illustrated in Fig. || where the full line shows the interpolated k{R) using a 
high damping parameter and the dashed line the run of k without damping. In this case 
a simple cubic spline seems to describe a better fit. The greatest deviations of k between 
two tabulated values is only about 3%, which is smaller than the quoted uncertainty for the 
opacities of about 10%. We would like to note here that neither of the chosen interpolation 
schemes can claim to reproduce the true values, as interpolation is always an estimation. 

To illustrate the influence of the opacity interpolation, two solar models were calculated, 
GARS0M3 with no and GARS0M2 with strong spline damping. The run of sound speed of 
both models is compared with a seismic mode from [^] in Fig. |2|. The deviation of GARS0M2 
(dashed line) from the seismic model just below the convective zone is about twice as large 
as compared with GARS0M3 (full line). At r = 0.67 Rq the temperature is approximately 
2 million Kelvin, logioi? ~ —1.7. Regarding Fig. || one notices that at this radius the model 
is in T-ii-regions of the opacity tables where the different interpolations schemes differ most. 



0.015 r 



quadratic sound speed 



0.000 - 




Figure 2: The run of quadratic sound speed of various models compared to seismic model 
from [7]. Full line: GARS0M2, dashed line: GARS0M3, dotted line: reference model from 
which the seismic model is inferred. 

As the possible error in the tabulated opacities (10%) is bigger than the difference due to 
interpolation (3%), neither of the models can be favoured from the theoretical point of view. 
Although we must admit that k{R, T) inferred from cubic splines without damping does look 
more reliable, we can not really rule out GARS0M2. It is therefore necessary to improve the 
grid density in the tables and the input physics in the opacities. 

To show that GARS0M3 is compatible with solar models from other groups, we have 
plotted in Fig. ^ also the run of sound speed of the reference model which was used to 
infer the seismic model of [^. The remaining discrepancies between the reference model and 
GARS0M3 may be due to slightly different ages or nuclear reaction rates. 
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With the successful completion of Gallex after six years of operation and the smooth tran- 
sition to Gno a milestone in radiochemical solar neutrino recording has been reached. The 
results from Gallex, 77±6±5SNU, and Sage, 74 lySNU, both being significantly be- 
low all solar model predictions, confirm the long standing solar neutrino puzzle and constitute 
an indication for non-standard neutrino properties. This conclusion is validated by the results 
of ^^Cr neutrino source experiments which have been performed by both collaborations and 
"^^As doping tests done in Gallex. 

GALLEX and SAGE: Radiochemical Solar Neutrino Recording 

The radiochemical gallium detectors, Gallex and Sage, have been measuring the integral 
solar neutrino flux exploiting the capture reaction ^^Ga -|- fe — > ^^Ge -|- e~ . The energy 
threshold being only 233 keV, this reaction allows to detect the pp-neutrinos from the initial 
solar fusion step which contribute about 90 % to the total solar neutrino flux. 

In a typical run the target, consisting of 30 tons of gallium in the form of 101 1 GaCls 
solution for Gallex and 55 1 of metallic gallium for Sage, is exposed to the solar neutrino 
flux for 3-4 weeks. In the following I will mainly focus on Gallex, as for Gallex and Sage 
the experimental procedure - apart from the chemical extraction of the neutrino produced 
''^Ge and the stable germanium carrier which is added at the beginning of each run - is 
rather similar. Both experiments use the signature provided by the Auger electrons and X- 
rays associated with the decay '^^Ge -|- e~ — > '^^Ga -|- for identification of ''^Ge during a 
several months counting time. Referring to |^, ^ for a detailed description of the detector 
setup and experimental procedure I concisely summarize the Gallex experimental program 
in table ||. 

Table 1: GALLEX experimental program. It comprises 4 periods of solar neutrino observations 
(Gallex 1 - Gallex 4), two chromium neutrino source experiments (Source I and Source II) and the 
arsenic test. 



date 


exposure period 


number of runs 


result 


5/91-4/92 
8/92-6/94 
6/94-10/94 
10/94-10/95 
10/95-9/96 
9/96-1/97 
1/97-3/97 


Gallex I 
Gallex II 
Source I (^^Cr) 

Gallex III 
Source II (^^Cr) 
Gallex IV 
As-test 


15 solar + 5 blank 
24 solar + 22 blank 

11 source runs 
14 solar + 4 blank 

7 source runs 
12 solar + 5 blank 

4 arsenic runs 


81±17±9SNU 
75 ± 10 +5 SNU 

54 ±11 ±3 SNU 
R = 0.84 toil 
118 ±19 ±8 SNU 
Y = 1.00 ±0.03 
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The combined result of all 65 Gallex solar runs is 76.4 ± 6.3 111 SNU. We note, however, 
that for the Gallex-IV period pulse shape information is used only for K-peak signals. The 
present result from Sage, 74 lySNU 0, is in perfect agreement with Gallex. 

The overall results from both experiments are only about 60% of the predictions from solar 
model calculations, which constitutes an indication for non-standard neutrino properties, even 
without considering the results from other solar neutrino experiments. 

The ^^Cr neutrino source experiments 

In order to validate this conclusion and check the reliability and efficiency of their detectors, 
both collaborations have prepared intense (ss 2MCi (Gallex) and 517kCi (Sage)) ^^Cr 
neutrino sources by neutron irradiating isotopically enriched Cr Q. ^^Cr decays by EC and 
emits neutrinos of energy 750 keV (decay to g.s., 90%) and 430 keV (decay to 320 keV excited 
level, 10%), nicely accommodating the neutrino energies of the solar ^Be branch. 

For Gallex, the source has been inserted in the central tube of the target tank. Alto- 
gether 18 extractions have been performed with the source in place, divided into two series of 
measurements with the source being re-activated in between. Figure |l| shows the individual 
run results. 



O Source 1: R = l.OliO.ll 
O Source 2: R = 0.84±0.12 
Combined: R = 0.93±0.08 




20 40 60 80 100 120 140 
TIME SINCE END OE BOMBARDMENT [d] 

Figure 1: Individual run results of the GALLEX chromium neutrino source experiments, 
normalized to the known source activity, as a function of time after the end of neutron 
irradiation of the source. The horizontal bars indicate the durations of the exposures, the 
exponential curve the expectation from the decay of the source. Source I and Source II runs 
are represented by circles and diamonds, respectively. 



An analysis of Source I and Source II yields R = 1.01 and R = 0.84 respectively, 
where R is the source activity as deduced from the neutrino measurement normalized to the 
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known true source activity. The combined analysis of both series results in R = 0.93 it 0.08, 
clearly demonstrating the absence of large systematic errors which could account for the 
observed 40% solar neutrino deficit. 

A similar experiment has also been performed by the Sage collaboration. They got a 
quantitative recovery ofi? = 0.95±0.12 

The fact that both experiments, though employing different chemistry, demonstrated in 
these experiments full efficiency, proves the trustworthiness of the radiochemical method and 
shows that the 40% solar neutrino deficit observed in the radiochemical gallium experiments 
cannot be ascribed to unknown systematic errors. In particular, the ^^Cr neutrino energies 
nicely accommodating those from the solar ^Be-branch, the full efficiency of the gallium 
experiments to ''Be-neutrinos is demonstrated. 

^^As experiments 

However, though the ^^Cr source in Gallex did outperform the sun by more than a factor 
15 after insertion into the target tank, the experiments still are low statistics, involving only 
several dozens of neutrino produced ''^Ge atoms. Therefore, at the very end of Gallex the 
collaboration has performed a large-scale test of potential effects of hot chemistry, which 
might lead to a different chemical behaviour of ^^Ge produced in a nuclear reaction compared 
to the stable Ge carrier isotope: The in-situ production of ''^Ge by /?— decay of ''^As. A known 
quantity of ^^As (O(IO^) atoms) has been added to the tank (t-sample), where it decayed with 
Ti/2 = 2.9 d to ^^Ge. Four runs have been made under different operating conditions (mixing, 
carrier addition, standing time), cf. table ^ For every spike a reference sample (e-sample) 
was kept aside, making possible to calculate the ratio of t- and e-sample which does not suffer 
from most of the systematic uncertainties associated with ^^Ge-counting. 



Table 2: Experimental conditions and results (ratio t-sample/e-sample) of the ''^As runs. 



run 


mixing conditions 
[h X m^/h] 


Ge-carrier 
addition 


standing 
time 


result 
(tank / external) 


Al 


22 X 5.5 


with As 


19.9 d 


1.01 ±0.03 




+0.17 X 60 








A2 


6 X 5.5 


no Gc-carrier 


19.9 d 


1.00 ±0.03 


B3-1 


24 X 5.5 


after As 


2.0d 


1.01 ±0.03 


B3-2 




after As 


22.0d 


1.00 ±0.03 



In all cases a quantitative recovery of 100% was achieved. This demonstrates on a 3%-level the 
absence of withholding effects even under unfavourable conditions like carrier-free operation]^. 

Gallium Neutrino Observatory 

With the ''^As-tests Gallex has completed its large-scale experimental program. However, 
solar neutrino measurements with a gallium target at Gran Sasso will be re-commenced in 
spring 1998 in the frame of the Gallium Neutrino Observatory (GNO) which is designed for 

^The stable germanium carrier is not only used for determination of the extraction yield, but also plays the 
role of an 'insurance' to saturate potential trace impurities which might capture '^^Ge in non- volatile complexes. 
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long-term operation covering at least one solar cycle. In its first phase GNO will use the 30 ton 
target of Gallex. However, for a second phase it is planned to increase the target mass to 
60 tons and later to 100 tons. In addition, as it is equally important to decrease the systematic 
uncertainty, effort is made to improve ''^Ge counting, both by improving the presently used 
proportional counters, and by investigating novel techniques like semiconductor devices and 
cryogenic detectors 
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Solar Neutrino Observation with Superkamiokande 

Yoshiyuki Fukuda (for the Superkamiokande Collaboration) 
Institute Cosmic Ray Research, University of Tokyo, Japan 

Introduction 

Superkamiokande, which is a second generation solar neutrino experiment, is an imaging 
water Cherenkov detector with 50,000 tons of pure water in the main tank. The detector is 
located 1000 m underground (2700 m water equivalent) in Kamioka Zinc mine in the Gifu 
prefecture of Japan, at 36.4°N, 137. 3°E and 25.8°N geomagnetic latitude. 




Figure 1: Schematic view of the Superkamiokande detector. 

The detector consists of a main inner counter and an outer anti-counter. A schematic view 
of the detector is shown in Fig. |l]. The inner counter is contained in a cylindrical stainless-steel 
tank and has a volume of 39.3 m in diameter x 42.0 m in height, containing 50,000 metric 
tons of water. A total 11,146 photo multipliers (PMTs) with 20 inch (j) photo cathode area 
cover 40% of the entire inner surface of the tank. The fiducial mass for the solar neutrino 
measurement is 22,500 tons, with boundaries 2.0 m from the inner surface. On the other 
hand, a 47r solid-angle anti-counter surrounding the inner counter is also a water Cerenkov 
counter of total mass 13,000 metric tons with 1850 PMTs to detect any signals coming from 
outside of the detector and to shield against gamma-rays and neutrons. 

Calibrations 

The timing calibration for all PMTs is done by the Xe lamp. We usually take those data for 
every 3 months and the maps of T-Q response for all PMTs are produced. This table is used 
for correct timing as a function of observed charge in real data analysis. 
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The variation of the water transparency is obtained by stopping muon spectrum (Michael 
spectrum). This is very important for obtaining the energy scale to be stable as a function 
of time. The energy is almost proportional to the number of hitted PMTs, however, the 
number maybe variable to the attenuation length of the water. Corrected number of hitted 
PMTs is very stable to the variation of water transparency. This is also confirmed by the 
7-ray source which is emitted by the reaction of Ni(n,7)Ni*. The absolute values of water 
transparency at various points of wavelength are measured by DYE-laser calibration. Recent 
water transparency itself is very stable. 

The performances of the detector, such as the absolute energy scale, the energy resolution, 
the vertex resolution and the angular resolution are mainly calibrated by LINAC system. The 
LINAC can generate an electron beam with 5 MeV to 16 MeV. This beam is induced via beam 
pipe and bent by magnetic coil into the several positions of the tank. The typical values for 
each resolution is 16%, 70cm and 22 degree for 10 MeV electrons. The absolute energy for 
electron beam is calibrated by Ge solid state detector at each calibration time. 

Monte Carlo simulation was tuned by the water transparency and other parameters to 
reproduce the energy scale in various tank positions within 1% difference between MC and 
real data. 



Solar neutrino analysis 

Superkamiokande has started from April 1996 and now processed about 400 days data for solar 
neutrinos measurement. Analysis procedures are (1) noise reduction, (2) vertex reconstruction 
and tight noise reduction, (3) spallation products cut, (4) fiducial volume (22.5 kton) cut and 
(5) gamma ray cut. Some of analysis techniques are similar to Kamiokande's ones, however, 
most of them have been newly developed. In first step, we eliminate /i and decayed electron 
and several electronics noise. The second step reconstructs the vertex using the timing of 
hitted PMTs within selected 50 n second window. Main backgrounds in the residual data are 
spallation products. In order to eliminate those events, we used the likelihood method using 
the time difference and the distance between induced muon and those spallation events as 
a function of muon energy. We can reject most of spallation products with 20% dead time. 
In the last reduction, we reject external 7-rays coming from outside of the detector. Most 
of those events are sitting at very close to the edge of fiducial volume and have an opposite 
direction with respect to the detector wall. Event rate of the final sample is obtained by 175 
events per day per 22.5 kton for 6.5 to 20 MeV. 

From 1st May 1996 to 22 Oct 1997, we obtained 374.2 days data from Superkamiokande 
measurement [j^]. Figure ^ shows the angular distribution to the solar direction and the 
heliograph for obtained final sample. In Fig. |^(a), best fit line which is expected by MC is 
also shown. Extracted number of solar neutrinos is obtained by this fit as 4951^]^|J 3 for 374.2 
days data. Observed solar neutrino flux is given by 

(t>{^B) = 2.37l°:°^(stat.)10:°?(syst.) x lO^cm'^sec-i 

or by taking ratio to the BP95 flux 

Data 



0.358l°:°°^(stat.)l°:°}^(syst.). 



SSMbP95 

The systematic error related to the energy scale has been reduced by LINAC calibration 
and is obtained as 2.3%. Observed neutrino flux is significantly deficit to the expectation 
of SSM(BP95) and it is consistent with the result from Kamiokande 0]. 
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Figure 2: (a) Angular distribution to the solar direction and (b) heliograph for 374.2 days 
data. 



The energy spectrum of observed neutrinos is shown in Fig. |3|(a). In this figure, expected 
spectra of MSW small angle parameter and just-so parameter are also shown. At first sight, 
small angle solution has better fit than flat (no oscillation) and just-so solution, but it is not 
significant within an experimental error. The day and night flux difference is obtained by; 

D-N , , , , 

= -0.031 ± 0.024(stat.) ± 0.014(syst.). 

If night data are divided into five bins, those differences are shown in Fig. ^(b). In this figure, 
typical day/night flux variation of the large angle and the small angle solutions are also shown. 
There is no significant difference in day/night fluxes in present observation. Also Fig. ^c) 
shows the seasonal variation of solar neutrino fluxes. Each season is pile up among different 
years. Solid line corresponds to the expected variation from an eccentricity of the Sun orbit. 
Within experimental error, there is no seasonal variation in present analysis. Those results 
are also same ones from Kamiokande. 



Two flavor neutrino oscillation 

For astrophysical solution, it is generally difficult to explain the solar neutrino problem with 
the modification of SSM including the observations from helioseismology. On the other hands, 
the elementary particle solution using MSW neutrino oscillation ||5| seems to be an excellent 
for explanation of the solar neutrino problem, because it can distort the spectra of solar 
neutrinos. Also MSW oscillation can give a effect in the day/night fluxes variation. As 
obtained by Fig. our observed spectrum can be seen slightly as distorted one, but not 
seen in variance between day/night fluxes. Using these results, we can obtained the excluded 
region at 95% C.L. in MS diagram as shown in Fig. |^(a). 

If we take the constraint of measured solar neutrino flux, allowed region of 95%, 90% 
and 68% C.L. as shown in Fig. §(b). From these figures, our results and the allowed region 
given by Ref. Q are consistent in MSW oscillation analysis. 
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Figure 3: (a) Energy spectrum, (b) Day/Night flux variation and (c) seasonal variation of 
obtained solar neutrino signal. 



Summary 

In summary, solar neutrino observation in Superkamiokande has started since May 1996 and 
has already taken 374.2 days data. The observed solar neutrino flux is about 37% of the 
prediction from SSM(BP95) and it is almost consistent with the result from Kamiokande. 
Using LIN AC calibration system, we can reduce the systematic errors related to the energy 
scale. Obtained energy spectrum is likely distorted and it is indicated that the new physical 
solution will be solved the solar neutrino problem. New challenge to lower threshold analysis 
(> 5 MeV) has been started. Present radioactive (^^^Rn) level is about 3 mBq/m"^, however, 
we will able to reduce that level to factor 1/5. Even though the present analysis, we have 
succeeded to extract solar neutrino signals from 5-6.5 MeV region. Those data will give a 
strong indication to the solution for the solar neutrino problem within a few years. 
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Figure 4: (a) 95% C.L. excluded region in MSW parameters. Inside of solid and dotted line 
are excluded by the observed energy spectrum and the day/night flux difference, respectively, 
(b) Allowed region in MSW parameters with measured ®B neutrino flux. Thick solid line, 
dashed line and thin solid line correspond to 95%, 90% and 68% C.L., respectively. Shaded 
region is allowed at 95% C.L. by Hata. 
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The Sudbury Neutrino Observatory 

M.E. Moorhead (on behalf of the SNO Collaboration) 

Particle and Nuclear Physics Laboratory, Keble Road, Oxford 0X1 3RH, UK 

The Sudbury Neutrino Observatory (SNO) Q is a 1,000 ton heavy water (D2O) Cherenkov 
detector in Sudbury, Ontario (Canada) which will start taking data in '98. The u reactions 
which occur in D2O and the extremely low background environment of the detector will allow 
the following measurements: i) the Ve and i/^ (flavour independent) fluxes, and their ratio, 
for '^B solar z/'s, ii) the energy spectrum of i/g's above 5 MeV, iii) time dependence in the 
^B flux, and iv) detailed studies of the v burst from a galactic supernova, including a search 
for and Uj- masses above 20 eV. For the ^B solar u measurements, the Ve/^^x flux ratio 
and the Ve energy spectrum constitute two separate tests of v oscillations which are both 
independent of solar model flux calculations ||2| . If the currently favoured MSW solution |^ 
of the solar neutrino problem is correct then the Ve/^x flux ratio should provide conclusive 
proof of oscillations with one year's data yielding a 17 sigma departure from unity (statistical 
error only). 

The detector is situated two kilometers underground in a dedicated laboratory that has 
been excavated in the Creighton nickel mine of INCO Corporation. This laboratory com- 
prises facilities for changing into clean-room clothing, a lunch room, a car wash for bringing 
equipment into the clean area, a utility room where the H2O and D2O systems are located, 
a control room and a 30 m x 23 m barrel shaped cavity for the detector itself. The walls 
of the cavity have been coated with concrete and low-activity Urylon, a water proof radon 
barrier. Inside this cavity is located a 12 m diameter spherical acrylic vessel (AV), recently 
completed, for containing the D2O neutrino target. Almost completely surrounding this AV, 
stands a 17 m geodesic sphere supporting 9,500 20-cm-Hamamatsu PMTs, each of which is 
equipped with a light reflecting concentrator to increase its effective photocathode area by a 
factor 1.7. Beginning in March '98, the detector will be filled with 7,000 tons of high purity 
H2O outside the acrylic vessel (to act as shielding for high energy gamma rays coming from 
the rock and the PMTs) and 1,000 tons of D2O inside the AV. The filling will take between 
3-4 months and hence data taking will begin in mid-98. 

The event rates for solar z/'s, assuming the full SSM ^B fiux and for a supernova 
(SN) at the center of our galaxy are given in Table 1. Apart from the neutral current (NC) 
reaction, all of the events are detected by the array of 9,500 PMTs via the Cherenkov radiation 
emitted by a single electron (or positron) of energy > 5 MeV, the detector's threshold. The 
NC reaction produces a free neutron in the D2O which can be detected, after thermalization, 
by observing a subsequent neutron capture reaction. There are three capture reactions of 
interest depending on what additives are placed in the D2O: 

i) Pure D2O: In the case of no additive there is a 30% probability of capture on deuterium, 
which produces a 6.25 MeV gamma. This gamma converts to electrons by Compton scattering 
and pair production, and the resulting Cherenkov light is detected by the PMTs. Five hundred 
events a year are expected above the detector's 5 MeV threshold. 

ii) MgCl2: Dissolving 2 tons of MgCl2 in the D2O gives an 83% chance of neutron capture 
on ^^Cl which produces an 8.5 MeV gamma cascade. The higher efficiency and Q-value of this 
capture (c.f. capture on deuterium in the pure D2O case) increases the number of detected 
events by a factor of 5 to 2,500 per year. 
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Table 1: Neutrino event rates in SNO (including detection efficiency). SSM refers to the event 
rate per year for solar neutrinos assuming the standard solar model [Q] and the small angle 
MSW solution j^. SN refers to a type II supernova at a distance of 10 kpc (the center of the 
galaxy). 



Neutrino reaction 






SSM 


SN 


Charged Current (CC): 




■ P + P + 


3000 


80 


Neutral Current (NC): 


v^ + d — 


> p + n + z^x 


2500 


300 


Electron Scattering (ES): 


Ve,x + e" 


Ve,x + e" 


400 


20 


Anti- neutrino CC in D2O: 


i^e + — * 


• n + n + 





70 


Anti-neutrino CC in H2O: 


Pe +?> ^ 


■ n + e"*" 





350 



iii) ^He Counters: An array of ^He proportional counters |^ (5 cm diameter tubes of 
800 m total length) placed vertically in the D2O in a square grid of Im spacing, gives a 42% 
chance of neutron capture on ^He. The energy and rise-time of the signals are used to separate 
n-capture (2,000 per year) from internal alpha and beta backgrounds. 

The dominant background for all these NC detection methods will probably come from 
photodisintegration of deuterium which produces free neutrons that are indistinguishable 
from NC neutrons. Thus, the detector components have been carefully selected for extremely 
low levels of thorium and uranium chain contamination so that the photodisintegration rate 
is small compared with the NC rate. This small residual photodisintegration rate must 
be measured, in order to subtract its contribution to the neutron capture signal. Several 
methods have been developed for this purpose: i) radiochemical extraction and counting of 
^^^Th, ^^^Ra, ^^^Ra, ^^^Rn and ^^^Pb, ii) analysis of low energy signals seen by the PMT 
array, iii) delayed coincidences between signals seen by the PMT array, and iv) prompt and 
delayed coincidences between signals seen by the PMTs and signals in the ^He proportional 
counters. 
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BOREXINO 

L. Oberauer (for the Borexino Collaboration) 

Technische Universitdt Miinchen, Physik Department E15, 8574.7 Garching, Germany 
and Sonderforschungsbereich 375 Teilchen-Astrophysik 

Physics Goals and Neutrino Detection with Borexino 

The aim of the solar neutrino experiment Borexino is to measure in real time the solar 
neutrino flux with low energy threshold at high statistics, and energy resolving via pure 
leptonic neutrino electron scattering v + e ^ v + e. 

Motivation for Borexino comes from the long standing solar neutrino puzzle. Data analysis 
of the existing experiments leads to the assumption of severe suppression of the solar ^Be- 
branch, which probably cannot be explained by modifications of the standard astrophysical 
model of the sun. The monoenergetic ^Be-neutrinos give rise to a compton like recoil spectrum 
in Borexino. Its edge will be at 660 keV, significantly higher than the aimed energy threshold 
of ca. 250 keV. Thus ^Be-neutrino detection is very efficient in Borexino. 

Assuming validity of the standard model a counting rate for '''Be-neutrinos, which would 
consist in this case purely as Vg, of roughly 55/day in Borexino is expected. 

In scenarios of total neutrino flavour conversion (i.e. for neutrino mass differences Am^ 
10~^— 10~^ eV^) a reduced flux of approximately 12/day would be measured due to the lower 
cross section of Uf^^j- scattering, which occurs only via neutral current interaction. 

In case of vacuum oscillations (i.e. for neutrino mass differences Am^ « lO"^'^ eV^) Borex- 
ino would see a distinct time dependent periodical neutrino signal due to the seasonal eccen- 
tricity of the earth's orbit around the sun. 

For neutrino mass differences in the range of Am^ ~ 10"'^ cV^ and for large mixing 
Borexino should see a 'day/night' effect due to electron neutrino regeneration during the 
path through the earth. 

Borexino also can serve for additional projects in neutrino physics. Search for a magnetic 
moment can be performed by means of terrestrial neutrino sources by investigating the elec- 
tron recoil shape at low momentum transfer. Via the inverse beta-decay i?e + p ^ e+ + n 
Borexino can look for signals from geophysical neutrinos as well as for neutrinos emitted by 
european nuclear power plants. The latter would serve as a long baseline neutrino oscillation 
experiment probing the so-called large mixing angle solution for the solar neutrino problem. 

The Detector and Background Considerations 

The detector is shielded successively from outer radioactivity by means of an onion-like struc- 
ture. Here the adjacent inner layer serves as shielding and has to provide an increased purity 
in terms of internal radioactivity. 

Borexino is placed in hall C of the underground laboratory at Gran Sasso, Italy. An 
overburden of ca. 3500 m.w.e. suppresses the cosmic muon flux to roughly 1/m^. The outer 
part of the detector consists of a steel tank with 18 m in height and diameter. Inside this 
'external' tank a stainless steel sphere will support 2200 phototubes on the inside and 200 
tubes at the outside. Most of the tubes inside the sphere will be equipped with light guides in 
order to increase the geometrical coverage and hence the energy resolution. Between external 
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tank and sphere high purity water will serve as shielding against external gamma rays and 
as active Cherenkov counter against cosmic muons. The steel sphere will be filled with a 
tranparent, high purity buffer liquid which itself holds a nylon sphere, filled with organic 
scintillator. The active scintillator mass will be around 300 t. By means of time of flight 
measurements event position can be reconstructed and a fiducial volume for solar neutrino 
interaction defined. The latter should be about 100 t, establishing a counting rate of 55 
neutrinos per day according to the standard solar model. The outer part of the scintillator 
sphere serves as active shielding. 

The demands on purity in terms of radioactivity in Borexino, especially for the scintillator 
itself, are very severe. In order to be able to extract a clear signal from background events 
also in case of total flavour conversion, an intrinsic concentration in Uranium and Thorium 
of ca. 10~^^ should not be exceeded significantly. The amount of ^^C/^^C must not be higher 
than ~ 10"^^. In order to test scintillating materials a large Counting Test Facility (CTF) 
has been built up in hall C of the underground laboratory at Gran Sasso, which resembles to 
a small prototyp (ca. 5 t of scintillator) of Borexino. From beginning of 1995 until summer 
1997 several tests about the feasibilty of Borexino including procedures to maintain the purity 
of the scintillator has been performed, which showed very encouraging results: ^^C/^^C = 
1.85•10-l^ 238u = (3.5±1.3)-10-i6 g/g, 232Th = (4.4±1.5)-10-^6 g/g. A complete discussion 
of the CTF results including experimental techniques for further background suppression is 
given in Q and [|[. Details about the experimental setup of the CTF can be found in Q. 

Highly developed neutron activation analysis of scintillation samples performed in Munich 
is now sensitive in the same regime. For uranium an upper limit of ^^^U < 2 • 10~^^ g/g (90% 
CL) has been obtained. In addition concentration values or limits have been measured by 
this method for a various amount of isotopes, including man-made nuclei. For details, see Q. 

Background studies for Borexino include also the interaction of cosmic muons. The direct 
detector response on muons has been determined by a coincidence measurement between the 
CTF and a muon telescope on top of it. The time distribution of such events can be used 
to discriminate between muon events and neutrino candidates at a level of 98%. However, 
to reach the sensitivity needed to reach the goals in Borexino, the leak rate for muons must 
not exceed a level of ~ 10~^. Our design of the muon veto system therefore is threefold: 
The outer region between external tank and steel sphere acts as Cherenkov counter, a special 
configuration of the tubes inside the sphere will act as an additional muon identification 
system, and finally the offline study of event topology like the time structure will help to 
suppress this kind of background sufficiently. 

Cosmogenic generation of radioactive nuclei has been sudied this fall at the 180 GeV muon 
beam at SPS in CERN. Most dangerous source of events will come from ^^C-production in 
the scintillator and surrounding buffer liquid. However, the energy spectrum of these events 
is between 1 MeV and 2 MeV since the decay mode is positron decay at 1 MeV endpoint 
energy. Thus the detection of ^Be-neutrinos is not affected, however that of pep-neutrinos. 

Prospects 

Borexino is an international collaboration of about 60 scientists. Approved funding already 
comes from INFN (Italy) and from BMBF and DFG (Germany). A substantial part should 
also be covered by NSF (USA) in the near future. Work on the external tank of Borexino is 
almost completed. We expect to finish with the inner steel sphere in 1999. Simultaneously 
the CTF will be upgraded. Finally it will serve as test facility for Borexino scintillator 
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procurement in batch mode. Given full funding also for our american collaborators first data 
taking may be expected at the end of the year 2000. 
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Measurements of Low Energy Nuclear Cross Sections 

M. Junker (for the LUNA Collaboration) 

Laboratori Nazionali Gran Sasso, Assergi (AQ), Italy 

and Institut fiir Experimentalphysik III, Ruhr- Universitdt Bochum, Germany 

The nuclear reactions of the pp-chain play a key role in the understanding of energy pro- 
duction, nucleosynthesis and neutrino emission of the elements in stars and especially in our 
sun Q. A comparison of the observed solar neutrino fluxes measured by the experiments 
GALLEX/SAGE, HOMESTAKE and KAMIOKANDE provides to date no unique picture 
of the microscopic processes in the sun [Q] . A solution of this so called "solar neutrino puz- 
zle" can possibly be found in the areas of neutrino physics, solar physics (models) or nuclear 
physics. In view of the important conclusions on non-standard physics, which might be de- 
rived from the results of the present and future solar neutrino experiments, it is essential to 
determine the neutrino source power of the sun more reliably. 

Due to the Coulomb Barrier involved in the nuclear fusion reactions, the cross sec- 
tion of a nuclear reaction drops nearly exponentially at energies which are lower than the 
Coulomb Barrier, leading to a low-energy limit of the feasible cross section measurements 
in a laboratory at the earth surface. Since this energy limit is far above the thermal en- 
ergy region of the sun, the high energy data have to be extrapolated down to the en- 
ergy region of interest transforming the exponentially dropping cross section (7{E) to the 
astrophysical S'-factor S{E) = a{E) E exp{2 7rr]), with the Sommerfeld parameter given by 
2 7rry = 31.29 Zi Z2{fi/E)^^'^ |||]. The quantities Zi and Z2 are the nuclear charges of the in- 
teracting particles in the entrance channel, ^ is the reduced mass (in units of amu), and E 
is the center-of-mass energy (in units of keV). In case of a non resonant reaction S{E) may 
then be parameterized by the polynomial S{E) = S{0) + S'{0)E + 0.5S"{0)E'^. 

As usual in physics, extrapolation of data into the "unknown" can lead onto "icy ground". 
Although experimental techniques have improved over the years to extend cross section mea- 
surements to lower energies, it has not yet been possible to perform measurements within the 
thermal energy region in stars. 

For nuclear reactions studied in the laboratory, the target nuclei and the projectiles are 
usually in the form of neutral atoms/molecules and ions, respectively. The electron clouds 
surrounding the interacting nuclides act as a screening potential: the projectile effectively sees 
a reduced Coulomb barrier. This leads to a higher cross section, as{E), than would be the 
case for bare nuclei, a]^{E), with an exponential enhancement factor fiahiE) = crs{E)/ah{E) ~ 
exp(7rr/ Ue/E) [^, ^ where Ue is the electron-screening potential energy (e.g. Uc — Zi-Z2-e'^ / R^, 
approximately, with i?a an atomic radius). For a stellar plasma the value of (Jh{E) must be 
known because the screening in the plasma can be quite different from that in laboratory stud- 
ies and (Tb(£^) must be explicitly included in each situation. Thus, a good understanding 
of electron-screening effects is needed to arrive at reliable iy\y{E) data at low energies. Low- 
energy studies of several fusion reactions involving light nuclides showed ||6|, |^, § indeed the 
exponential enhancement of the cross section at low energies. The observed enhancement 
(i.e. the value of C/g) was in all cases close to or higher than the adiabatic limit derived from 
atomic-physics models. An exception are the ^He-l-^He data of Krauss et al. [jl^, which 
show apparently no electron screening down to E=25 keV, although the effects of electron 
screening should have enhanced the data at 25 keV by about a factor 1.2 for the adiabatic 
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limit C/c=240 eV. Thus, improved low-energy data are particularly desirable for this reac- 
tion. 

The low-energy studies of thermonuclear reactions in a laboratory at the earth's surface 
are hampered predominantly by the effects of cosmic rays in the detectors. Passive shielding 
around the detectors provides a reduction of gammas and neutrons from the environment, 
but it produces at the same time an increase of gammas and neutrons due to the cosmic- 
ray interactions in the shielding itself. A 47r active shielding can only partially reduce the 
problem of cosmic-ray activation. The best solution is to install an accelerator facility |Tl| ] in 
a laboratory deep underground [^. The worldwide first underground accelerator facility has 
been installed at the Laboratori Nazionali del Gran Sasso (LNGS) in Italy, based on a 50 kV 
accelerator. This pilot project is called LUNA and has been supported since 1992 by INFN, 
BMBF, DAAD-VIGONI and NSF/NATO. 

The major aim of the LUNA project is to measure the cross section of ^He(^He,2p)^He 
which is one of the major sources of uncertainties for the calculation of the neutrino source 
power. It has been studied previously [10| down to about -Ecm=25 keV, but there remains the 



possibility of a narrow resonance at lower energies that could enhance the rate of path I at 
the expense of the alternative paths of the pp-chain that produce the high-energies neutrinos 
{El, > 0.8 MeV). The LUNA-facility allows to study this important reaction over the full 
range of the solar Gamow Peak, where the cross section is as low as 8 pbarn at £'cm=25 keV 
and about 20 fbarn at Scm=17 keV. 

Figure || shows the results obtained in the energy region between £'cm=25 keV and 
£'cm=20.7 keV. The lowest counting rate was of 3 events per day at ii^cm=20.7 keV. At 
this energy about 1000 Cb of ^He"*" have been accumulated on the target. The data obtained 



at higher energies (450 kV accelerator in Bochum) with the LUNA setup [11| are included 
for completeness. Previous data obtained ||T^ at £'cm=25 are also shown in figure |l|. The 
LUNA data have been obtained at energies within the solar Gamow Peak, i.e. below the 21 
keV center of this peak, and represent the first measurement of an important fusion cross 
section in the thermal energy region. 



The astrophysical 5- factor has been extrapolated to zero energy |12]: 



S{0) = 5.40 ± 0.05(stat) ± 0.30(sys) ± 0.30(screen). 



The first error contains the statistical error (one standard deviation) including counting statis- 
tics and apparative variations of pressure, beam power and temperature measurement. The 
second error is the systematical error (one standard deviation) including uncertainties in pres- 
sure, beam energy and power, efficiency and energy loss data. The 10% error of the energy 
loss data transforms to an uncertainty of 0.1 to 0.2 keV in the effective energy. This in turn 
leads to an error of 1 to 3.5 % for the astrophysical S-factor. This is due to the exponential 
decrease of the cross section at low energies jl^. The third component of the error results 
from the lack of understanding of electron screening: Fitting all existing data sets with S{0), 
S"(0), S"{0) and Uc as free parameter gives an Electron Screening Potential Ue of 323 eV with 
5(0) = 5.30 MeVb. Fixing 5(0), 5'(0) and 5"(0) at energies higher than 100 keV (where the 
electron screening effect is neglectible) gives S{0) = 5.1 MeVb. In turn fitting Ue in a second 
step including also the low energy data gives Ue = 432 eV. A fit to the data assuming no 
electron screening gives S{0) = 5.7 MeVb. The adiabatic limit for Ue based on theoretical 
atomic physics is 240 eV. The first two methods give Ug values higher than the adiabatic limit 
(240 eV), consistent with observations in other fusion reactions. The difference between ob- 
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Figure 1: The S{E) factor of ^He^He from [p^ j and the present work. The dashed and 
soHd curves represent S\,{E) and Ss{E) respectively. The solar Gamow Peak is shown in 
arbitrary units. 



served and predicted Ue values is not understood at present giving a considerable uncertainty 
on 5(0). 

Based on the success of the LUNA project and on the fruitful international collaboration 
(in particular the support provided by LNGS), further experiments with the 50 kV machine for 
a better understanding of the electron screening effect are planned. The proposed reactions are 
D(p,7)^He, D('^He,p)^He, ^Li(p,a)^He and ^^B(p,a)^Be(Q!)a. In addition, some interesting 
aspects in nucleosynthesis of the early universe are related to the low energy cross section of 
the first reaction |13|. 

While studying these reactions an enlarged LUNA II collaboration will install a 200 kV 
accelerator facility at Gran Sasso to study the pp-reactions ^He(a,7) ''Be and '^Be(p,7)^B and 
the key reaction of the CNO-cycles, ^^N(p,7)^^0, at energies far below the present limits. 
All these reactions are critical to the solar neutrino puzzle. The reaction rate of ^^N(p,7)^^0 
is also one of the ingredients needed to determine the theoretical scenario used to constrain 
both the age and the distance of the oldest stellar system in our galaxy, namely the Globular 
Clusters. In such a way, more stringent limits to the cosmology will be obtained. 

In addition, the new 200 kV accelerator will afford the study of many (p,7) reactions of 
the NeNa and MgAl cycles below an incident proton energy £'p=200 keV. Experimental data 
about these channels, very important for the understanding of nucleosynthesis processes, are 
up today still missing or very uncertain. For example the NeNa cycle may play a role in 
understanding the almost pure ^^Ne abundance found in meteorites samples, while the MgAl 
cycle may provide the mechanisms for production of ^^Al, the decay of which gives rise to 
the ■^^Mg/^''A1 anomaly found in some meteorites. All the involved (p,7) cross sections of 
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these cycles are scarcely known at low energies: a continuation of the underground experi- 
mental program of the LUNA experiment could be devoted in the future to the measurement 
of these channels for proton incident energies below 200 keV, where many unmeasured reso- 



nances are present |14]. For example the strength of the low lying resonances in the reaction 



^^Mg(p,7)^^Al, which is crucial for the production of ^^Al, could be experimentally deter- 
mined for the first time. Also the other reactions of these cycles can be investigated using 
stable targets and accelerator and detector systems of the LUNA experiment. 
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Solar Neutrinos: Where We Are and What Is Next? 

G. Fiorentini 

Dipartimento di Fisica dell'Universitd di Ferrara and Istituto Nazionale di Fisica Nu- 
cleare, Sezione di Ferrara, 1-44100 Ferrara, Italy 



What has been Measured? 

All five experiments report a deficit of solar neutrinos with respect to the predictions of 
Standard Solar Models (SSMs), see Ref. |]. 

What Have We Learnt on Solar Neutrinos, Independently of SSMs? 

Actually, the solar neutrino puzzle (SNP) is not just the discrepancy between experimental 
results and the predictions of standard solar models. Rather, experimental results look incon- 
sistent among each other with the only assumption that the present total neutrino flux can 
be deduced from the present solar luminosity, unless something happens to neutrinos during 
the trip from Sun to Earth, see Fig. I and Refs. @, |, i. 

What Has Been Calculated? 

Accurate predictions of solar neutrino fluxes are anyhow extremely important. If neutrino 
masses (differences) are as small as suggested by several proposed solutions to the SNP, then 
the only way to measure neutrino masses is through the interpretation of future solar neutrino 
experiments, which requires accurate theoretical predictions of solar properties. 

Refined solar models are thus necessary. All these have to be confronted with the pow- 
erful helioseismic constraints, see Ref. @, particularly for a quantitative (and conservative) 
determination of the accuracy of solar properties as deduced from helioseismology. 

Recent SSM calculations, using accurate equations of state, recent opacity calculations 
and including microscopic diffusion, look in agreement with heliosesimology, see Figs. I and I 
and Ref. Alternative solar models should be as successful as these are [^, ^. 

Actually, one can exploit helioseismology within a different strategy. One can relax some 
assumptions on the most controversial ingredients of solar models (e.g. opacity and metal 
abundance) and determine them by requiring that helioseismic constraints are satisifed. These 
helioseismically constrained solar models (HCSM) all yield the same central temperature 
within about one percent 0. The main uncertainties for the determination of solar neutrino 
fluxes arise now from nuclear physics measurements. After the succesful LUNA experiments 
at LNGS the main uncertainties are now from the ^He-|-^He and ^Be-|-p reactions. 

What Is Missing? 

In a prophetical paper of 1946 @ Bruno Pontecorvo wrote: "direct proof of the existence of 
the neutrino . . . must be based on experiments the interpretation of which does not require 
the law of conservation of energy, i.e. on experiments in which some characteristic process 
produced by free neutrinos ... is observed." 
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Table 1: The proposed solutions, their fingerprints and the experiments looking at them. 



Proposed solutions Signatures 







Oscillation 


spectral 


Day-night 


Seasonal 


CC/NC 






at reactor 


deformation 


variation 


modulation 


events 


MSW small an| 


^le 


NO 


TINY 


TINY 


NO 


YES 


MSW large ang 


jle 


NO 


TINY 


YES 


NO 


YES 


JUST-SO 




NO 


YES 


NO 


YES 


YES 


Universal oscil. 




YES 


NO 


NO 


NO 


YES 


Experiment 




CHOOZ 


SUPERKAM. 


SUPERKAM. 


BOREXINO 


SNO 


Data 




now 


now 


now 


2000 


2000 



The situation now looks very similar, just change existence with nonstandard properties, in 
that the strongest argument for a particle physics solution to the SNP arises from energy con- 
servation (the luminosity constraint) and actually we need a direct footprint of some neutrino 
property, not predicted within the minimal standard model of electroweak interactions. 

The four most popular particle physics solutions (small and large angle MSW effect, just 
so oscillations and universal oscillations) all predict specific signatures which are being or 
will be tested by the new generation of experiments (Superkamiokande, Borexino, SNO, . . .), 
see Table |. 

The hypothesis of universal oscillation proposed in Ref. llO] has just been falsified by the 



recent negative result of CHOOZ |11]. This nice and small (in comparison with the gigantic 
solar neutrino devices) experiment at a nuclear reactor is cleaning some of the fog in the air. 

Let us wait and wish that (at least) one of the fingerprints of neutrino oscillations is 
detected by the new experiments. 
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Figure 1: The and ^Be+CNO neutrino fluxes, consistent with the luminosity constraint 
and experimental results for standard neutrinos. The dashed (solid) lines correspond to the 
central (±lcr) experimental values for chlorine |13|, gallium |jl^ and v-e scattering exper- 
iments p5[ . The dashed area corresponds to the region within 2a from each experimental 
result. The predictions of solar models including element diffusion (full circles) ||l^, |l^, 17, |l^ 
and neglecting diffusion (open circles) [ |19| are also shown. The dotted lines indicate the be- 
haviour of nonstandard solar models with low central temperature O. 



44 




0,1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

R/Ro 



Figure 2: The difference between U as predicted by selected solar models, f/mod, and the 
helioseismic determination, Uq, normalized to this latter (thin line). The allowed area is that 
within the thick lines. BP95 is the model with metal and He diffusion of Ref. JCD96 is 
the "model S" of Ref. pO|; D96 is the "model 0" of Ref. pT|; CDR97 is the "best" model with 



He and heavier elements diffusion of Ref. |17] 



o 





Figure 3: Helioseismic determinations and solar model predictions of properties of the con- 
vective envelope. The box defines the region allowed by helioseismology. Open circles denote 
models without diffusion, squares models with He diffusion, full circles models with He and 
heavier elements diffusion Q. 
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Phenomenology of Supernova Explosions 

Wolfgang Hillebrandt 

Max- Planck- Institut fiir Astrophysik, D-85748 Garching, Germany 
Abstract 

An attempt is made to match observed properties of supernovae, such as spectra, hghtcurves 
etc. with theoretical ideas and predictions. It is demonstrated that neither do observational 
data constrain the models in a unique way nor does theory allow for an unambiguous inter- 
pretation of the data. Possible improvements on the theoretical side are briefly discussed. 

Some Observational Facts 

In general, supernovae are classified according to their maximum light spectra. Those showing 
Balmer lines of H are called Type II's, and all the others are of Type I. Those Type I's which 
show a strong Si absorption feature at maximum light are named Type la, and the others are 
Ib's or Ic's, depending on whether or not they have also He I features in their spectra [p. At 
later times, several months after the explosion, when the supernova ejecta become optically 
thin. Type II spectra are dominated by emission lines of H, O, and Ca, whereas Type la's 
have no O, but Fe and Co. Type Ib,c's, on the other hand side, show emission lines of O and 
Ca, just like the Type II's However, the spectral classification is not always as clear. For 
example, SN 1987K started out as a Type II with H lines, but changed into a Type Ib,c like 
spectrum after 6 months |||]. 

As far as the light curves are concerned, Type II's seem to be more complicated than 
Type I's. Type II-L are characterized by a peak lasting for about 100 days, followed by a 
"linear" decay in the blue magnitude vs. time diagram. In contrast. Type II-P's have a 
somewhat wider peak, followed by a "plateau" phase and an occasionally rather complicated 
tail decaying not like a single exponential. Typically, Type II-L's are brighter than Type II- 
P's in the blue. In addition, there are objects like SN 1987A which possess a very complicated 
light curve, with an early narrow peak, a first minimum, followed by a broad hump after a 
few months, and a final nearly exponential decay with indications of some flattening at late 
times. SN 1987A-like objects are much dimmer than all other Type II's [j^. 

In contrast, all Type la light curves are quite similar, making them good candidates for 
standard candles to measure the cosmic distance scale. Moreover, since they are the brightest 
among all supernovae, they can be observed even at high redshifts and, in fact, a Type la 
supernova at a redshift of about 1 has recently been observed Although their absolute 
peak luminosity may vary by about one magnitude, an observed correlation between the 
luminosity and light curve shape (the brighter ones have broader light curves) allows one 
to correct for the differences. So recently observations of Type la's have become a powerful 
tool in attempts to determine cosmological parameters Of course, one has to make sure 
that the supernovae one is observing are indeed of Type la, and not of Type Ib,c which are 
intrinsically fainter, but this can be achieved if a spectrum near maximum light is available. 

Other observational information on supernovae is, in general, less solid. With the excep- 
tion of SN 1987A, we do not have direct information on the properties of the progenitors, 
the energetics, or the masses of the ejecta and of the (compact) remnants, if there are any. 
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Indirect information can be obtained from the light curves and the spectra leading, however, 
to an "inverse" problem as far as theoretical interpretations are concerned. 



Theoretical Classification 

The theoretical classification of supernovae is usually done according to the suspected progen- 
itor stars and the explosion mechanism, and it dates back to a classic paper of Fred Hoyle and 
Willy Fowler in 1960 |^]. Based on very few observational facts available at that time, they 
postulated that Type II supernovae are the consequence of an implosion of non-degenerate 
stars, whereas Type I's are the result of the ignition of nuclear fuel in degenerate stars, and 
today this is still believed to be true in general, if "Type I's" are substituted by "Type la's" . 



Theoretical classification 



Progenitors 
massive stars /white dwarfs 




with without Mchan s\ab-Mchaii 

H-env. H-env. (C+O) (C+0;He) 




SN la 



Figure 1: Theoretical classification scheme of supernovae according to their progenitors. 

To be more specific, the present theoretical classification schemes based on progenitor 
properties and the energetics of the explosion are given in Figs. ^ and ^, respectively. It is 
now generally believed that Type II and Ib.c supernovae stem from collapsing massive stars 
with and without hydrogen envelopes, respectively, and that SN la originate from explosions 
of white dwarfs, although the way to explosion (accretion vs. merging of two white dwarfs) as 
well as the mass of the white dwarf (Chandrasekhar vs. sub-Chandrasekhar mass) just prior 
to the explosion are still heavily disputed (see Ref. [|[ and references therein). 

The classification with respect to the energetics is even more complex, and a variety of 
possible models is available with no unique and well accepted answer. For example, in the 
case of thermonuclear explosions, it is not even clear what the mode of propagation of the 
burning front is, and observations supply poor constraints only. A (fast) deflagration wave 
in a Chandrasekhar-mass white dwarf can explain the spectra and light curves of Type la's, 
but so can deflagrations changing into detonations at low densities with or without pulsations 
(so-called delayed detonations), or even pure detonations in stars with low enough densities 
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Figure 2: Theoretical classification scheme of supernovae according to their energy sources. 

are not completely ruled out. As far as the physics of core collapse supernovae is concerned, 
the situation is not significantly better. Whether or not rotation is important, whether or not 
magnetic fields play any role, whether neutrinos are diffusively or convectively transported 
from the cooling proto-neutron star to the stellar mantle, and many other questions are largely 
unanswered. 

Merging Theory and Observations 

Of course, a simple way to try to improve the rather unsatisfactory situation outlined in 
the previous section is to call for more observations. Since computing realistic theoretical 
light curves and synthetic spectra for a given supernova model has become feasible recently 
(although computing reliable spectra for Type la's still requires novel techniques) high quality 
data might help to rule out certain possibilities. But as long as most theoretical models contain 
a large number of more or less free parameters one may doubt the success of this approach. 

A few examples may serve as illustrations. As long as for a thermonuclear explosion the 
propagation velocity of the burning front and the degree of mixing can still be treated as free 
functions, it is easy to fit a given light curve and spectrum of a Type la supernova (see, e.g. 
Ref. but one does not prove that the model is correct. Similarly, in the core collapse 
scenario, neutrinos can transfer momentum to the mantle of a star and cause a supernova 
explosion, provided their luminosity is sufficiently high. If this luminosity results from low 
opacities or from convection does not matter. Also, only thousands of neutrino events from 
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a single supernova seen in a neutrino detector could possibly tell the difference. Finally, a 
similar momentum transfer can, in principle, be mediated by appropriate combinations of 
rotation and magnetic fields. Again, a lucky accident of a galactic supernova may supply 
enough data, including detections of gravity waves, to eliminate several models, but only for 
one special event and not for the entire class. 

Therefore, it appears to be more promising to search for ways to reduce the freedom one 
still has in building theoretical models. An obvious possibility is to replace parameter studies 
by first principle calculations, whenever this is possible. Again, a few examples are given 
to demonstrate that there is plenty of room for improvements. A first and obvious example 
is again neutrino transport in core collapse supernovae. One will note that very primitive 
(and likely incorrect) approximations for neutrino interactions with dense nuclear matter 
are commonly used in all hydrodynamic simulations, and only very recently first attempts 
have been made to compute those cross sections on the basis of microscopic theories 0. 
Similarly, well-developed methods to calculate the properties of dense nuclear matter are 
available, but the equations of state used in numerical studies are still computed ignoring 
nucleon-nucleon correlations, despite the fact that the stiffness of the equation of state is very 
important ingredient. As far as thermonuclear models are concerned, numerical methods to 
handle turbulent combustion are presently developed for combustion in engines and for reactor 
safety, and there is no principle problem in applying them also to supernovae. Moreover, 
simple microscopic models can be calibrated to laboratory combustion experiments, which 
again will help to reduce the number of parameters. 

In conclusion, not underrating the importance of new and better observational data in 
supernova research, I still think that most of the progress in the near future has to come and 
can come from theory. 
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Supernova Rates 

Bruno Leibundgut 
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Apart from stellar winds supernovae are the only mechanism which releases chemically pro- 
cessed material from stars into the interstellar and intergalactic gas. They are almost exclu- 
sively responsible for the chemical enrichment of galaxies and the universe as a whole. As 
one of the end stages of stellar evolution and one of the production channels of pulsars and 
black holes supernovae are also placed at a major link of stellar and non-stellar matter. They 
further inject kinetic energy into the interstellar gas which may be important for gas heating, 
cosmic rays, and possibly the kinematics of galaxies. 

Knowing the frequency of supernovae means to know the rate with which these processes 
occur. Having further knowledge of the temporal evolution of this frequency can provide the 
chemical history of galaxies. Supernova statistics are thus a main ingredient in the study of 
formation and transformation of matter in the universe. 

Recent reviews of supernova rates has been published by van den Bergh &: Tammann 
(1991), Tammann (1994), and Strom (1994). The visual searches have been summarized by 
van den Bergh and McClure (1994). A series of papers by Cappellaro et al. (1993a, b, 1997) 
have explored the derivation of the supernova rate in detail. SN statistics starting from the 
available SN catalogs and exploiting various indirect routes were presented by Tammann, 
Loffler, & Schroder (1994). 

Supernova statistics make use of minimal information about the explosions themselves. 
The main ingredients are the type of the explosion and a crude description of the parent 
population in which the explosion occurred. This is normally restricted to the morphological 
type of the host galaxy. The last piece of information is the time when the explosion took 
place. The frequency of supernovae fsN is described as the ratio of known SNe ATgN over the 
total galaxy luminosity Lg^i per unit time t 

NSN 
i^SN = J 7- 

The SN rate is normally expressed in Supernova Units (SNu) which corresponds to one su- 
pernova per 10^° (blue) solar luminosities per 100 years. 

Although in principle easy to determine, the derivation of supernova rates suffers from a 
number of problems. Supernovae are very rare objects. The number of supernovae in galaxies 
with known parameters (morphological type, luminosity, color) is rather small. There are 
about 60 to 80 supernovae per year, but only about 10 bright supernovae in nearby galaxies 
(cf. Timmes & Woosley 1997). The surveyed galaxy sample has to be well described and 
defined. Many supernovae occur in galaxies which are not part of the current catalogs, thus 
being lost for the derivation of the rates. Another problem is the control times, i.e. the time 
the galaxies have actually been surveyed for supernovae. Large corrections can be incurred 
when the galaxies are not observed frequently. 

The small number statistics becomes even more apparent when we consider that the known 
supernovae have to be split into a number of subsamples to become meaningful. Typically 
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there is a subdivision into three SN types and about five galaxy morphologies, which imme- 
diately creates 15 bins to group the supernovae. With typical numbers of about 200 SNe 
(e.g. van den Bergh & McClure 1994, Cappellaro et al. 1997) in the sample the rates are not 
immune against small number statistics. 

Over 1200 supernovae been discovered to date. Many of these objects are not classified 
and can not be used for the statistics. However, about 80% of all supernovae since 1989 have 
a type associated. The distribution is roughly 55% SNe la (thermonuclear explosions of white 
dwarfs) and 41% SNe II and SNe Ib/c (core-collapse in massive stars). Since SNe la are 
on average about 2 magnitudes brighter than other supernovae, they can be detected over a 
larger volume and thus are present at a higher percentage. Excluding all the supernovae from 
searches specifically targeted at very distant supernovae (Perlmutter et al. 1997, Leibundgut 
& Spyromilio 1997), which mostly are Type la, we arrive at a near equipartition between 
thermonuclear and core-collapse SNe (48% SNe la and 49.5% SNe II and SNe Ib/c since 1989). 
With the same argument as above (SNe la more luminous) we can immediately conclude that 
core-collapse SNe are more frequent than thermonuclear SNe in the nearby universe. This 
result is further amplified by dust obscuration which affects core-collapse SNe stronger than 
thermonuclear SNe. 

Several corrections enter the determination of supernova rates. The best known is the 
Shawn effect (Shawn 1979) which describes the effect of supernovae lost in the glare of centers 
of galaxies. Visual and CCD searches are less affected than older photographic searches which 
lose the contrast in bright regions. Absorption can obscure supernovae in the discs of galaxies 
and if, as currently believed, core-collapse supernovae stem from massive stars and arc thus 
more affected by dust, the derived SN statistics are skewed. Several proposals how to correct 
for dust have been put forward (e.g. van den Bergh & Tammann 1991). The luminosity of 
the various subtypes of supernovae differs by factors of about 10 and the differences in light 
curve shapes introduces varying control times for the SN classes. These have to be factored 
in when the SN rate is derived. 

There are two distinct approaches to derive SN statistics. The 'search technique' uses 
only supernovae which have been discovered in controlled searches. The most recent and 
complete description of this method is given by Cappellaro et al. (1997). The advantage here 
is that the control times and the galaxy sample are defined very well. The drawback is the 
small number of supernovae in the sample. Cappellaro et al. (1997) had to discard about half 
of the discovered supernovae, because they appeared in galaxies which are not part of the 
catalogs. They were left with 110 supernovae from five searches which have lasted for almost 
20 years and amounted to a total of 42500 years x IO^'^Lq control time (sometimes referred 
to as 'galaxy years'). 

The 'catalog technique' makes use of all supernovae discovered in a known galaxy sample 
(Tammann et al. 1994). This maximizes the number of supernovae as all supernovae discov- 
ered in galaxies of the sample are included in the statistics, but strong assumptions on the 
control time have to be made. The latter quantity can be derived very badly as serendipitous 
SN discoveries enter the catalogs. 

Despite all these problems, there is a fair agreement for the rates of thermonuclear su- 
pernovae (Type la) in star-forming (spiral) galaxies. The distribution is independent of the 
detailed morphological type of the galaxy at about 0.2 SNu. A large discrepancy exists for 
old stellar systems (elliptical galaxies) for which Tammann et al. find a rate increasing to 
a level three times higher, while Cappellaro et al. actually find indication of slightly lower 
SN la rate in ellipticals. There is only a very small number of supernovae known in elliptical 
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galaxies, however, and the statistics are very uncertain. 

Core-collapse supernovae (Types II and Ib/c) have not been observed in elliptical galaxies. 
Their rate increases with the fraction of young stars in galaxies, i.e. the morphological type. 
The most prolific supernova producers are Sc and Sd galaxies with about 1.2 SNu. There 
is a discrepancy of about a factor of two in the absolute rates between Tammann et al. and 
Cappellaro et al. for all galaxy types. The relative rates, however, agree quite well. The SN 
rates in irregular galaxies are largely undetermined due to the small numbers of known SNe. 

The Galactic SN rate is about 20 it 8 per millennium. This number depends on the mor- 
phological type and the total blue luminosity of the Galaxy. There are 6 historical supernovae 
on record (van den Bcrgh & Tammann 1991) which implies that roughly 2/3 of all Galactic 
supernovae have been hidden. Supernovae in the Local Group have been observed in M31 
(S And: SN 1885A) and in the Large Magellanic Cloud (SN 1987A). The prediction for the 
Local Group is slightly above 2 (excluding the Galaxy), which appears consistent. The small 
observational baseline, however, excludes any firm conclusion from such a comparison. 
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Pulsar Velocities and Their Implications 
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In several respects, the Galactic distribution of pulsars mimics the distribution of popula- 
tion I species, such as young stars, HII regions and supernova remnants. In particular, the 
galactocentric radial distributions are very similar. However, the widths of the distributions 
in distance from the Galactic plane, Z, differ by nearly a factor of 10. Gunn and Ostriker 
(1970) Q realised that this difference could result from high velocities which may be given 
to pulsars at birth, possibly arising in the violence of their formation in supernova collapse. 
It was several years later before the first measurement of the proper motion of a pulsar gave 
direct evidence for this hypothesis |Q. Transverse velocities are now known for more than 
100 pulsars and their high values are generally well established. 

In this review, I describe briefly the techniques used in the measurement of pulsar ve- 
locities, summarise the main features of the distribution of the velocities and discuss their 
implications for the physics of the formation events and for their spatial distribution. 

There are two main techniques which have been used for the determination of pulsar 
velocities, in both cases providing the transverse components of the velocity only. These 
techniques respectively involve the measurement of proper motion and of the velocity of the 
interstellar scintillation pattern across the Earth. Proper motions are determined from high- 
resolution radio interferometry or from timing measurements. Like optical techniques, the 
former involves measurement of the position of a pulsar relative to a nearby source over a 
period of time. With modern instruments, position accuracy of 10 milliarcseconds (mas) is 
readily achievable, so that only a few years are required to give proper motion errors of only a 
few mas/yr |^, ^, ^, Timing measurements can give similar precision in principle, but are 
usually limited by irregularities in the rotation rate of the pulsar, known as "timing noise," 
which is predominant in young pulsars. However, timing is particularly useful for millisecond 
pulsars which have very stable rotation and whose narrow pulses allow very high accuracy. 

Estimates of pulsar transverse velocities Vt are then obtained from the proper motion ^ 
(mas/yr) and the distance D (kpc): Vt = 4.74 ^uL) km/s. The distance D is obtained from 
the dispersion measure and an electron density model Q which is based upon a number of 
independent distance measurements. Nearly 100 such measurements of transverse velocity 
are available. 

The velocities of pulsars estimated from observations of the speed of the interstellar scin- 
tillation patterns [0, |l3] are considered rather less reliable than the more direct proper motion 



measurements, mainly because it has recently been shown |11] that the values are systemati- 
cally low by an average factor of 2 because of a localisation of the scattering medium close to 
the galactic plane. There are about 71 pulsars for which scintillation speeds are available ||lO|| , 
although most of these pulsars have rather better proper motion measurements, leaving 27 
which have only scintillation measurements. Combining both techniques gives a sample of 
about 100 values of 14, of which about 8 are upper limits. 

From this sample, the mean transverse velocity is found to be (V^) = 300 it 30 km/s. 
With a small number of exceptions, the velocity vectors shown in Fig. || demonstrate a clear 
movement away from the galactic plane, consistent with pulsars being formed from massive, 
young, population I stars close to the plane and receiving velocities of a few hundred km/s at 
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Figure 1: The motions of pulsars relative to the galactic plane. Pulsars are represented as 
filled circles in galactic longitude and Z-distance, the tail representing its approximate motion 
in the past 1 Myr. 



the same time. How large are the typical space velocities? Unfortunately, there is a selection 
effect visible in the data in which high velocity pulsars quickly move away from the plane of the 
Galaxy and become undetectable, leaving an excess of low- velocity pulsars. This effect can be 
minimised by considering the velocity distribution of only the youngest pulsars, which cannot 
have moved far since birth. The 29 pulsars with characteristic age of less than 3 Myr give 
{Vt) = 350 lb 70 km/s. Monte Carlo simulations show that the mean space velocity required 
to give such an observed transverse velocity distribution is (Vg) = 450 it 100 km/s [^]. 

There are two other pieces of evidence which support such high velocities: the increase 
of the mean distance of pulsars from the plane as a function of characteristic age and the 
movement of a small number of pulsars from the centres of the putative supernova remnants 
associated with their birth 

Similar analyses for the 13 millisecond pulsars with proper motion measurements give 
(Vt) = 88 lb 19 km/s and {Vs) = 130 it 30 km/s, showing them also to be a rather high-velocity 
population. Not surprisingly, no migration from the galactic plane is seen, because of their 
large ages. These velocities probably represent the low- velocity end of the "kick" spectrum — 
higher velocity kicks are likely to disrupt any binary system, preventing any millisecond pulsar 
formation by accretion spin-up in a binary system. 

Clearly, the high velocities explain why only a few percent of pulsars are in binary systems, 
compared with the high proportion of the likely progenitor stars in binaries. However, the 
high velocities pose something of a problem for the retention of pulsars in globular clusters: 
the velocities mostly exceed the cluster escape velocities and yet there is a substantial pulsar 



56 



population which has remained gravitationally bound |[l3| , [T^. It is also worth noting that 
about half of the pulsar population will escape the gravitational potential well of the Galaxy, 
ending up in intergalactic space, the remainder being in a large galactic halo. 

There has been much speculation on the origin of the high velocities, the three main 
possibilities being the disruption of binary systems, an electrodynamic "rocket" effect in the 
early life of the pulsars and asymmetry in the supernova explosions. The first two have 
great difficulty in explaining the size of the velocities and the most likely origin seems to 
lie in the supernova explosion itself. The high momentum asymmetry of the neutron star 
after the explosion must be matched by that of either the ejecta or neutrinos. The origin 
of the asymmetry is not clear, although both convection and the magnetic field have been 
cited as the possible determining agents. It is worth noting that up to now, the magnetic 
properties of the neutron stars do not seem to implicate the magnetic field in the production 
of the high velocities. 
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Introduction 

Observations of SN 1987A instigated recent work on multi-dimensional numerical modeling 
of Type-II supernova explosions. Two-dimensional (2D) simulations |2|, |3|, ^ indicate that 
overturn instabilities and mixing in the i/-heated region between PNS and supernova shock 
may indeed help i^-driven explosions to develop for conditions which do not allow explosions 
in spherical symmetry. Nevertheless, a critical investigation of the sensitivity of the explosion 
to the value of the v flux from the neutron star j^, Q cannot confirm claims that a "convective 
engine" ensures the robustness of the explosion for a wide range of conditions. Instead, the 
explosion turns out to be extremely sensitive to the u fluxes from the inner core. Explosions 
can occur in 2D as well as in ID, provided the v luminosity is large enough. Similarly, for 
too small core v fluxes, strong convective overturn in the z^-heated region cannot develop 
and the explosion fizzles. These findings are supported by recent simulations of Mezzacappa 
and collaborators |^]. In addition, all currently successful numerical models of supernova 
explosions produce nucleosynthesis yields that are in clear contradiction with observational 
abundance constraints. These facts indicate severe problems of the numerical modeling of 
the explosion. The first 2D simulations of the evolution of the nascent neutron star |^ 
suggest that long-lasting quasi-Ledoux convection in the PNS can significantly raise the v 
luminosities and may have a positive impact on the supernova explosion and conditions for 
explosive nucleosynthesis. 

Results of Two-Dimensional Hydro dynamical Simulations 

Our simulations were performed with an elaborate description of the nuclear equation of 
state and the v transport and were started with a PNS model of Bruenn Q at ~ 25 ms after 
core bounce. These calculations showed that driven by negative lepton fraction and entropy 
gradients convective activity can encompass the whole PNS within ~ 1 s and can continue for 
at least as long as the deleptonization of the PNS takes place. Because of the importance of 
the diffusive v transport this convection is not an ideal, adiabatic Ledoux convection. Also 
specific thermodynamical properties of the nuclear equation of state have to be taken into 
account. In differentially rotating PNSs stabilizing angular momentum distributions seem to 
suppress convection near the rotation axis. 

The convective pattern is extremely non-stationary and has most activity on large scales 
with radial coherence lengths of several km up to ~ 10 km and convective "cells" of 20°-30° 
angular diameter, at some times even 45° (Fig. ||a). The maximum convective velocities are 
~ 5 X lO^cm/s, but peak values of ~ lO^cm/s can be reached. Relative deviations of the 
lepton fraction l^ep from the angular mean can be several 10% in rising or sinking buoyant 
elements (Fig. ||b), and for the entropy {S) they can reach 5% or more. Rising flows always 
have larger l^ep Oind S than their surroundings. 
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Figure 1: Panel a shows the absolute values of the convective velocity for a non-rotating 
PNS model 0.7 s after the start of the simulation in units of 10^ cm/s. The computation 
was performed in an angular wedge of 180° assuming axisymmetry. Panel b displays the 
relative deviations of the lepton fraction Y|ep from the angular means (llep) at each radius for 
t = 0.7 s. The maximum deviations are of the order of 50%. Lepton-rich matter rises while 
deleptonized material sinks in. 



Since convective activity continues for more than 1 s in a large region of the PNS, convec- 
tive I' transport shortens the cooling and deleptonization timescales of the PNS compared to 
results of ID simulations and enhances the calculated v luminosities of the star by up to 65%. 
The latter can have important influence on the interpretation of the measured i' signal from 
SN 1987A and might change limits of various quantities in nuclear and elementary particle 
physics which have been derived from these measurements. Moreover, the enhancement of the 
v fluxes can also have important consequences for the explosion mechanism of the supernova, 
because it aids v-driven explosions. The faster deleptonization modifies the luminosity ratio 
of fe and z^e in such a way that the electron fraction Y^^ in the z/-heated SN ejecta will be 
raised during the first few 100 ms after core bounce, but it will be lowered for t ^ 1 s com- 
pared to ID simulations. This might help to solve the severe problems of the nucleosynthesis 
in current models of Type-II supernovae which disregard long-lasting convective activities in 
the PNS. Finally, convection in the PNS causes stochastical asymmetries of the i' flux. In 
our calculation the resulting recoil accelerates the PNS to ~ 9 km/s during 1.2 s (Fig. ^). 
This is too small to explain measured proper motions of pulsars of a few 100 km/s. Both 
convective mass motions and anisotropic v emission are a source of gravitational waves for 



59 



ft ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' fl 




0.2 0.4 0.6 0.8 1 1.3 10 100 1000 10" 

t [s] f [Hz] 

a b 

Figure 2: Panel a shows the recoil velocity of our PNS model in km/s as a function of 
time resulting from stochastical asymmetries in the v flux. Panel b displays the spectral 
energy density of the emitted gravitational radiation caused by mass motions and anisotropic 
V emission. 

more than 1 second (Fig. ^). Most of the gravitational radiation is emitted at frequencies of 
300-1400 Hz. 
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Abstract 

Hydrodynamic instabilities during the supernova explosion and neutrino cooling phase lead 
to stochastic acceleration of the nascent neutron star and are hardly able to account for recoil 
velocities >100km/s. It is argued that an internal magnetic field of ~ 10^^ G can define a 
preferred direction of the neutrino emission with an anisotropy of several per cent which is 
sufficiently large to produce kicks even in excess of 1000 km/s and thus can explain the fastest 
motions of observed pulsars. 

Introduction 

Pulsars are observed with large mean velocities of several 100 km/s (see A. Lyne, this con- 
ference, for references to observations), and a few neutron stars are measured with velocities 
of more than 1000 km/s. Similar velocities can also be deduced from associations of neutron 
stars with nearby supernova remnants. Binary evolution arguments strongly suggest that 
neutron stars experience a recoil already during their birth in the supernova explosion of a 
massive progenitor star (e.g., ^). The reason for this acceleration is not finally under- 
stood and several mechanisms were proposed or investigated, among them anisotropic core 
collapse and mass ejection during the supernova explosion (e.g., [^), large-scale hydrodynamic 
instabilities during the explosion (e.g., |Q, ^), or anisotropic neutrino emission of the cool- 
ing proto- neutron star for example associated with convective processes (e.g., [§, 0)) 
magnetic fields (e.g., |^, or neutrino oscillations (e.g., |jl^, |ll|]). A 1% anisotropy of the 
neutrino emission during the neutrino cooling phase of the newly formed neutron star would 
be sufficient to accelerate the remnant to about 300 km/s. 

While some of the proposed kick mechanisms are highly speculative and involve one or 
even more components of new or non-standard physics (e.g., neutrino masses or magnetic 
moments, superstrong magnetic fields), other mechanisms are unsatisfactory in the sense 
that they rely on certain initial conditions in the collapsing stellar core or in the nascent 
neutron star which are very specific and cannot be easily justified from stellar evolution 
calculations (e.g., special magnetic field configurations or density distributions to produce 
global anisotropics). More "natural" anisotropics which develop in self-consistent simulations 
from hydrodynamic instabilities without the use of specially chosen initial conditions, on the 
other hand, turned out to be sufficient to account for small velocities of a few 10 km/s but 
seem hardly able to explain the very fast motions of pulsars. 

Hydrodynamic Instabilities and Neutron Star Convection 

For example, anisotropic mass motions due to convection in nascent neutron stars lead to 
gravitational wave production and anisotropic emission of neutrinos (see W. Keil, this con- 
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ference) . The angular variations of the neutrino flux determined by 2D simulations are of the 
order of 5-10% ||l^, [l^. With the typical size (10-30 degrees angular diameter) and short co- 
herence times (a few milliseconds) of the convective structures, however, the global anisotropy 
of the neutrino emission from the cooling proto-neutron star is certainly less than 1% (more 
likely only 0.1%) and observed high kick velocities in excess of 300 km/s can definitely not be 
explained. In the simulation described by Keil, the determined recoil velocity is disappoint- 
ingly low, only about 10 km/s! There is no reason to expect this case to be an especially 
unfavorable one. Even worse, in three dimensions the convective elements tend to become 
smaller than in 2D because of fragmentation and turbulent energy transport from larger to 
smaller scales which is forbidden in 2D due to angular momentum conservation. This makes 
large recoil velocities of neutron stars by the described mechanism even more unlikely. 

Convective motions in the neutron star have important consequences for the magnetic field 
structure. Specific kinetic energies of the order of ^v"^ ~ 10^^ erg/g could lead to equipartition 

fields of more than B ~ 10^^ pj^^ G. Such strong fields increase the neutrino opacities Q 
and could affect the anisotropy of the neutrino emission. In the dense medium of the neutron 
star, where the electrons are degenerate, the influence of the magnetic field on the charged- 
current reactions becomes essential when the characteristic energy of an electron on a Landau 
level reaches the size of the electron chemical potential, i.e. when B^Bc • (/ig/w-eC^) ~ 
5 X 10^^ ... 2 X 10^^ G where Be = mlc^/{he) = 4.4 • 10^^ G. The presence of such strong 
fields changes the phase space distribution of the electrons and thus the cross section for 
absorption on neutrons Whether this results in sizable neutron star kicks or not will 

depend on the strength and structure of the magnetic fields inside the neutron star. Turbulent 
convective motions will certainly produce a much more irregular magnetic field distribution 
than the idealized mirror asymmetry assumed by Bisnovatyi-Kogan who suggested that 
an initial, strong toroidal field component is amplified by a wound-up poloidal field in one 
hemisphere, whereas both toroidal field components (the initial one and the one produced 
by winding the poloidal field) are superposed destructively in the other hemisphere. A more 
irregular field evolution will imply a much smaller global anisotropy of the neutrino emission 
than estimated in reference [^. 

Rapid rotation of the collapsed stellar core has very important effects on the develop- 
ment and strength of convection. Two-dimensional, self-consistent hydrodynamic simula- 



tions 1 12, 13 1 reveal that a positive gradient of the specific angular momentum in the nascent 
neutron star stabilizes hydrodynamically unstable entropy or lepton number stratifications. 
This positive gradient of the specific angular momentum suppresses convective motions near 
the rotation axis and allows strong convection to occur only close to the equatorial plane. 
Therefore, the neutrino emission is convectively enhanced around the equatorial plane while 
it is essentially unchanged near the poles. For the same reason, the magnetic field structure 
may retain an initial dipole component and irregular fields must be expected to be produced 
by convective mass motions only in a belt around the equator. 

Parity Violation and Cumulative Asymmetry 

According to a recent suggestion by Horowitz and Li jl^ , a global anisotropy of the neutrino 
emission of the nascent neutron star could originate from the cumulative effect of a large num- 
ber of neutrino scatterings off polarized nucleons in the magnetized stellar medium. Although 
the polarization corrections to a single neutrino-lepton interaction are small, typically of the 



order of only 10 of the cross section [15, 16 1, Horowitz and Li recognized that a macrO' 
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scopic asymmetry can build up because the anisotropy of the neutrino flux increases with the 
average number of scatterings per neutrino, i.e. with the optical depth of the medium. The 
polarization term in the cross section represents a slightly different chance for a neutrino to be 
scattered into the direction of the external magnetic field than into the opposite direction, or 
means a tiny enhancement of the scattering probability of neutrinos moving into (or opposite 
to) the field direction. Formally, the neutrino flux is not only described by a diffusive propa- 
gation mode (proportional to the gradient of the neutrino energy density), but an additional 
"advective" vector component in the magnetic field direction (see [ p!7| , ^]). The ratio between 
advective and diffusive component grows approximately linearly with the scattering optical 
depth r. Physically, this can be understood |15| by recalling that a diffusing particle has 



travelled a mean distance d = \/ NX after N steps with (constant) scattering mean free path 
A. If the particle has a tiny chance P to be scattered into the field direction, the distance 
it has propagated along the field after scatterings is / = PNX. Therefore the relative 
anisotropy increases as l/d oc Py/N oc r (for Py/N ^ 1). The second proportionality holds 
because the optical depth is defined as r = d/X, thus r = ^/N. 

The discussed process would be highly interesting if the magnetic fields required to create 
a few per cent anisotropy of the neutrino emission of a neutron star were much smaller 
than the very strong fields ( > several 10^^ G) where the neutrino opacity is affected by the 



change of the electron phase space distribution. Horowitz and Li [15| have only considered 
very simple situations of scattering media to show the fundamental characteristics of the 
cumulative process and have disregarded essential complications of neutrino transport in 
neutron stars, e.g., neutrino absorption, transport of different neutrino types with opposite 
signs of the polarization terms, and the feedback of the transport on the energy distribution 
in the neutron star. The question must be asked whether the proposed cumulative process 
leads to a global anisotropy of the neutrino emission of the nascent neutron star and how 
large this anisotropy can be for a given magnetic field. 



The neutron star situation was more realistically investigated by Lai and Qian |18] who es- 
timated a maximum recoil velocity of about 200km/s for average magnetic fields of ~ 10^^ G, 
provided the magnetic field in a nascent neutron star possesses a strong dipole component. 
Their analysis, however, is based on the assumption that neutrinos and antineutrinos, in par- 
ticular Ue and z/e, have the same opacity for interactions with the stellar medium. Therefore 
their kicks are only produced by the deleptonization neutrinos which carry away only a minor 
fraction (~ 10%) of the gravitational binding energy of the neutron star which is released by 
neutrino emission. This assumption, however, leads to a large underestimation of the recoil 
acceleration p9| . 

Monte Carlo simulations |19| which include neutrino production and absorption, different 
opacities and polarizations for different neutrino and antineutrino flavors (i^e; i^e; and heavy 
lepton neutrinos Vx and v'x), and the evolution of the stellar background in a simplified but rea- 
sonable, approximative description^, suggest that a polarization term in the neutrino-nucleon 
scattering opacity of Ks,poi/'^s ~ 10~^ can lead to a global asymmetry q of the integrated 
neutrino luminosity of about 16% for the one-dimensional (plane) case (Fig. 1). Taking 
into account projection effects by multiplying this result with 1/2, the corresponding three- 
dimensional asymmetry qsn is estimated to about 8%, i.e. qsd ~ 0.08- [{ks,poi/i<-s) /10~^] . This 



^Monte Carlo simulations become extremely expensive and CPU-time consuming for high optical depths. 
Instead of considering a neutron star with neutrino optical depth of 10000-100000, a toy model with optical 
depth r — 100 and correspondingly rescaled polarization terms of the neutrino interactions was investigated. 
The scaling of the results was tested by comparative computations for t — 30 and for r — 300. 
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Figure 1: Neutrino luminosity and emission anisotropy as functions of time step for two Monte Carlo 
simulations of the neutrino cooling of a one-dimensional toy model of a neutron star. L = L+ + L- is 
the total neutrino luminosity, AL = h-\. — L- AL = L+ — L- the luminosity difference between the two 
hemispheres (both normalized to the maximum luminosity Lmax), E — E-\- + E- the energy radiated in 
neutrinos, and AE = E+ — E^ the energy difference (both normalized to the final energy i^max)- The 
integral asymmetry of the neutrino emission is given by the parameter q = AE/E which approaches 
nearly 16% towards the end of the simulations. The three-dimensional asymmetry in the integrated 
neutrino luminosity, (J3d, would be smaller by a factor 1/2 (for small anisotropies) compared to the 
displayed one-dimensional results. The Monte Carlo simulations were performed for models with total 
optical depths of r = 30 (left) and r — 100 (right), where the opacity polarizations were rescaled such 
that the results should be identical if the anisotropy grows linearly with the optical depth. The very 
good agreement between both calculations (left and right) thus confirms the scaling of the cumulative 
anisotropy with the optical depth of the star. 



can be related to the anisotropy parameter for the momentum transfer to the neutron 
star: 030 = §930 ~ 0.05 • [(ks,poi/'*s) /10~^]. The effective polarization of neutrino-nucleon 
scatterings depends on the magnetic field strength B, the plasma temperature T, and the 
abundances Yn and Yp of neutrons and protons, respectively: 

- ■ ^6.4 X 10-5 - 13.3 X 10-5^) • f 1 + 1.134 ^V' • , (1) 



YnJ \ YnJ Tio 

where -B14 is measured in 10^^ G and Tiq in lOMeV. Using representative values for tem- 
perature and composition during the phase when most of the gravitational binding energy 
of the neutron star is radiated away in neutrinos, T = 10...30MeV, Yn = 0.8. ..0.9, and 
Yp = 1 — Yn = 0.2 ... 0.1, a space-time average of Kg poi/^s ~ 10"^ can be expected for a mean 
interior magnetic field {Bz) of the neutron star of ~ 10^^ G. This leads to an estimated recoil 
velocity of 



0.05/ V / V3 X 10^3 erg 

The reason for a nearly 10 times larger anisotropy compared to the result found by Lai 



and Qian |18l is twofold. During deleptonization and neutronization the neutron star plasma 



becomes increasingly asymmetric, i.e. the abundance of protons decreases and the plasma 
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becomes neutron-rich. On the one hand, the unequal concentrations of neutrons and protons 
cause different opacities for Vf, and Pg; on the other hand the growing abundance of neutrons 
imphes an increase of the total polarization of the nucleonic medium because the contributions 
from neutrons and protons have opposite signs and partly cancel. For a typical post-collapse 
composition, Yp 0.35 and Yn ~ 0.65, this cancellation is severe and reduces the polarization 
to roughly 1/8 of the value for a neutronized medium with Y^ ~ 0.1 and Y^ ~ 0.9. The latter 
composition is representative of the late Kelvin-Helmholtz cooling phase after the deleptoniza- 
tion of the proto-neutron star. Most (60-90%) of the neutrino energy is radiated during this 
phase which is therefore decisive for a "rocket engine" effect by asymmetric neutrino emission. 

Of course, the estimated velocity in Eq. (^) is an upper limit which could be realized only 
if the interior magnetic field were uniform in the z-direction throughout the neutron star. 
Convective processes that reach deep into the star will certainly destroy the ordered structure 
of the field. In case of rapid rotation of the newly formed neutron star (rotation periods of 
a few milliseconds), however, convection can develop only near the equatorial plane but is 
suppressed near the rotation axis where an ordered field structure could persist in a large 
fraction of the neutron star volume. 



Conclusions 

If parity violation of weak interactions plays the crucial role to kick pulsars and if convection 
is important during the early phases of a neutron star's life, one would therefore conclude 
that the observed velocities of fast pulsars between several lOOkm/s and more than 1000 km/s 
require interior magnetic fields from a few 10^^ G to ~ 10^^ G, and that the largest velocities 
might indicate rapid rotation of the newly formed neutron star. 
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Abstract 

Neutrinos emitted from the supernovae which have ever occurred in the universe should make 
a diffuse neutrino background at present. This supernova reUc neutrino background (SRN) 
is one of the targets of the Supcrkamiokandc (SK) detector which wiU be constructed in 
this year and the SRN, if at all detected, would provide a new tool to probe the history of 
supernova explosions in the universe. The SRN spectrum is calculated by using a realistic, 
time-dependent supernova rate derived from a standard model of galaxy evolution based on 
the population synthesis method. The expected event rate at the SK is also calculated. The 
SRN spectrum we show here is the most realistic at present, because the largest uncertainty 
in previous theoretical predictions has come from unrealistic assumptions of the supernova 
rate so far made. Our major results include: (1) the supernova rate is much higher in the 
early phase of evolution of galaxies and there appears a hump in the SRN spectrum in the 
low-energy region of ^ 5 McV, (2) the SRN flux depends on the Hubble constant {Hq) in 
a way approximately proportional to Hq and only weakly on the density parameter of the 
universe (fio) and a cosmological constant (Aq), (3) the uncertainty in the star formation 
history of spiral galaxies affects the resulting SRN flux by about a factor of 3, and (4) the 
plausible event rate at the SK is 1.2 yr^^ in the observable energy range of 15-40 McV. Such 
a low event rate is due mainly to a quite low supernova rate at present which is averaged over 
the morphological types of galaxies. The most optimistic rate in our model is found to be 4.7 
yr~^ in the same energy range, and if more events are detected, we will have to reconsider 
our current understanding of collapse-driven supernovae and evolution of galaxies. 

Summciry and Conclusions 

We have calculated the SRN spectrum and also the expected event rate at the SK detector by 
using a realistic, time-dependent supernova rate from a reasonable model of galaxy evolution 
based on the population synthesis method. This model of galaxy evolution has been used to 
investigate the geometry of the universe from comparison with the number counts of faint 
galaxies, and the predicted number counts are well consistent with the recent observations. In 
our analysis the absolute SRN flux in the model is determined from the luminosity function 
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and the mass-to-luminosity ratio of nearby galaxies, but not from the present supernova rate 
of our Galaxy which is subject to a large margin of uncertainty. 

In the model we used, the supernova rate is much higher in the early phase of elliptical 
galaxies, and more than half of total supernovae explode during the initial 1 Gyr from the 
formation of galaxies. Thereafter, the supernova rate in spiral galaxies becomes dominant 
and the total number of supernovae until present is consistent with the nucleosynthesis re- 
quirement. A constant rate of supernova explosions, which was assumed in most of previous 
papers, is not justified from a viewpoint of galaxy evolution, and therefore the SRN spectrum 
we show here is the most realistic at the present time. 

Supernovae produced at a very high rate in the early phase of ellipticals are responsible for 
the appearance of a hump in the low-energy part 5 MeV) of the SRN spectrum, because 
the energy of neutrinos from those supernovae is degraded by the cosmological redshift effect. 
This energy degradation makes early neutrinos unobservable and the high supernova rate 
causes almost no enhancement in the event rate in the observable energy range of 15-40 MeV 
at the SK. 

We used a low-density, flat universe with a non-zero A as a standard cosmological model 
{O,Q=0.2, Aq = 0.2, and h = 0.8), and investigated the SRN spectrum for other cosmological 
models relative to the standard model. The SRN flux strongly depends on Hq as roughly 
proportional to Hq. The SRN flux changes with and A in a way that smaller and/or 
larger A suppresses the flux in the high-energy region lOMeV), however, these changes 
are rather small. The shape of the SRN spectrum is almost insensitive to either of the above 
three parameters. 

We investigated the uncertainty of the SRN spectrum from the model of galaxy evolution. 
The simple, one-zone and n = 1 (SI) model, which we use as a standard evolution model 
of spiral galaxies, gives an intermediate SRN flux among the four different models of spiral 
galaxies considered in this paper (for explanation of the models see Refs. [|^, Q)- The difference 
between the highest flux from the II model and the lowest flux from the S2 model is within 
a factor of 3. In the case of elliptical galaxies, however, the different SFRs hardly change the 
SRN spectrum in either shape or intensity. The effect of changing zp on the SRN is almost 
negligible above 15 MeV, and from zp = 3 to 5 the event rate at the SK changes only by a 
factor of 1.3. 

The total SRN flux in the whole energy range becomes 44 cm^^s^^ in our 'standard' model 
(ilo = 0.2, Aq = 0.8, h = 0.8, zp = 5, the SI model for spirals, the standard SFR model for 
ellipticals, and Cb = 1-93 x lO^h LbqMpc~^), which is a moderate value compared with the 
previous estimates, however, the SRN flux in the observable energy range of 15-40 MeV is 
as small as 0.83 cm~^s~^ corresponding to a low end of the SRN flux estimated before. This 
comes from the very small present supernova rate which is averaged over morphological types 
of galaxies with the weight of the mass-to-luminosity ratio. This small flux directly leads to 
a low event rate, and makes the SRN detection at the SK more difficult. 

The event rates at the SK detector are calculated for the various models considered in this 
paper (table |l]), and in the case of the 'standard' model, the expected rate is as small as 1.2 
yr"'^ in the observable range of 15-40 MeV. The event rate depends mainly on the evolution 
model of spiral galaxies, the Hubble constant, and the luminosity density of galaxies, while 
only weakly on the geometry of the universe (i.e., on and Aq) and the epoch of galaxy 
formation. These properties are in contrast with those of the number counts of faint galaxies, 
and the SRN observation may be complementary to the number counting of faint galaxies, 
if the SRN is detected. From the dependence of the flux on the Hubble constant, the SRN 
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signal would provide a new test for measuring the Hubble constant. However, because of 
presumably low event rates, such a test may not be possible in practice. 

Table 1: Expected Event Rates at the Superkamiokande Detector. 



Event Rate (yr ^) 
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" The standard SFR model is used in common to represent the E/SO-type evolution. 
^ The luminosity density of galaxies in units of 10® /i x Lbq Mpc~^ where 

h = iJo/lOOkms-^Mpc^^ 

Kinetic energy of recoil positrons. 
Note. — We refer the 'standard' model to the case with [VIq, Aq, h, zp, S-type evolution, C-b)= 
(0.2, 0.8, 0.8, 5, SI, 1.93), and the 'most optimistic' model to the case given in the last row 
with (0.2, 0.8, 1.0, 3, II, 1.93+0.8). 

Changing the model parameters within the allowable range in a way of increasing the 
event rate {h I, zp ^ 3, Sl^ II, and Cb (1.93 + 0.8) x lO^h LbqUpc'^), we predict 
the 'most optimistic' event rate at the SK, which is 4.7 yr~^ in the energy range of 15 -40 
MeV. The very small rate of expected SRN events at the SK implies the difficulty of gaining 
some insight on the star formation history in galaxies from the SRN observation in the SK 
experiment. More events exceeding 4.7 yr~^, if detected, may enforce a serious revision in 
our current understanding of collapse-driven supernovas and evolution of galaxies. 
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Figure 1: Time variation of supernova rates in a galaxy with the mass of 10 for various 
galaxy types derived from the standard model of galaxy evolution based upon the population 
synthesis method. The time is elapsed from the epoch of galaxy formation. The dashed lines 
show the rates for the individual galaxy types, and the thick line for the weighted average 
over the types. The SI model is used to represent the evolution of spiral galaxies, and the 
standard SFR model for elliptical galaxies. 



70 



E/SO 




0.01 0.1 1 10 



Evolutionary Time [Gyr] 



Figure 2: Cumulative numbers (i.e., the time integration of Fig. 1) of supernovae in a galaxy 
with the mass of 10^^ Mq as a function of time elapsed from the epoch of galaxy formation. 
The dashed lines correspond to the cases for the individual galaxy types, and the thick line 
to the weighted average over the types. The SI model is used to represent the evolution of 
spiral galaxies, and the standard SFR model for elliptical galaxies. 
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Figure 3: Energy spectra of the supernova relic neutrino background (SRN) for various cos- 
mological models. All lines represent the differential number flux of electron antineutrinos 
(i^e's). The thick Unes show the SRN spectrum for a low-density, flat universe with ^Iq = 0.2, 
Aq = 0.8, and three different values of the Hubble constant; h = 0.5 (bottom line), h = 0.8 
(middle line), and h = 1.0 (top line). The dashed line shows the SRN spectrum for a low- 
density, open universe with J7o = 0.2, Aq = 0.0, and h = 0.8, and the dot-dashed line for an 
Einstein-de Sitter universe with f)o = 1-0, Aq = 0.0, and h = 0.5. The SI model is used to 
represent the evolution of spiral galaxies, and the standard SFR model for elliptical galaxies. 
The epoch of galaxy formation is set to be zp = 5. For an illustrative purpose of comparison, 
two SRN spectra based on a model with a constant supernova rate are also shown by the 
dotted lines. The supernova rates of the two curves are chosen as follows. In the case of the 
upper curve, the time-integrated number of supernovac is the same with that of our standard 
model (thick, middle line), and in the case of the lower curve, the present rate of supernovae 
is the same with that of the standard model. 
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Figure 4: Expected event rates at the Superkamiokande (SK) detector for various cosmological 
models, as a function of kinetic energy of positrons produced by the reaction VeP e^n. This 
figure corresponds to Figure 3 of the SRN spectrum. The thick lines show the event rate for 
a low-density, flat universe with f^o = 0.2, Aq = 0.8, and three different values of the Hubble 
constant; h = 0.5 (bottom line), h = 0.8 (middle line), and h = 1.0 (top line). The dashed line 
shows the event rate for a low-density, open universe with Qq = 0.2, Aq = 0.0, and h = 0.8, 
and the dot-dashed line for an Einstein-de Sitter universe with = 1.0, Aq = 0.0, and 
h = 0.5. The SI model is used to represent the evolution of spiral galaxies, and the standard 
SFR model for elliptical galaxies. The epoch of galaxy formation is set to be zp = 5. For an 
illustrative purpose of comparison, two event rate curves based on a model with a constant 
supernova rate are also shown by the two dotted lines. The used rates of supernovae are the 
same with those used in Fig. 3. 
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The collapsed cores of supernovae are so dense and hot that neutrinos are trapped. There- 
fore, the transport of energy and lepton number is governed by the neutrino transport coeffi- 
cients, i.e. by their effective scattering and absorption cross sections. While electrons cannot 
be entirely ignored, it is the interaction with nucleons which is thought to dominate the neu- 
trino opacities. The cross sections can be dramatically modified by medium effects; strong 
magnetic fields which can also be important will not be considered here. 

Final-state blocking by degenerate nucleons is one obvious modification of the vacuum 
cross sections Even this trivial effect depends sensitively on the nucleon dispersion rela- 
tion; a reduced effective nucleon mass increases the degree of degeneracy. Moreover, depending 
on the equation of state the composition may deviate from a naive proton-neutron mixture in 
that there may be a significant fraction of hyperons. Studies of opacities with self-consistent 
compositions and dispersion relations reveal large modifications |2|. 

The neutrino interaction with nucleons is dominated by the axial vector current, i.e. it 
is a spin-dependent phenomenon. At supernuclear densities one thus expects a dominant 
role for spin-spin correlations |Q; the use of uncorrelated single-particle states for an opacity 
calculation amounts to an uncontrolled approximation. 

A less intuitive effect is the cross-section reduction by the nucleon spin fluctuations which 
are described by the dynamical dependence of the spin-density structure function [Q, ^, |^. 
In Fig. f the reduction of the average neutrino-nucleon axial-vector cross section is shown 
as a function of the assumed nucleon spin fluctuation rate T„. A perturbative calculation 
yields Tu/T ~ 20-50 for a SN core, but the true value is not known. Still, spin fluctuations 
(dynamical correlations) certainly have a significant impact on the neutrino opacities. 




t-J I I I I I I I I I I I I i-U 

2 4 6 8 

Figure 1: Variation of the thermally averaged axial-current neutrino scattering cross section 
on nucleons as a function of the nucleon spin-fluctuation rate. The curve marked "Lorentzian" 
is obtained from a "resummed" dynamical structure function [Q. 
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A controlled determination of the neutrino opacities requires a calculation of the full 
dynamical spin and isospin density structure function of the correlated nuclear medium at 
finite temperature, a task which has not been accomplished in spite of the recent progress 
in shell-model Monte Carlo calculations of nuclear structure 0. Currently the most reliable 
information of SN neutrino opacities arises from the SN 1987A neutrino pulse duration 
It indicates that the true opacities cannot be too different from the "naive" ones, but this 
conclusion is based on the few late events. The observation of a high-statistics neutrino light- 
curve from a future galactic SN would go a long way towards an empirical calibration of the 
neutrino opacities in hot nuclear matter! 

The fluctuating nucleon spins allow for a more effective transfer of energy between the 
medium and the neutrinos than is possible by elastic collisions alone. This implies that fi 
and T neutrinos may lose more energy between their "energy sphere" and the surface of the 
SN core than had been thought previously, causing their spectra to become more similar to 

|6|, p. If the magnitude of this differential effect is as large as estimated in Refs. |^ 
a number of interesting consequences would obtain, notably for scenarios of SN neutrino 
oscillations which depend on the spectral differences between different flavors. 
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Introduction 

It has been recognized for some years that a neutrino propagating in a plasma acquires an 
induced charge Q. This charge is due to the forward scattering interactions between the 
neutrino and the electrons in the plasma. As these scatterings are electroweak interactions, 
the induced charge is rather small, for example, 10~^^ C, for the dense plasma near the core of 
a type II supernova (SN). Though weak, this charge allows neutrinos in a plasma to undergo 
processes which would usually be restricted to charged particles, such as Cherenkov emission 
or absorption of a photon. 

The electromagnetic interactions of a neutrino in a plasma have been considered in a 
variety of physical and astrophysical scenarios (for a review, see Q). Some recent work by 
Bingham et al. Q proposed a "neutrino beam instability" where an intense flux of beamed 
neutrinos propagating through a plasma leads to the production of an exponentially growing 
number of photons in the plasma through stimulated Cherenkov emission. This form of in- 
stability, caused by electron or photon beams, is well known in plasma physics. The proposed 
application of the neutrino process was as a reheating mechanism for the plasma behind the 
stalled shock of a type II SN. While such an instability is possible in principle, it has recently 
been shown |Q] that this does not occur in type II SNe. 

More recently, Tsytovich et al. Q have proposed an alternative mechanism whereby the 
neutrinos from a type II SN may diffuse slightly in momentum space by propagating through 
a pre-existing saturated thermal distribution of photons. The diffusion mechanism, known as 
quasilinear diffusion, is based on the averaged effect of the individual interactions that occur 
between the photons and the neutrinos. Again, the analogous effect involving electrons is well 
known is plasma physics In their initial calculation, Tsytovich et al. obtained a timescale 
for angular diffusion of the neutrinos from the core of the SN of Tang ~ 10~^ s, this corresponds 
to a scattering length of approximately 30 km, independent of energy, which would be of great 
interest for neutrino transport near the shock of a type II SN. The calculation reported here 
represents a more rigorous calculation of this process and application of the results to a model 
calculation of the plasma properties of a type II SN. It is concluded that this process is only 
likely to be of importance to low energy neutrinos (below 10 keV) and is unlikely to have any 
bearing on the explosion of a type II SN. 

Quasilinear Diffusion Rate 

Given the nature of the plasma behind the shock of a type II SN, it is reasonable to assume a 
high level of plasma turbulence. Within the weak turbulence approximation, this turbulence 
is represented by a distribution of longitudinal photons in the plasma. The strongest level of 
plasma turbulence allowed would be where the energy density associated with the longitudinal 
photons is equal to the energy density associated with the thermal motion of the plasma. 
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Assuming this maximum turbulence gives a natural limit on the strength of the plasma- 
neutrino interaction. 

One may describe the photons in the plasma and the neutrinos from the SN core through 
distributions functions, A^(k), and /(q), respectively. Here, k = {k,9,(l)') denotes the photon 
momentum, and q = (g, a, (p) the neutrino momentum, both given in spherical polar coordi- 
nates. The neutrino distribution evolves through the absorption and emission of longitudinal 
photons, a process which occurs with probability per unit time denoted by u)(q, k). Given 
these definitions the equation governing the time evolution of the neutrino distribution is 
(for the electron equivalent, see @) 



where 
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In equation (|l]), only the largest component of the quasilinear equation has been kept. 
Equation (0) may be evaluated for a saturated thermal distribution of waves and one obtains 



Daa{q) 



IGvrmeC-^ (qc)'^ ' 



(5) 



where Wl is the energy density in the waves, and ujp is the plasma frequency. 

To obtain a simple estimate of the importance of this effect, a few simplifying assumptions 
about the nature of the neutrino distribution are made. It is assumed that the distribution 
is separable, and that the angular distribution is given such that 



with 
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Now, assuming that the angle subtended by the core of the SN on the sky is small, a <C 1, 
one has 
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Thus, 
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where Tang represents the timescale on which the angular scale of the neutrino distribution 
grows through one e-folding. 
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Assuming a saturated thermal distribution, and Wl ~ NgmeV^^ and that VTe ~ c/2 
leads to 

Tang 167r p'^c'^ ' 

which may be parameterized, 

Figure [l| shows this relation plotted for a neutrino spectrum with temperature of 1 MeV 
propagating through the interior regions of a model of a 25 solar mass star which was followed 
through core collapse and bounce by S. Bruenn. Note that the small-angle approximation 
made in producing the angular diffusion timescale break down close to the SN core, thus the 
upper values of this plot are overestimated. Near the shock, at around 10^ cm, the scattering 
timescale is of the order of tens of milliseconds, which is several orders of magnitude too 
slow to be of importance. However, given the inverse square dependence on the energy of 
the neutrino, low energy neutrinos could be scattered quite effectively. Unfortunately, this 
does not represent a reheating mechanism for the post shock plasma, as a thermal neutrino 
distribution of 3 MeV (a realistic temperature) has 99% of its energy in neutrinos above 
2.5MeV, and only a tiny fraction at energies below lOkeV, where the neutrinos are scattered. 

In conclusion it is noted that, while an intrinsically interesting physical effect, the scatter- 
ing of neutrinos off plasma turbulence in a type II supernova is not strong enough to modify 
the bulk of the distribution of neutrinos in the plasma surrounding the SN core. 



hverse scattering timescale vs. r 
10^1 ' ' ' 




10^ 10^ 10^ 10^ 

radius (cm) 



Figure 1: Inverse timescale for angular scattering for 1 MeV neutrino. Calculation of collapse 
and bounce by S. Bruenn, from a 15 solar mass progenitor, t= 0.225. 
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Introduction 

The idea that a large-scale magnetic field can be responsible for a significant neutrino emission 
asymmetry in a supernova explosion goes back to Chugai Q in 1984. Since then there 
have been a number of suggestions of how to implement this idea in more detail. Lately, in 
particular after the Ringberg workshop, several papers appeared dealing with how to calculate 
neutrino cross sections in a magnetic field more accurately and with more realistic neutron 
star parameters |3| . The various approaches are either to use the exact Landau levels for 
the charged particles which is necessary for very strong fields, or to include the field effects 
only through the polarization of the medium. The disadvantage with the full Landau level 
approach is that the number of Landau levels that have to be included grows quadratically 
with the Fermi energy. If the field is S ~ 1 (MeV)^ ~ 1.6x10^^ Gauss and Fermi energy 
Ep = ^/ig — m-g ~ 100 MeV approximately Ep/2eB ~ 15000 Landau levels are filled, and 
with all the transition matrix elements that have to be calculated the problem becomes a 
numerical challenge. For weak fields, where the strength should be compared with Ep rather 
than ml, it seems logical to rely on a weak-field approximation. 

To explain the observed peculiar velocities of neutron stars by asymmetric neutrino emis- 
sion an asymmetry of a few percent is needed. It is found here, as in other recent papers 
1^], that the asymmetry in the damping rate is itself much too small for fields of the order 
of 10^^-10^^ Gauss. On the other hand, it is certain that the final asymmetry cannot be 
esitmated simply by the asymmetry in the damping rate. One step towards a more complete 
analysis was taken in where a cumulative parity violation was suggested to enhance the 
asymmetry. For a reliable estimate it would be desirable to run a full magneto-hydrodynamics 
simulation. 

In this presentation I shall discuss the asymmetry from the electron polarization in 
neutrino-electron scattering and URCA processes. In addition one will have to add the asym- 
metry from proton and neutron polarization in URCA processes which may be as important 
as the electron polarization. The reason is that polarization of the degenerate electron gas is 
suppressed by the large chemical potential even though the electron magnetic moment is much 
larger than the nuclear magnetic moment. Furthermore, I am here using free propagators for 
the nucleons which may have to be improved to account for rapid spin-fiip processes Q. 

Neutrino Damping and the Imaginary Part of the Self-Energy 

The simplest formalism for computing damping coefficients in thermal field theory is, in my 
opinion, from the imaginary part of the self-energy The relation to be used for a 

neutrino with momentum Q is 



r(Q) 
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The damping coefficient consists of two parts F = Tabs + ^cm corresponding to processes tliat 
decrease and increase tfie number of neutrinos in the state {^{Q)). The precise meaning of 
r as defined in Eq. (|l]) becomes clear from the master equation for the neutrino distribution 
function 

^^ = -rabs/.(g) + rcre(i-/.(g)) , (2) 

at 

from which one can show that the neutrino equihbration rate is exactly F. For bosons the 
equilibration rate would have been F = F^bs ~ Tcre which also comes out automatically from 
the imaginary part of the self-energy. For URCA processes Fabs is exactly the damping from 
absorption of v while for I'-e processes Fabs includes both annihilation with V and scattering 
with e, i.e. anything that removes a neutrino from {^{Q)). 
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Figure 1: The two-loop diagram that contains all tree- level v-e processes in its imaginary 
part. 

Let us start with the i/-e-scattering process; the rules for calculating F are given in and 
correspond to a cut setting-sun diagram (Fig. 1): 

(F|ImS|i/) = — ^ J — ^^-yg— p(go,A;o,Po;^e,/^i^)(2vr)^ 

6{K' - ml)5{{K - Pf - ml)5{{Q - Pf){v\^^{l - f )7^(1 - 75)!^) 

tr b'^igv - 9A1^)[W + me)Y{gv - gAi^){W- F + me)] , (3) 

where the distribution functions for initial/final state factors enter through 

P = (^0{qo-po)-U{qo-Po)^I^O{ko)-fe{ko)^l^-0{-ko+po) + fe{ko-po)^ 
+ (-0{-qo+Po) + Uqo-Po)) [-0{ko) + feiko)^ (-0{ko-po) - /e(fco -Po)) , (4) 

where f{po) = (exp[/3(|po| — sign(po)^) + The general rules in ^ for calculating F 

include all possible processes that change the number of neutrinos in a certain state, and the 
thermal factors for scattering/absorption/creation are all automatically included. It is not 
difficult to separate out Fabs and F^e from F. Absorption corresponds to the first term of p in 
Eq. (0) and creation is the second one. It should be noted that F represents the integrated rate 
at which the incoming neutrino changes state. In a Boltzmann equation the differential rate 
of scattering into another momentum would be needed, and that is easily obtained by simply 
omitting the integration over the neutrino inside the loop. It is also easy to separate Fabs and 
Fere by noting that detailed balance requires Fabs/rcre = (1 — fu)/fu = exp[/3(go — Pu)] (for 
go > 0) in equilibrium, which is automatically satisfied with Eq. (Q). 
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So far the damping rate corresponds to the zero field case, but the generalization to a 
background magnetic field is straightforward after identifying the different factors in Eq. (|3|). 
The 5-functions come from the imaginary part of the usual electron and neutrino propagators 
(or rather from the thermal part of the real-time propagators). We now have to use the 
electron propagator in an external field instead, and here it is convenient to use Schwinger's 
formulation 



as exp 



p2 _ 77j,2 _|_ Jq cos{eBs 



. , „9 tan(ei?s) „o o . ^ 

" eBs 



^— ieBsa:, 



where a • fty = aobb — a-zhz and a-b± = —axbx — a-yby. This expression is quite difficult to use in 
general but the complicated parts, cos(ei?s) and tan{eBs)/eBs, have no linear dependence 
on eB and can therefore be approximated by 1 for weak fields. The s-integral can then be 
performed and after taking the imaginary part we end up with the replacement rule 

(f + me)5(p2 _ ml) ^ + me)eScT,sign(p2 - ml)5{(P^ - mlf - (eBf) . (6) 

To obtain this rule I have kept the full eB dependence in the exponentials in Eq. (^) since 
they are necessary for the IR convergence of the energy integrals. In the final result, after 
performing all integrals, the damping rate has a very clear linear dependence on eB for 
weak fields. 

To compute the damping from absorption and creation of neutrinos in URCA processes 
one can use the same formalism but with neutron and proton propagators in the loop. In the 
present case absorption is dominted hy v + n^p + e and creation by the inverse process. 
The nucleons are also essentially nonrelativistic which simplifies the matrix elements. On the 
other hand, since the general expression comes out for free and the trace is easily evaluated 
with a symbolic algebraic program, we might as well keep the full relativistic expression and 
all possible channels when we perform the numerical integrals. 



Results 

The physical conditions in a supernova vary considerably during the short time of the explo- 
sion. I have as an illustration evaluated the damping asymmetry for two more or less typical 
conditions. The first case has a high electron chemical potential (/Xg = 200 MeV) and neutrino 
energy qq = 100 MeV, and the other one has fie = 50 MeV and qo = 30 MeV. In both cases 
I used T = 10 MeV and Ye = 0.3. Local thermal equilibrium was assumed and the proton 
chemical potential was determined from charge neutrality. These two cases correspond to 
the densities 2x10^^^ g/cm^ and 4x10^^ g/cm^. For simplicity the incoming neutrino has its 
momentum parallel to the magnetic field (the linear term in eB is of course antisymmetric 
in Qz). Writing the damping factor as F = F^'^^ -|- eBT^^\ I found the relative asymmetry for 
URCA and u-e processes to be (i?i4 is the field strength in units of 10^^ Gauss) 

~ 9.8xl0~'^Si4 , (7) 
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for He = 200 MeV and 



eBT 



(1) 

URCA 



URCA 



1.0x10-^^14 



eBT 



(1) 



7.3x10-^5 



14 



(8) 



for /ig = 50 MeV. These relative asymetries are as large as the ones calculated in Q and one 
can expect the nucleon polarization to give a comparabel contribution. I therefore conclude, 
based on the results in that neutrino absorption is an important process for the neutrino 
opacity (see Fig. 4 in [Q) and that the asymmetry in URCA processes needs to be taken into 
account in any reliable estimate of the total neutrino propagation asymmetry in a neutron 
star. However, it is also my impression that the final word has not yet been said about the 
influence from collective behavior of nucleons on the neutrino propagation. 
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Abstract 

We discuss the Cherenkov process u ^ in the presence of a homogeneous magnetic field. 
The neutrinos are taken to be massless with only standard-model couplings. The magnetic 
field fulfills the dual purpose of inducing an effective neutrino-photon vertex and of modifying 
the photon dispersion relation such that the Cherenkov condition uj < |k| is fulfilled. For a 
field strength -Bcrit = m^/e = 4.41 x 10^^ Gauss and for E = 2me the Cherenkov rate is about 
6 X 10"^^ s"^. 

In many astrophysical environments the absorption, emission, or scattering of neutrinos occurs 
in dense media or in the presence of strong magnetic fields |jl|. Of particular conceptual 
interest are those reactions which have no counterpart in vacuum, notably the decay 7 — > Pi/ 
and the Cherenkov process ly — > lyj. These reactions do not occur in vacuum because they are 
kinematically forbidden and because neutrinos do not couple to photons. In the presence of a 
medium or i?-field, neutrinos acquire an effective coupling to photons by virtue of intermediate 
charged particles. In addition, media or external fields modify the dispersion relations of all 
particles so that phase space is opened for neutrino-photon reactions of the type 1 — > 2 -|- 3. 

If neutrinos are exactly massless as we will always assume, and if medium-induced modifi- 
cations of their dispersion relation can be neglected, the Cherenkov decay — > z^7 is kinemat- 
ically possible whenever the photon four momentum k = {uj, k) is space-like, i.e. — > 0. 
Often the dispersion relation is expressed by |k| = noj in terms of the refractive index n. In 
this language the Cherenkov decay is kinematically possible whenever n > 1. 

Around pulsars field strengths around the critical value i3crit = m'^/e = 4.41 x 10^^ Gauss 
are encountered. The Cherenkov condition is satisfied for significant ranges of photon frequen- 
cies. In addition, the magnetic field itself causes an effective i/-7-vertex by standard-model 
neutrino couplings to virtual electrons and positrons. Therefore, we study the Cherenkov 
effect entirely within the particle-physics standard model. This process has been calculated 
earlier in but we do not agree with their results. 

Our work is closely related to a recent series of papers who studied the neutrino 
radiative decay v i^'j in the presence of magnetic fields. Our work is also related to the 
process of photon splitting that may occur in magnetic fields as discussed, for example, in 
Refs. ii. 

Photons couple to neutrinos by the amplitudes shown in Figs. 1(a) and (b). We limit our 
discussion to field strengths not very much larger than Bent = m1/e. Therefore, we keep only 
electrons in the loop. Moreover, we are interested in neutrino energies very much smaller than 
the W- and Z-boson masses, allowing us to use the limit of infinitely heavy gauge bosons and 
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Figure 1: Neutrino-photon coupling in an external magnetic field. The double line represents 
the electron propagator in the presence of a i?-field. (a) Z- ^-mixing, (b) Penguin diagram 
(only for Vf.)- (c) Effective coupling in the limit of infinite gauge-boson masses. 

thus an effective four-fermion interaction [Fig. 1(c)]. The matrix element has the form 

M = -^Ze,v^,{l - 75)z^ (ffyn^'^ - 5Anf ), (1) 
v2e 

where e is the photon polarization vector and Z its wave-function renormalization factor. 
For the physical circumstances of interest to us, the photon refractive index will be very 
close to unity so that we will be able to use the vacuum approximation Z = 1. Further, 
gv = 2siv? 9w + |, 9A = h for i^e, Qv = 2sin^ Ow - \, and qa = -\ for v^^r- 
Following Refs. |, g we find 



where F'^'^ = ^e^'^^'^Fpa- and F12 = —F21 = B. The || and _L decomposition of the metric is 
5|j = diag(— , 0, 0, -|-) and g^_ = g — g\\ = diag(0, 0) and k is the four momentum of the 
photon. A'o, A^j_,A''||, C_\_ and Cy are functions of B, fcy and k\. They are real for uj < 2me, 
i.e. below the pair-production threshold. 

Four-momentum conservation constrains the photon emission angle to have the value 



cos 9 = — 

n 



l + („2_i)^ 



(3) 



where 9 is the angle between the emitted photon and incoming neutrino. It turns out that for 
all situations of practical interest we have |n — 1| <C 1 @, This reveals that the outgoing 
photon propagates parallel to the original neutrino direction. 

It is easy to see that the parity-conserving part of the effective vertex (11'^'') is proportional 
to the small parameter (n — 1)^ <C 1 while the parity- violating part {H'^'^) is not. It is 
interesting to compare this finding with the standard plasma decay process ^ ^ vv which 
is dominated by the Yl^^ . Therefore, in the approximation svc? 9^/ = ^ only the electron 
flavor contributes to plasmon decay. Here the Cherenkov rate is equal for (anti)neutrinos of 
all flavors. 

We consider at first neutrino energies below the pair-production threshold E < 2me. For 
UJ < 2me the photon refractive index always obeys the Cherenkov condition n > 1 [§, |^. 
Further, it turns out that in the range < u < 2me the functions Cm ,Cj_ depend only weakly 
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on uj so that they are weU approximated by their value at u; = 0. For neutrinos which 
propagate perpendicular to the magnetic field, the Cherenkov emission rate can be written 
in the form 

^ ' > h{B) = 2.0xl0-^s-^ - — - — h{B), (4) 



135(47r)4 VBcrit/ ' ' \2mJ V^crit/ 

where 

/j^^) = / (4/25) (5/Bcrit)' for 5 « Sent, 
^ ^ I 1 for S » Bent. 

Turning next to the case E > 2me we note that in the presence of a magnetic field the 
electron and positron wavefunctions are Landau states so that the process v — > i^e+e" becomes 
kinematically allowed. Therefore, neutrinos with such large energies will lose energy primarily 
by pair production rather than by Cherenkov radiation; for recent calculations see |[To[] . 

The strongest magnetic fields known in nature are near pulsars. However, they have a 
spatial extent of only tens of kilometers. Therefore, even if the field strength is as large as the 
critical one, most neutrinos escaping from the pulsar or passing through its magnetosphere 
will not emit Cherenkov photons. Thus, the magnetosphere of a pulsar is quite transparent 
to neutrinos as one might have expected. 
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In astrophysical plasmas the dispersion relation of photons is usually dominated by the elec- 
tronic plasma effect. It was recently claimed Q, however, that in a supernova (SN) core the 
dominant contribution is caused by the photon interaction with the nucleon magnetic mo- 
ments. This contribution can have the opposite sign of the plasma term so that the photon 
four-momentum K could become space-like, allowing the Cherenkov processes ^ u and 
^ 'ji'. Because of a numerical error the results of Ref. are far too large but the 
correct effect is still large enough to call for a closer investigation. 

Therefore, we will estimate the different possible contributions to the refractive index, 
which is taken to be a real number, 

nrefr = ^e/i = /c/w. (1) 

It relates the wavenumber k and frequency uj of an electromagnetic excitation and is given in 
terms of the dielectric permittivity e and the magnetic permeability Q . We may also write 
the dispersion relation in the form uP' — 1? = rrigg, where m^g can be positive or negative. 

We first assume that the medium constituents (protons, neutrons, electrons) interact with 
each other only electromagnetically. This leaves us with the usual plasma effect which is much 
larger for the electrons than for the protons. It is approximately described by n^g£j. 1— Wp/w^, 
where = ^TraUp/mp for the protons (nondegenerate, nonrelativistic) and uJ^ = (4a/37r);Ug 
for the electrons (degenerate, relativistic) with //e the electron chemical potential. 

If the nuclei are not completely dissociated, nuclear bound-free transitions might be sig- 
nificant. However, at least for the simple example of deuteron states this is not the case. For 
a medium consisting of deuterons alone we employ the usual relation between refractive index 
and forward-scattering amplitude /o 0] n^eix = 1 + (27r/u;^) Re/o(u;). From the values for 
Re/o of Ref. |5| we conclude that the contribution from deuteron-like states never exceeds the 
single-proton plasma effect and is thus negligible. 

The main point of our work is to calculate the contribution from correlated spins. We 
study the simplest model, a medium consisting purely of neutrons. In the nonrelativistic 
case they do not respond to an applied electric field so that e = 1. We write the magnetic 
permeability /x = 1 + x in terms of the magnetic susceptibility X = x' + i'x" which implies 
'^refr ~^ ~ x'i^j k). The fluctuation dissipation theorem Q establishes a relationship between 
the absorptive part x" ^ind spin fluctuations. The neutron spin fluctuations can be described 
by the dynamical spin-density structure function Sa{io,k) 0. We use a simple model for its 
actual shape § in terms of the spin fluctuation rate Fg- which can be estimated from the 
nucleon interaction potential. We calculate the refractive index (or rather x') from x" via a 
Kramers-Kronig relation and express it as an effective mass m^g = (1 — n'^^i^)uj'^. 

In Fig. [I we show the different contributions to m^g relative to the one by the electrons 
for typical conditions of a SN core (p = 7 x lO^^gcm"^, electron to baryon ratio 1^ = 0.3, 
T = 40MeV, Tu = SOT). All the approximations made in our calculation work in the same 
direction of overestimating the magnetic-moment contribution. 
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Figure 1: The contributions of electrons, protons and correlated neutrons for p = 7 x 
lOi^gcm-3, Ye = 0.3, T = 40MeV and = SOT. 

The contribution of the correlated spins has the opposite sign relative to the plasma terms 
only in a small to range and thus cannot change the sign of the total rn^g. Thus the refractive 
index is always less than 1 and Cherenkov effects remain forbidden. 
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Abstract 

The transients following GRB970228 and GRB970508 showed that these (and probably all) 
GRBs are cosmological. However, the host galaxies expected to be associated with these 
and other bursts are largely absent, indicating that either bursts are further than expected 
or the host galaxies are under luminous. This apparent discrepancy does not invalidate the 
cosmological hypothesis, but instead host galaxy observations can test more sophisticated 
models. 



Absence of Expected Host Galaxies 

Observations of the optical transients (OTs) from GRB970228 pl and GRB970508 H have 



finally provided the smoking gun that bursts are cosmological. In most cosmological models 
bursts occur in host galaxies: are these galaxies present, and conversely, what can we learn 
from them? Underlying any confrontation of theory and data must be a well defined model. 
Here we show that the host galaxy observations are not consistent with the expectations 
of the simplest cosmological model, and that these observations can be used to test more 
sophisticated models. 

In the simplest ("minimal") cosmological model the distance scale is derived from the 
intensity distribution logN-logP assuming bursts are standard candles which do not evolve 
in rate or intensity. Bursts occur in normal galaxies at a rate proportional to a galaxy's 
luminosity. This model predicts the host galaxy distribution for a given burst. Are the 
expected host galaxies present? 



For GRB970228 an underlying extended object was found [|19|], but its redshift and nature 
have not been established. If the observed "fuzz" is indeed a galaxy at z ~ 1/4, it is ~5 
magnitudes fainter than expected for a galaxy at this redshift. For GRB970508 no obvious 
underlying galaxy was observed flj] and the nearest extended objects have separations of 
several arcseconds, but spectroscopy with the Keck telescopes [|^] led to the discovery of 
absorption and emission lines giving a GRB redshift of z > 0.835. The HST magnitude limit 
R\im ~ 25.5 |l^] for a galaxy coincident with the transient again suggests a host galaxy fainter 
than expected. Similar conclusions follow from the inspection of IPN error boxes ll],[l^, |l3l, 



but see also M, 20 1. This absence of sufficiently bright host galaxies is often called the "no- 



host" problem, which is a misnomer. The point simply is that if galaxies such as the Milky 
Way provide the hosts to most bursters, and if their redshifts are less than unity, as predicted 
by the minimal model, we expect to find bright galaxies inside a large fraction of the smallest 
IPN error boxes. 

To demonstrate this quantitatively, consider the apparent magnitude of a typical host 
galaxy, which we assume has M^{B) = —20 (approximately the absolute magnitude of an L^, 
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galaxy — see discussion below). Using Peebles' notation ||T3[, the apparent magnitude is 

m = 42.38 + M + 5log[y{z){l + z)] + K{z) + E{z)+A{n,z)+x{z) , (1) 

where K(z) is the usual K-correction, E(z) corrects for the possible evolution of the host 
galaxy's spectrum, A is the sum of Galactic foreground (position dependent) and intergalactic 
extinction, and x{^) represents any corrections that apply in hierarchical galaxy formation 
scenarios, where galaxies are assembled through the merger of star forming subunits. The 
commonly found term 51og(/i) is already absorbed in eq. (1). Neglecting potentially large 
corrections from the K, E, A, and xi^) terms, a host like the Milky Way with M ~ —20 
would have an apparent magnitude m ~ 22 for redshifts of order unity. Several small IPN 
error boxes have no galaxy of this magnitude or brighter. Our simplified treatment agrees 
with Schaefer's conclusion [16, 17] that typical galaxies at the calculated burst distance are 



absent from burst error boxes. 

Thus bursts are further than predicted from the logN-logP distribution without evolution, 
or they occur in underluminous galaxies; an extreme limit of the latter alternative is that 
bursts do not occur in galaxies. 

Host Galaxies as a Probe of Cosmological Models 

The search for host galaxies is a powerful test of cosmological burst models. From the two 
above mentioned OTs we conclude that GRBs are cosmological, but the observations have not 
fixed the distance scale quantitatively, nor have they determined the energy source. While the 
x-ray, optical, and radio lightcurves (for GRB970508 only) are consistent with the predictions 
of the basic "fireball afterglow" picture, the fireball's central engine could be the merger of a 
neutron star binary, the collapse of a massive, rotating star, or the jet produced by accretion 
onto a massive black hole residing at the center of an otherwise normal galaxy. 

The host galaxies found within burst error boxes are a powerful discriminant between 
different models for the burst energy source. Note that the region within which a host would 
be acceptable surrounding the sub-arcsecond localizations of an OT is effectively the error box 
for the host galaxy. Almost all models assume that bursts are associated with galaxies; the 
issue is the relationship between the burst and the host. In models such as the momentary 
activation of a dormant massive black hole the burst rate per galaxy is constant. On the 
other hand bursts are an endpoint of stellar evolution in most models, and therefore to first 
order we expect the burst rate per galaxy in these models to be proportional to the galaxy's 
mass and thus luminosity. These two model classes have different host galaxy luminosity 
functions ^{M) with different average values of M (the absolute magnitude). In the first 
case, ip{M) is proportional to the normal galaxy luminosity function, while in the second case 
■ip{M) is proportional to the normal galaxy luminosity function weighted by the luminosity 
L oc 10~^'^^. We approximate the normal galaxy luminosity function with the Schechter 
function: 

cI>(M) = Acl00-^(*^*-*^)("+i) [exp(-10°-4(*^*-^^))l (2) 



where k is the normalization, a is the slope of the faint end, and M^, is the absolute magnitude 
of an galaxy. Here we use a = —1. In the B band M^{B) = —19.53, which corresponds 
to L^{B) = 1.8 X lO^^L© ~ 3 X lO^^ Lq{B) hj^. In Fi gure 1 we show the cumulative 
distributions for the host galaxy magnitudes for the two model classes. As can be seen, the 
average host galaxy magnitude (i.e., at 0.5) differs by ~ 1.75 magnitudes. 
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Figure 1: The cumulative distribution of host galaxy magnitudes if their luminosity function 
is weighted by number or luminosity. If GRBs trace light, the 50% point is near M^,. If it is 
proportional to galaxy number a typical host would be ~ 2 magnitudes fainter. 



However, we can make better predictions about the host galaxies in cosmological models 
where bursts are a stellar endpoint. In such models, the burst rate should be a function of 
the star formation rate (SFR). If there is a substantial delay (e.g., of order a billion years 
or more) between the GRB event and the star forming activity that created the progenitor, 
then the burst rate integrates over a galaxy's SFR, and we would not expect the host galaxy 
to display the signatures of recent star formation. Furthermore, if the progenitor is given a 
large velocity, then it may travel a large distance from the host galaxy before bursting, and 
it may become impossible to associate a galaxy with the burst. 

In many models the burst occurs shortly after its progenitor star forms (e.g., within a 
hundred million years or less). We would then expect that on average bursts would occur in 
galaxies showing evidence of recent star formation. The burst rate should be proportional to 
the SFR, both for individual galaxies and for a given cosmological epoch. 

In particular, the burst rate and the SFR should have the same history, as was recently 
considered by several groups [|l2|, 15, |l^, 21 1. Extensive redshift surveys and data from the 
Hubble Deep Field have reliably determined the cosmic star formation history to z ~ 5 
1^, ^, 1^, ^, ^- The data clearly suggest a rapid increase in the comoving SFR density with 
increasing redshift, SFR oc (1 -|- z)^, reaching a peak rate (at z ~ 1.5) about 10-20 times 
higher than the present-day rate, and decreasing slowly to the present value by z ~ 5. This 
evolution function, r/(z), enters the differential rate vs. (bolometric) peak flux 



dpR oc ^(^) (1 + ^)-3 [(1 + z)d,y{z) + y{z)] 



-1 



(3) 



where 'E{z) and y{z) are defined in ||13[. For small redshifts the logarithmic slope of this 
function is Euclidean, i.e. —5/2. The solid curve of Figure 2 shows the effects of geometry 
(bending of logN-logP) and the dashed curve demonstrates how r]{z) compensates for the 
geometry out to the redshift at which the cosmic SFR peaks. At larger redshifts the effects 
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Figure 2: The slope of the differential GRB brightness distribution. Non-evolving sources 
(solid curve) quickly show significant deviation from the Euclidean value —5/2 with increasing 
redshift. If the GRB rate is proportional to the SFR (see text) the apparent Euclidean slope 
extends to z ~ 1 (dashed curve). For greater z the geometry of the universe together with a 
now decreasing burst rate cause the slope to deviate rapidly from —5/2. 



of geometry and decreasing SFR then combine and the slope flattens quickly. Comparison 
with BATSE data suggests that this SFR model deviates from the pseudo-Euclidean slope 



too abruptly. While several studies [|15|, [l^, ^ report that the observed SFR generates a 



brightness distribution consistent with BATSE data, our findings support the different result 



of Petrosian & Lloyd [12|, who suggest that other evolutionary effects must be present in 
addition to the density evolution described by 7?(z). While a good fit to the data requires a 
more sophisticated model of source evolution, the basic message is likely to be the same: the 
logN-logP distribution does not exclude GRB redshifts much greater than unity. 

Therefore, the absence of the host galaxies predicted by the "minimal" cosmological model 
does not call the cosmological origin of bursts into question. Instead, host galaxy observations 
will teach us where bursts occur. 
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Introduction 

Thirty years after their discovery and over 140 models later, we are still struggling to under- 
stand the nature of gamma-ray bursts. There has been considerable evolution in our views 
during the last few years, especially the last six months, due largely to the recent observa- 
tional data summarized elsewhere in these proceedings by Dieter Hartmann. A consensus, by 
no means universal, has developed that gamma-ray bursts occur at cosmological distances, at 
least Gpc's, far enough away for their distribution to be very isotropic, yet sampled with suf- 
ficient sensitivity that the edge of the distribution is being seen. This requires burst energies 
of over 10^^ erg (times a beaming factor that could be <C 1) in an event that occurs about 
1000 times annually, has median energy around 600 keV, and typical duration from a small 
fraction of a second to hundreds of seconds (biniodal with means 0.3 and 20 s). The large 
77-opacity above the pair threshold requires that a burst originate in a large volume, perhaps 
10^^ cm or more in radius, and the short duration at this large size requires large relativistic 
factors, r <; 100. According to the current paradigm (e.g., Meszaros k. Rees 1993), this is 
accomplished by having a small amount of matter, ^ 10~^ M,:,, move with the requisite V and 
produce the gamma-rays either by internal shocks or by interaction with the circum-source 
material, which may be the interstellar medium. 

The quest in current gamma may burst "models" — there is still no complete example 
of a natural event modeled thru to the production of gamma-rays — has thus turned to an 
intermediate goal of delivering to some small amount of matter the necessary tremendous 
energy to make it move so close to the speed of light. There have been many attempts. Most 
can be ruled out. 

Models that don't work 

a) All Forms of Accretion Induced Collapse in White Dwarf Stars 

In its simplest form (e.g., Dar et al. 1992), an unmagnetized white dwarf collapses to a neutron 
star and produces a pair fireball by neutrino annihilation above its surface. The problem is 
that the same neutrinos blow a wind from the neutron star of about 0.01 M,:, s^^. The energy 
from neutrino annihilation is deposited deep within this wind and adiabatic expansion leads 
to a very cool photosphere. No gamma-rays are produced and the mass ejected has V much 
less than one. 

A variation on this model is the rotating magnetic white dwarf that collapses to a very 
energetic pulsar (Usov 1988). The pulsar converts a portion of its rotational energy into 
gamma-rays and pairs — another pair fireball model. The problem here is again the neutrino 
driven wind, which since it has a much greater energy density than the magnetic field, blows 
it away before any burst can be generated. 
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b) Phase Transitions in Neutron Stars 

An example might be the transition of a neutron star to a strange star after accreting a critical 
mass in a binary system. If the phase transition results in a more compact core forming in 
less than a sound crossing time, the bounce forms a shock that propagates to the edge of 
the neutron star with energy that may be in excess of 10^^ erg. Steepening of the shock 
in the density gradient near neutron star surface might concentrate this energy in a small 
amount of mass. The problem here, besides the obvious issue of whether such supersonic 
phase transitions ever occur, is that the shock hydrodynamics has been calculated and shown 
inadequate to produce sufficient high energy ejecta. Fryer k. Woosley (1998) show that, even 
in an optimal model, less than 10^^ erg is concentrated in matter having V ^ 40. Most of 
the shock energy goes into heating and expanding the neutron star and what remains, about 
10^^ erg, ejects not IQ-^ Mg, but 0.01 Mq. Most of this matter is at low V. 

c) Supermassive Black Holes 

Such objects are thought to power AGN's including gamma-ray blazars; why not gamma-ray 
bursts (Carter 1992)? One major problem is the long orbit time near the event horizon, 
~5000 (M/10^ Mq) s, implies an event duration much longer than observed in gamma-ray 
bursts. We presume that matter must orbit the hole at least a few times before building up 
large magnetic fields or other means of dissipating energy and forming jets. 

Models That Might Work 

In general, the models that might work all involve accretion at a very high rate into a newly 
formed black hole whose mass is between two and a few hundred solar masses. The obser- 
vational counterpart to neutron star formation is a supernova. Perhaps a gamma-ray burst 
is the signal that a black hole has been born. There are four ways of setting up the requi- 
site conditions: a) merging neutron stars forming a black hole with a residual disk of a few 
hundredths of a solar mass; b) neutron star, black hole mergers which might leave a disk of 
a few tenths Mq; c) collapsars in which the core of a massive star collapses to a black hole 
and accretes the rest of the star (several solar masses endowed with sufficient rotation) ; and 

d) black hole, white dwarf mergers. All of these events should occur chiefly where there are 
stars, i.e., in galaxies. It would be difficult for binaries consisting of compact stars and BH's 
to remain bound when kicked out of the galaxy. WD-BH mergers might occur not only in 
ordinary stellar systems, but in binaries formed in dense galactic nuclei (Quinlan and Shapiro 
1990). Collapsars in particular occur only in star forming regions. 

The accretion rate depends on the viscosity of the disk, or in the case of the coUapsar 
model upon the time required for the star to collapse to the disk (446/p^/^ s). Because the 
star has a range of densities, accretion would continue at a diminishing rate for a very long 
time after the gamma-ray burst perhaps providing an enduring x-ray transient. Typically the 
accretion rate in all these models starts at ~ 0.1 Mq s~^ and declines. For such accretion 
rates one expects densities in the inner disk of around 10"^^ g cm^^ (Popham &; Gammie 1997; 
or simply seen by equating irR'^hpvdnit ~ M). At this density the inner disk is marginally 
thick to neutrinos. Typical drift velocities are ~1000 km/s implying a residence time in the 
disk of order 0.1 s (this number is viscosity (a) dependent, and uncertain), so the duration 
of bursts from merging compact objects could be this short (longer burst time scales are 
allowable depending on the circumstellar interaction, but none shorter than the cycle time 
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of the engine itself). A longer time scale would characterize the collapsars where the disk is 
continually fed by stellar collapse. Perhaps the long complex bursts (20 s mean) are due to 
collapsars and the short bursts are merging compact objects. 

For an accretion rate of 0.1 Mq s~^, the energy dissipated, mostly in neutrinos, is 1.8 x 10^^ 
F erg/s where F is 0.06 for Schwarzschild geometry (perhaps appropriate for merging compact 
objects) and 0.4 for extreme Kerr geometry (perhaps appropriate for collapsars). Coupled 
with their short time scale this implies that merging compact objects may give bursts of much 
less total energy than collapsars, but the efficiency for getting the jet out of a coUapsar may 
modulate this considerably (see below). In all models, a portion of the energy goes into MHD 
instabilities in the disk that may result in an energetic wind, or perhaps the rotational energy 
of the black hole is tapped as in some models for AGN. But the simplest physical solution for 
ejecting relativistic matter, if it works, is neutrino annihilation along the rotational axis of the 
hole (Woosley 1993). The efficiency for this is uncertain, but numerical simulations (see the 
paper by Ruffert, these proceedings) suggest an efficiency for converting neutrino luminosity 
into pair fireball energy of Thus one has a total energy for black holes merging with 

neutron stars of order 10^° erg (Ruffert et al. 1997) and a total energy for collapsars that 
might approach 10^^ erg. 

The chief uncertainty in all these models — aside from exactly how the relativistic matter 
reconverts its streaming energy into gamma-rays — is the baryon loading. One expects this 
to be relatively small in neutron star and black hole mergers, but still perhaps more than 
10~^ Mq, depending upon the mingling of the wind from the disk, residual matter from the 
merger, and the relativistic outflow. Detailed calculations remain to be done. 

The problem is much more severe in the coUapsar model where the jet that is formed has 
to penetrate the overlying star — or what is left of it — a little like making the 7-ray burst at 
the center of the sun. Still 10^^ erg is a lot of energy and we also expect that the jet from 
collapsar models might be tightly beamed. This is initially the case because the disk makes a 
transition from very thick (geometrically) to thin at the point where neutrino losses start to 
occur on an accretion time scale. This happens at a definite radius, ~100 km, because of the 
T^ dependence of the neutrino rates (Bob Popham, private communication). Calculations in 
progress by Miiller and colleagues at the MPA will show how effective this jet is at tunneling 
thru the star. Our own recent 2D simulations suggest a solid angle of only about 1% for 
the jet. A lot of energy gets used up ejecting all the mass above about 45 degrees from the 
equator. A supernova is one consequence; a collapsar is not a "failed supernova" after all, but 
this mass ejection at large angles is non-relativistic. Along the axis the jet clears a path for 
more energetic matter to follow, and since the burst continues for many jet-stellar crossing 
times, F rises. Only ~0.01 Mq initially lies in the path of the jet for the expected small 
opening angle. 

The energy of the ejecta and presumably the hardness of the transient that is seen there- 
after thus depends upon viewing angle. Along the axis, high F will overtake low F material 
before the observable event actually commences. The start of the burst occurs when the jet 
sweeps up 1/F times its rest mass, or about 10~^ Mq. This happens after the jet has gone 
about 10^^ cm and interacted with the stellar wind ejected prior to the explosion (about 
10~^ Mq assumed). 

In summary we expect about lO^'^ erg, ^ 0.1% of the total energy dissipated in the disk, 
to go into a jet with F ^ 100 and an opening angle of about 10 degrees. The event rate could 
be 10~^ to 10~^ per bright galaxy. 

Below we show pictures of our first 2D calculation of rotating stellar collapse about 10 s 
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after the collapse began (nb., no explosion, contrary to Bodenheimer & Woosley, 1983). A 
black hole and accretion disk have clearly formed. The assumed angular momentum here was 
about 10^^ cm^ s~^ in the region that formed the disk. Stellar evolution calculations suggest 
that a value of a few times 10^^ may be more appropriate to the stellar mantle, but larger 
values exist farther out. This will merely delay the formation of the disk. The figure on the 
right shows a first crude attempt at running a jet thru this stellar model. The code employed 
lacked relativistic hydrodynamics and the jet was artificially injected at the inner boundary 
with a kinetic energy of 4 x 10^*^ erg and a speed of 2/3 c corresponding to a mass loss rate 
of 0.001 Mq s-^ For now this is just a suggestive picture of what a more realistic calculation 
might give. In a relativistic code, the same energy injection with T ~ 100 would correspond 
to a (rest) mass loss rate 100 times less, comparable to what one needs in a 7-ray burst. 




1.00x10'° 



X (cm) 



Figure 1: Left: Velocity field and density contours of a preliminary rotating collapsar model 
at t = 10 s after collapse. The dashed circle in the center represents a completely absorbing 
inner boundary at r = 5 x 10^ cm. The core has grown to 2.85 Mq. Right: Jet emerging from 
stellar surface at 2 x lO^'' cm. The longest arrow is a velocity of 1.7 x lO'''^ cm s"^. Clearly 
given enough momentum, the jet can penetrate the star. 



The material that is ejected at lower F at angles farther away from the rotational axis may 
also be of considerable interest. For intermediate values of F, say ~1-10, one may still get a 
bright circumstellar interaction, but because it is more massive and must intercept a larger 
fraction of its rest mass before giving up its energy, the display from this interaction may 
be delayed and prolonged. We expect that this is the origin of the optical, x-ray, and radio 
afterglows detected in the recent gamma-ray bursts of 022897 and 050897. As much or more 
energy could be in these displays than in the 7-ray burst itself. The softer radiation would 
also be beamed to a larger part of the sky suggesting the possibility that there are many more 
optical, radio, and x-ray transients from 7-ray burst sources than GRB's themselves. 

Further details of our models and other issues discussed in this short paper will appear in 
a paper to be submitted soon to the ApJ. 
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Abstract 

We have simulated a variety of models, in order to shed light on the accessible physical 
conditions during the mergers, on the amount of matter dynamically ejected during the 
merging, on the timescales of mass accretion by the (forming) black hole, on the conversion of 
energy into neutrino emission, on the amount of energy deposited by i/z^-annihilation, and on 
the baryon loading of the created pair-plasma fireball. To this end, we varied the masses and 
mass ratios as well as the initial spins of the neutron stars, changed the impact parameter 
to consider spiral-in orbits and direct collisions, included a black hole (vacuum sphere) in 
our simulations, and studied the dynamical evolution of the accretion torus around the black 
hole formed after the neutron star merging until a (nearly) stationary state was reached. 
While the neutrino emission during the dynamical phases of the mergings is definitely too 
small to power gamma-ray bursts (GRBs), we find that the masses, lifetimes, and neutrino 
luminosities of the accretion tori have values that might explain short (0(0.1-1 s)) and not 
too powerful (~ 10^^/(47r) erg/(s • sterad)) gamma-ray bursts. 

Summary of Numerical Procedures and Initial Conditions 

The hydro dynamical simulations were done with a code based on the Piecewise Parabolic 
Method (PPM) developed by Colella & Woodward Q. The code is basically Newtonian, but 
contains the terms necessary to describe gravitational wave emission and the corresponding 
back-reaction on the hydrodynamical fiow (Blanchet et al. Q). The terms corresponding to 
the gravitational potential are implemented as source terms in the PPM algorithm. In order to 
describe the thermodynamics of the neutron star matter, we use the equation of state (EOS) of 
Lattimer & Swesty ||^ . Energy loss and changes of the electron abundance due to the emission 
of neutrinos and antineutrinos are taken into account by an elaborate "neutrino leakage 
scheme" Q. The energy source terms contain the production of all types of neutrino pairs by 
thermal processes and of electron neutrinos and antineutrinos also by lepton captures onto 
baryons. Matter is rendered optically thick to neutrinos due to the main opacity producing 
reactions which are neutrino-nucleon scattering and absorption of neutrinos onto baryons. 
More detailed information about the employed numerical procedures can be found in Ruffert 
et al. Q. The following modifications compared to the previously published simulations (Q] 
and Q) were made in addition, (a) The models were computed on multiply nested and refined 
grids to increase locally the spatial resolution while at the same time computing a larger total 
volume. Our grid handling is described in detail in Section 4 of Ruffert 0. (b) An entropy 
equation instead of the equation for the total specific energy (specific internal energy plus 
specific kinetic energy) was used to calculate the temperature of the gas in order to be able 
to determine low temperatures more accurately, (c) The equation of state table was extended 
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to 100 MeV < T < 0.01 MeV and 2.9 x lO^^g/cm^ > p > 5 x lO'^g/cm^ because in extreme 
cases very high temperatures can be reached and the density of ejected gas decreases to low 
values at large distances. 

We started our simulations with two cool (temperatures initially some 0.01 MeV to a few 
MeV) Newtonian neutron stars with baryonic masses between 1.2 Mq and 1.8 Mq (depending 
on the particular model), a radius of approximately 15 km, and an initial center-to-center 
distance of 42-46 km. The distributions of density p and electron fraction Yg were taken from 
one-dimensional models of cold, deleptonized neutron stars in hydrostatic equilibrium. We 
prescribed the orbital velocities of the coalescing neutron stars according to the motions of 
point masses, as computed from the quadrupole formula. The tangential components of the 
velocities of the neutron star centers correspond to Kepler velocities on circular orbits, while 
the radial velocity components result from the emission of gravitational waves that lead to the 
inspiral of the orbiting bodies. For an observer not corotating with the system, the neutron 
stars were given either no spins or the two stars were assumed to be tidally locked or to have 
spins opposite to the direction of the orbital angular momentum 

neutrino radiation, Model D 

3.0 L ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' ' ' I ' . 
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Figure 1: Luminosities of individual neutrino types (fe, i^e, and the sum of all Vx) and of the 
total of all neutrinos as functions of time for the merging of a I.2M0 and a I.SMq neutron 
star. 

Some Results for Merging Neutron Stars 

From the long list of our models we present here some results only for one exemplary case in 
more detail: the coalescence of a 1.8Mq primary with a 1.2Mq secondary (baryonic masses), 
both initially corotating (solid-body like). The results of this model were rather typical and 
not the most extreme in any direction. 

The neutrino luminosities during the neutron star merger are shown in Fig. |l|. During the 
phase of dynamical merging, a constantly rising neutrino luminosity is produced. It saturates 
after about 7 ms, when the secondary has effectively been tidally shredded and wrapped up 
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Figure 2: Contour plots of the mass density distribution in the orbital plane for the merging 
of a 1.2Mq and a 1.8Mq neutron star. The outer density contour corresponds to lO"*^*^ g/cm^, 
the inner to 10^^ g/cm^. In the top right corner of each panel one finds the time elapsed since 
the beginning of the simulation. 

along the surface of the primary. This evolution of the neutrino luminosity is very similar to 
what has been published for equal-mass neutron stars (Ruffert et al. Q), although somewhat 
higher luminosities of 2-5x10^^ erg/s were found in the recent simulations, because the use of 
the larger grid (out to about 200 km distance from the center of mass) allowed to keep track of 
the mass flung to very large radii and then falling back towards the center and contributing to 
the neutrino emission. Since the geometry and neutrino luminosities are similar for equal-mass 
and nonequal-mass neutron star mergers, also the energy deposition rate by i/z^-annihilation 
is of the same order. We calculated values of several 10^^ erg/s, about a factor of 10 larger 
than the ones given in |^, mainly because of the larger neutrino luminosities of the more 
recent models. However, the dynamical phase of the merging lasts only a few milliseconds 
before one must expect the collapse to a black hole of the merged object with significantly 
more than 2Mq. Therefore the energy pumped into a e^7 fireball during this phase hardly 
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exceeds a few 10 erg. Even worse, by far most of this energy is deposited in the surface-near 
regions of the merged stars and will drive a mass loss ("neutrino wind") which will pollute 
the fireball with an unacceptably large baryon load. 

The neutrino emission from the coalescence of two neutron stars is very different from the 
case of head-on collisions in which neutrinos are emitted in two very short (about 1 millisec- 
ond) and extremely luminous bursts reaching peak values of up to 4 x 10^'^erg/s Q! This 
gigantic neutrino luminosity leads to an energy deposition of about lO^'^ erg by i/j^-annihilation 
within only a few milliseconds. However, the surroundings of the collision site of the two neu- 
tron stars are filled with more than 1O~^M0 of ejected matter and the maximum values of 
the Lorentz-factor F are only about 10~^, five orders of magnitude lower than required for 
relativistic fireball expansion that could produce a GRB. 

Snapshots of the density contours in the orbital plane for the considered 1.2Mq-1.8Mq 
merger can be seen in Fig. ||. Initially (0ms < t < 3ms) the secondary is tidally elongated by 
the primary and a mass transfer is initiated. The top left panel of Fig. ^ nicely shows how 
matter is concentrated to flow through the Ll-point onto the surface of the primary. As more 
and more matter is taken away from the secondary it becomes ever more elongated (bottom 
left panel). Most of the matter of the secondary finally ends up forming a rapidly rotating 
surface layer of the primary, while a smaller part concentrates in an additional, extended thick 
disk around the primary. This matter has enough angular momentum to stay in an accretion 
torus even after the massive central body has most likely collapsed to a black hole. 

Neutron Tori around Black Holes and GRBs 

If the central, massive body did not collapse into a black hole, it would continue to radiate 
neutrinos with high luminosities, like a massive, hot proto-neutron star in a supernova. But 
instead of producing a relativistically expanding pair-plasma fireball, these neutrinos will 
deposit their energy in the low-density matter of the surface and thus will cause a mass flow 
known as neutrino-driven wind (e.g., [^). Most of this energy is consumed lifting baryons in 
the strong gravitational potential and the expansion is nonrelativistic. 

This unfavorable situation is avoided if the central object collapses into a black hole on 
a timescale of several milliseconds after the merging of the neutron stars. We simulated the 
subsequent evolution by replacing the matter inside a certain radius by a vacuum sphere 
(black hole) of the same mass. The region along the system axis was found to be evacuated 
on the free fall timescale of a few milliseconds as the black hole sucks up the baryons. Thus an 
essentially baryon-free funnel is produced where further baryon contamination is prevented by 
centrifugal forces. This provides good conditions for the creation of a clean e^7 fireball by z/P- 
annihilation. Also the thick disk closer to the equatorial plane loses matter into the black hole 
until only gas with a specific angular momentum larger than j* = 3i?s^^k(3-Rs) = V^GM/c 
{Rs = 2GM/cP is the Schwarzschild radius of the black hole with mass M, t;k(3i2s) the Kepler 
velocity at the innermost stable circular orbit at 3Rs) is left on orbits around the black hole. 
We find torus masses of up to Mt « 0.2-0.3 Mq at a time when a quasi-stationary state is 
reached. The temperatures in the tori are 3-10 MeV, maximum densities a few lO^^g/cm^. 
Typical neutrino luminosities during this phase are of the order of « lO^^erg/s (60% v'e, 
35% z/e). With a maximum radiation efficiency of « 0.057 for relativistic disk accretion 
onto a nonrotating black hole we calculate an accretion rate of Mt = L,^/ (c^Sr) ~ 1 -^o/s and 
an accretion timescale of tacc = Mt/Mt ~ 0.2-0.3 s. This is in very good agreement with the 
analytical estimates of Ruffert et al. . 
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From our torus models we find that the energy deposition rate by z^P-annihilation near the 
evacuated system axis is about Eyp ~ 5 x 10^'' erg/s. The corresponding annihilation efficiency 
ffa is of the order of Eyu/Ly ~ 5 x 10~^ which is rather large because of the large and Ve 
luminosities and because the torus geometry allows for a high reaction probability of neutrinos 
with antineutrinos. In order to account for an observed GRB luminosity of ~ 10^^ erg/s the 
required focussing of the z/P-annihilation energy is moderate, 50/(47r) = E^p/(2L^) ~ 1/4, 
corresponding to an opening angle of about 60 degrees. This value does not seem implausible 
for the geometry of a thick accretion torus. Taking into account general relativistic effects 
into the treatment of zvi/-annihilation reduces the above estimates only by about 10-50%. 

In summary, our numerical simulations show that stellar mass black holes with accretion 
tori that form after the merging of neutron star binaries have masses, lifetimes, and neutrino 
luminosities that might provide enough energy by i/i^-annihilation to account for short and 
not too powerful GRBs. The longer bursts and very energetic events, however, would require 
an alternative explanation, e.g., failed supernovae (or "collapsars" ) Q where a stellar mass 
black hole could have a more than 10 times more massive accretion torus than in the binary 
neutron star scenario. 
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Physical Processes Near Black Holes 

R.A. Sunyaev 
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The talk was devoted to the new information about black holes in the Galaxy and in the 
nuclei of the galaxies: 

• optical observations and the mass function of the invisible components of transient 
X-Ray binaries, 

• differences in X-Ray spectra of accreting neutron stars and black hole candidates 
(GRAN AT and MIR-KVANT data), 

• superluminal radiosources in two galactic X-Ray transients, 

• quasiperiodic oscillations of X- Rays from GRS 1915+105 (RXTE data), 

• scaling in the radiojet lengths in galactic and extragalactic superluminal radiosources, 

• X-Ray fluorescent lines in SGR B2 molecular cloud (ASCA data) and information about 
the X-Ray luminosity of the Center of our Galaxy in the recent past, 

• scattering of X-Ray lines on neutral and molecular gas (the spectrum of recoiled photons 
and the Lyman gap). 



High-Energy 
Neutrinos 



109 



Astrophysical Sources of High-Energy Neutrinos 

K. Mannheim 
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This contribution reviews currently discussed astrophysical sources of high-energy (>TeV) 
neutrinos. Particle astrophysics with high-energy neutrinos allows to study electroweak in- 
teractions, neutrino properties, particle acceleration theories, cosmology, dark matter candi- 
dates, and the enigmatic sources of cosmic rays |jl|. The detection of high-energy neutrinos 
is greatly facilitated by the rising neutrino cross section and muon range which cause the 
detection probablity for neutrinos to increase with energy. From an astrophysical point of 
view, the existence of sources of high-energy neutrinos seems very likely, if not guaranteed, 
for two strong reasons. 

First of all, there are numerous cosmic synchrotron sources (radio, optical. X-ray) for 
which simple arguments give electron energies in the GeV-TeV range [^. Since most elec- 
tromagnetic particle acceleration mechanisms predict accelerated baryons along with the ac- 
celerated electrons, hadronic energy losses due to pion production will lead to high-energy 
neutrino emission. If the maximum energies of particles emerging from cosmic accelerators 
are determined by the balance between energy gains and energy losses, baryons must reach 
much higher energies than electrons since their energy losses are much weaker than those of 
the electrons at the same energy. Therefore, pion energies may be expected to reach ener- 
gies even above TeV. Neutrinos from pion decay will, of course, also be very energetic, since 
each flavor carries ~ l/4th of the decaying pion energy. Baryon acceleration could only be 
avoided in a plasma exclusively composed of electrons and positrons. Although pair plas- 
mas probably exist in pulsars, gamma-ray bursts, and active galactic nuclei, baryon pollution 
(from the neutron star surface or due to entrainment of environmental plasma) changes the 
picture quite radically: As the pair plasma expands most of the internal energy is quickly 
converted into kinetic energy of the polluting baryons which then has to be tapped in a 
second stage dissipative process to give rise to the observed radiation. It has been widely 
discussed in the context of pair fireball models for gamma-ray bursts |Q that the observed 
non-thermal gamma-ray spectra generally require emission from a baryonic bulk flow at a 
large distance to the energizing compact object rather than from a pair plasma close to the 
compact object itself. 

Secondly, we know baryonic cosmic rays exist (up to energies of 10^ TeV) and they must 
come from somewhere. Supernova remnants are still among the prime candidates for the ac- 
celeration of the cosmic rays below ~ 10^ TeV, they should therefore be visible as gamma-ray 
and neutrino sources at comparable flux levels. This flux level is beyond reach for the first 
generation of high-energy neutrino experiments. Giant molecular clouds could be somewhat 
more important, since they provide most of the target mass in the vicinity of star forming 
regions where supernovae occur. Cosmic rays propagating through the Milky Way produce 
pions in hadronic interactions with the interstellar material and can be traced in gamma-rays. 
A diffuse flux of neutrinos at a flux comparable to the gamma-ray flux should be seen from the 
Galactic plane corresponding to ~ 150 upward events above one TeV per year in a detector 
with an effective volume of 1 km^ (compared with 6.5 x 10^ atmospheric background events). 
One can also speculate about enshrouded sources which only show up in neutrinos but not 
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in gamma-rays. However, the acceleration of particles to high energies generally requires a 
low-density plasma. The acceleration of cosmic rays above ~ 10^^'^ TeV in sources belonging 
to the Milky Way seems impossible. Among the prime candidates for these ultra-high energy 
cosmic rays are (radio and gamma-ray emitting) active galactic nuclei ||5| and gamma ray 
bursts p. In those sources cooling of accelerated protons (ions dissociate over cosmological 
distances) is dominated by interactions with low-energy synchrotron photons (from the accel- 
erated electrons). Due to electromagnetic cascading most gamma-rays are shifted to below 
the TeV range whereas the neutrinos remain unaffected from such reprocessing. Hence the 
7/i/-ratio decreases with energy in spite of being a constant function (from decay kinematics) 
at the production site (if proton-matter collisions dominate the energy loss, the j/u-iatio in 
the TeV range may be constant and of order unity [0]. The total electromagnetic power of the 
gamma-ray emitting active galactic nuclei can explain the observed extragalactic gamma-ray 
background and is of the same order as the power in an extragalactic E~'^ differential cosmic 
ray spectrum dominating above 10^'^ TeV. Therefore, neutrino flux predictions are generally 
bounded by ~ 10~^ GeV cm~^ s~^ st~^ (Fig. [l|). The models [^, |9| yield ~ 300 upward events 
above TeV per year per km^. 
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Figure 1: Overview of theoretical neutrino flux predictions as compiled by R.J. Protheroe 
during a workshop on high-energy neutrino astrophysics in Aspen, 1996. Labels refer to (a) 
atmospheric neutrinos (the solid part indicates the part of the spectrum measured by the 
Frejus experiment), (b) galactic disk in center and anti-center directions (due to Domokos), 
(c) AGN Q, (d) AGN |], (e) AGN [|, (f) charm production upper limit [0], (g) GRBs 
(due to Lee), (h) UHE CRs (due to Stecker), (i) TDs |o| (courtesy of W. Rhode). 
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A third, more speculative class of high-energy neutrino sources are topological defects 
producing gauge bosons at the grand-unified-theory scale of ~ 10^^ TeV |jT^. The unstable 
gauge bosons generate fragmentation jets with ultra-high energy protons, gamma-rays, and 
neutrinos which may be responsible for the extended air showers with primary energies above 
the so-called Greisen cutoff at 10^ TeV above which protons from sources further way than 
~ 30 Mpc lose their energy in interactions with the microwave background. 
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Atmospheric Muons and Neutrinos Above 1 TeV 

p. Gondolo 
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The ultimate background to high energy neutrino astronomy are the neutrinos produced in 
the interaction of cosmic rays with the Earth atmosphere. A cosmic ray nucleus penetrates on 
average ~ 60 g/cm^ of atmosphere before colliding with an air nucleus. The secondary mesons 
produced in the collision can propagate an additional ~ 100 g/cm^, but if their lifetime is 
short enough they may decay before interacting with air. Their semileptonic decays then 
produce an atmospheric neutrino flux. 

The competition between interaction with air and decay is at the basis of the energy 
and angular dependence of the atmospheric neutrino fluxes. Let us compare the mean free 
path of mesons in air to the distance they travel in a mean lifetime. In a cascade developing 
vertically downwards, the primary interaction occurs at a height of ~ 20 km, and the mean 
free path is ~ 10 km; in a cascade developing horizontally, the primary interaction is ~ 500 km 
from the observer, and the mean free path is « 100 km. The distance travelled in a mean 
lifetime increases linearly with the meson momentum: d = pr/m, where p, r, and m are the 
meson momentum, lifetime and mass. For p = 1 TeV, pions and kaons can travel kilometers 
= 55.9 km, dj^o = 31.2 km, dx+ = 7.52 km), but charmed mesons only centimeters 
{d£)0 = 6.65 cm, du+ = 1.70 cm, d^+ = 7.1 cm, = 2.7 cm). 

As a consequence: (1) vertical pions and kaons do not have time to decay before interact- 
ing with air, while horizontal pions and kaons decay readily: the subsequent "conventiona" 
neutrino flux is stronger in the horizontal than in the vertical direction; (2) charmed mesons 
decay promptly whether vertical or horizontal: the subsequent "prompt" neutrino flux is the 
same in the horizontal and in the vertical direction; (3) as the energy increases, the mesons 
live longer, and fewer and fewer pions and kaons have time to decay within the available mean 
free path, while all charmed mesons have time to decayj^ the energy spectrum of conventional 
neutrinos is steeper than that of prompt neutrinos by one power of energy. 

The different energy dependence means that beyond a certain energy, the atmospheric 
neutrino flux is dominated by decays of charmed mesons. The energy at which this hap- 
pens depends on the relative strength of charmed and non-charmed meson production in 
the primary nucleus-nucleus collision. Unfortunately, laboratory data do not extend up to 
the relevant energies. One could either look for help from measurements of the atmospheric 
muons produced in association with the neutrinos or extrapolate laboratory data with the 
help of a theoretical model for charm production. 

In a search for muons, prompt and conventional components could be separated by their 
different angular dependence. Several attempts have been made with underground exper- 
iments and with air shower detectors. In underground experiments, high energy muons, 
prompt muons in particular, come only from directions close to the horizon, where unfortu- 
nately there is an intense muon background generated by low-energy atmospheric neutrinos 
in the rock surrounding the detector. Subtraction of this neutrino-induced background is 
a source of big systematic errors [|l[] which spoil the search for prompt muons. Air shower 

^At high enough energy (~ EeV), also charm decays are fewer and fewer for an analogous reason. 
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detectors search for horizontal showers generated by a bremsstrahlung photon emitted by a 
high energy muon. The AKENO array has seen none above 100 TeV, and has therefore put 
only an upper limit to the prompt muon flux at high energy. 

Because of these experimental difficulties, we must resort to theory. And here lie many 
uncertainties, because charm production occurs half-way between the perturbative and non- 
perturbative regimes. The theoretically preferred model, perturbative QCD, was thought to 
be inadequate because to first order in Og it could not account for several aspects of the early 
data on open charm production, data which by the way were in conflict with each other |2|. 
So, even if a couple of QCD calculations had appeared ||3|, studies of atmospheric fluxes 
have traditionally favored, besides semi-empirical parametrizations of the cross section (see 
e.g. [Q, m, ^]), non-perturbative sources of charm beyond basic QCD: for example, the dual 
parton model j^, ^, based on Regge asymptotics, and models incorporating the assumption 
of an intrinsic charm component in the nucleon 



B 



^ 1^ 
"7 > 



10^ 



10' 



blazars 
py,pp — Tt— V 



atmospheric 
pA u,K— V 




ZHVa 
atmospheric 

pA ^ AAp^ V 
VFGS 



BNSZ 

GGHVZ 

ZHVe 

TIG 



TeV 



PeV 



EeV 



Figure 1: Selected predictions for high energy neutrino fluxes. 
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Today, however, charm production experiments give a consistent set of data, and are com- 



patible with QCD calculations to next-to-leading order |11|. This has removed the impetus 



to look for unconventional mechanisms of charm production, and has motivated a detailed 



study of atmospheric neutrinos within perturbative QCD [12|. 

In Fig. H I have collected some of the expected astrophysical neutrino signals and several 
predictions for the high energy atmospheric background. The astrophysical signals shown 



are the diffuse fluxes from blazars estimated by Dar and Shaviv ||l^, Protheroe |14], and 
Mannheim [^] , and the neutrino emission from the galactic interstellar medium as calculated 
by Domokos et al. fl^ . The conventional atmospheric background is from Lipari |l^, and 
the prompt neutrino fluxes are from Zas et al. (ZHVa for model A, ZHVe for model E) |4|, 
Volkova et al. (VFGS) [|, Bugaev et al. (BNSZ) fiol, Gonzalez-Garcia et al. (GGHVZ) Q, 



and Thunman et al. (TIG) [|12|. Also indicated are the upper limits on the prompt neutrino 



flux from MACRO [18| and AKENO [19|. Although these limits exclude the highest prediction 
ZHVa, the remaining uncertainty in the prompt neutrino background is still a factor of 100. 

Since the high-energy atmospheric neutrino background may be a nuisance to the detection 
of astrophysical diffuse neutrino fluxes, it may be of value to decrease the uncertainty in the 
predictions. The best way would be a detection of the associated prompt muons, maybe 
through larger air shower arrays. A second option would be to increase the precision of the 
theoretical inputs, for example through better measurement of the charm production cross 
section and the charm semileptonic decay rates, which is under way/under project in high 
statistics charm production experiments. At last, it may happen that we have to wait for the 
neutrino telescopes themselves to know the high energy atmospheric neutrino background. 
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Abstract 

The measurements of atmospheric neutrino interactions during 326 days of the Super-Kamio- 
kande detector operation are reported. The measured ratio of muon and electron neutrinos 
has been found smaller than expected from theoretical models. The angular distribution of 
muon neutrinos is also inconsistent with expectations. 

Introduction 

Atmospheric neutrinos are produced in a layer of about 15 km as a result of hadronic cascades 
originating from interactions of cosmic rays. A number of models have been developed to 
calculate the fluxes of atmospheric neutrinos [||, ^ at energies below 5 GeV. The absolute 
neutrino fluxes are predicted with uncertainties of 20%. However the ratio of i'^ + to 
+ fluxes is known to better than 5% and therefore experiments measure the double ratio 
R = (/-i/e)DATA/(^/e)MC5 where fi/e denote the ratio of /i-like to e-like neutrino interactions 
correspondingly in the data and simulated event samples. 

The flavor composition of atmospheric neutrinos was studied in the earlier underground 
experiments. The Kamiokande Q and 1MB Q experiments, using water Cherenkov detectors, 
and the Soudan Q experiment, using iron calorimeters, have reported the double ratio to be 
smaller than one. A dependence of R on the zenith angle, and hence on the neutrino path 
length, was also observed in Kamiokande ||6|. Those results have been often interpreted by 
neutrino oscillations. 

Here the preliminary results of measurements of atmospheric neutrino fluxes are presented 
for 20.1 kton-year exposure of the Super-Kamiokande detector. 

Super-Kamiokande Detector 

Super-Kamiokande is a 50 kton water detector located in the mine near Kamioka town in 
Japan. It is situated at a mean overburden of 1000 meters below the peak of Mt. Ikeno. 
Cherenkov photons emitted by relativistic particles inside the inner cylindrical volume of 
water, 16.9 m in radius and 36.2 m high, are recorded by 11146 photomultiplier tubes (PMTs) 
of 50 cm diameter (40% photocathode coverage). Each inner PMT signal provides the first 
photon arrival-time in 1.2 fis range at 0.3 ns resolution and the collected charge at 0.2 pC 
resolution (equivalent to 0.1 p.e.). The outer layer of water, about 2.7 m thick, is instrumented 
with 1885 outward facing PMTs of 20 cm diameter with wavelength shifting plates. The two 
detector regions are optically separated. The outer detector tags events with incoming or 
exiting particles. 

Water transparency is measured using a dye laser and is about 100 m at wavelenth of 
420 nm. The neutrino interactions were collected inside the fiducial volume 2 m from the 
inner detector walls, comprising 22.5 ktons of water. The accuracy of the absolute energy 
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calibration is estimated to be ±2.4% based on the detailed analysis of cosmic-ray muons, 
muon-decay electrons and electrons from a linac. The estimated momentum resolutions for 
electrons and muons are 2.5%/-\/-E'(GeV) + 0.5% and +3%, respectively. 

Results 

The data reduction, single-ring selection and lepton identification have been done by 2 groups 
independently and the results were found to agree very well. Three simulated samples have 
been compared with the data. The simulations used 2 different flux calculations, 3 different 
neutrino interaction models and 2 codes for particle and light propagation in the detector. 
All the simulated samples were processed by the same procedures as the data. The differences 
were used in estimates of systematic errors. The results discussed below come from one data 
analysis and Monte Carlo simulations based on flux calculations of Ref . fll . 

Out of about 300 million triggers recorded during 320 days of detector operation between 
May 1996 and June 1997, about 5000 events were classified as fully-contained events, i.e. 
events with no significant energy deposit outside of the inner detector. Using the photon 
time-arrival information recorded by PMTs, 2708 events with energy Svis > 30 MeV had 
vertices reconstructed in the fiducial volume. The visible energy E^is is defined as the energy 
of an electron which would produce the same number of Cherenkov photons. 

From Monte Carlo simulation it is estimated that 83.2% of all charged-current neu- 
trino interactions in the fiducial volume are retained in the current sample. The remain- 
ing fraction of 16.8% consists mostly of events with exiting particles (9.3%) or energy lower 
than 30 MeV (5.8%). 

To study the neutrino flavor composition we select a sample of single-ring events. Quasi- 
elastic interactions, for which the information about neutrino flavor is provided by the charged 
lepton identification, consist 75 it 3% of the sample. Particle identification is based on the 
pattern of Cerenkov rings. Diffuse, showering patterns associated with electromagnetic cas- 
cades (e-like events), are separated from sharper, non-showering rings caused by stopping 
muons (/x-like events). The misidentification probability for single-ring events is estimated 
to be 0.8 ±0.1%. 

The results for events with E^is < 1-33 GeV (Sub-GeV sample) are presented in Table |l|. 
For single-track events the electron and muon momenta are greater than 100 and 200 MeV/c, 
respectively. From these data one obtains: R = 0.63 ± 0.03(stat) ± 0.05(sys). The systematic 
error of 8% comes mostly form uncertainty in the calculation of the f^/fe ratio (5%), uncer- 
tainties in interaction cross sections and nuclear effects in ^^O (4.1%), single-ring selection 
(4%) and particle identification (2%). 

An independent signature of muon neutrino interactions is provided by electrons from 
muon decays. The detection efficiency for muon decay is estimated to be 0.80 for and 
0.63 for /i^. The fraction of /x-like events with a decay electron is 67.5 ± 1.7% for the data 
and 68.2 ± 1.0% for MC. For e-like events the corresponding fractions are 9.3 ± 1.1% and 
8.0 ± 0.3%. 

Angular distributions of the reconstructed track directions for single-ring events have also 
been studied. Figure |l| displays the distributions of zenith angles for /u-like and e-like events. 
The double ratio R dependence on the zenith angle is shown in Figure ^. 

Preliminary results of the analysis of the Multi-GeV sample of events with Ey[s > 1.33 
GeV is shown in Table The resulting double ratio is ii = 0.60 ± 0.06(stat) ± 0.07(sys). 
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Figure 1: Zenith angle distributions for: (a) e-like events, (b) ^u-likc, where cos0 = 1 for 
down- going particles. Histograms with hatched error bars show the MC prediction with its 
statistical error. Dotted histograms show the 25% systematic uncertainty on the absolute 
normalization, which is correlated between ^u-like and e-like events. 
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Figure 2: Zenith angle distribution of R. 
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Data 




Monte Carlo 






total 


l/e CC 


CC 


NC 


single ring 


1453 


1551.5 


539.7 


927.8.0 


84.0 


e-like 


718 


593.8 


537.2 


8.9 


47.7 


/x-like 


735 


957.7 


2.5 


918.9 


36.3 


multi ring 


533 


596.6 


148.4 


241.3 


206.9 


total 


1986 


2148.1 


605.6 


1076.5 


290.9 



Table 1: Summary of the Sub-GeV sample compared with Monte Carlo estimation, for 20.1 
kton-years of Super-Kamiokande. 





Data 


Monte Carlo 


single ring 


288 


304.8 


e-like 


149 


120.5 


//-like 


139 


184.3 


2 rings 


207 


221.0 


3 rings 


110 


133.0 


total 


605 


658.8 



Table 2: Summary of the Multi-GeV sample compared with Monte Carlo estimation, for 20.1 
kton-years of Super-Kamiokande. 

Finally the preliminary analysis of partially contained events, i.e. interactions with par- 
ticles exiting the inner detector, provided 156 data and 210.6 MC events. According to the 
simulations charged current and z^^ interactions constitute 97% of this sample. 

Conclusions 

The deficit of muon neutrinos reported from earlier experiments have been confirmed during 
the first year of data collection in Super-Kamiokande. The double rato values oi R = 0.63 it 
0.03(stat) ± 0.05(sys) for the Sub-GeV sample and R = 0.60 ± 0.06(stat) ± 0.07(sys) for the 
Multi-GeV sample are consistent with results from Kamiokande, 1MB and Soudan. 

The angular asymmetry observed with much smaller statistical weight in Ref. is ob- 
served in the present sample for both Sub-GeV and Multi-GeV cuts. Those results strongly 
suggest the neutrino oscillation interpretation. However, more data are needed to exclude 
any possibility of detector related biases. 
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Introduction 

The AMANDA experiment plans the construction of a neutrino detector of about 1 km^ size 
(ICECUBE). Deep (> 2 km) holes are melted into the thick ice-cover of Antarctica at the 
South Pole. Large (8 inch) photomultiplier tubes, embedded in glass pressure housings, are 
attached to an electro-optical cable (vertical separation of 10-20 m) and lowered into the hole 
before re-freezing starts. By deploying several of these "strings" a 3 dimensional matrix of 
optical sensors is achieved. 
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Figure 1: The AMANDA- A and AMANDA-B detectors. 
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Key signatures for nucleon interactions below the detector are upward going muon 
tracks traversing the detector, i/g, Ut can be identified by large electromagnetic or hadronic 
cascades not associated with a down-going muon track. Muon tracks emit Cerenkov light with 
a cone of fixed angle. By measuring the arrival time and number of photons the direction and 
energy of the muon is reconstructed. The muon direction is pointing close into the direction of 
the initial neutrino, thus enabling us to operate the detector as a neutrino telescope. The main 
background comes from down-going atmospheric muons, which may be falsely reconstructed 
as up- going tracks. In case of showers the vertex of the interaction can be reconstructed from 
the times of the spherically propagating Cerenkov photons. 

The rate of neutrinos produced in the Earth's atmosphere is falling off more rapidly 
with higher energies than expected for cosmic sources. Therefore the AMANDA detector is 
optimised for muon detection above energies of 1 TeV. However, muons are detected down to 
energies of a few GeV with decreasing efficiency. 

A sketch of the currently operating AMANDA-A (4 shallow strings) and AMANDA-B 
(10 deep strings) detectors is shown in figure ||. In a next stage, AMANDA-II, 11 additional 
1 km long strings are being installed around the AMANDA-B detector. The first three strings 
of 1.2 km length are being deployed]^ during the winter season 97/98. The results of these 
strings will outline the cubic km scale technology 

Besides the main purpose — the search for point-sources of high energy neutrinos (e.g. from 
active galactic nuclei, AGN) — a large variety of research topics are covered, like the measure- 
ment of the total neutrino fluxes (from all AGN), neutrinos in coincidence with gamma ray 
bursts, atmospheric neutrinos, or neutrinos from decays of exotic dark matter. Due to the 
low temperatures and clean environment the noise rates of the photomultipliers are small. 
A special DAQ system continuously monitors these noise rates in different time windows, 
allowing to search for bursts of low energy (MeV) neutrinos. Already in the present detector 
a supernova within our galaxy would yield a statistically significant excess in the summed 
count-rate of all sensors. Another possible source of low energy neutrino bursts could be due 
to gamma-ray bursts g, 1,3. 

The South Pole Site 

Besides an excellent local infrastructure and good transportation logistics provided by the 
American Amundsen-Scott South-Pole station the construction and long-term operation of a 
large neutrino telescope in the deep Antarctic ice is suggested by a variety of advantages Q. 

A continuous 3 month access (and year round maintenance) is possible. Long good weather 
periods, 24 h daylight and large available space allow complex deployment operations and 
preparations on the stable ice surface. This includes the use of heavy equipment and parallel 
work e.g. at two holes. 

Relatively short cables are sufficient to connect the 2 km deep strings to the surface data 
acquisition system. While only simple and robust components are buried deep into the ice, 
more vulnerable electronic components are located at the surface, reducing the risks of failure. 
The short distance also allows to connect each PMT by its own cable to the surface electronics. 
The analog PMT signals are transmitted to the surface via electrical and optical fibres. In 
case of failure of an optical connector the electrical one still has sufficient accuracy to serve as 
a backup. The PMT anode current drives an LED directly without using complex electronics. 



^By the deadline of this report all 3 have been successfully installed, covering a depth from 1200m to 2400m. 
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Thus an Amanda optical module consist of a bare PMT with an LED and a voltage divider 
in a pressure housing. It is connected to a high voltage cable and an optical cable via two 
connectors. 

Though increasing the price per channel, the non-hierarchical approach of individual cables 
provides great advantages, when operating in a non-laboratory environment. First, the risk 
of single-point failure, e.g. the failure of essential components is minimised — single PMTs 
may fail but not larger parts of the detector as in the case of several PMTs connected to one 
cable or "digitisation module." Secondly, an evolutionary approach during the construction 
of AMANDA-B allowed to improve certain components year by year, based on practical 
experience. Older installations are still operated, without the need to redesign the detector. 
Despite of the fact that details of a later km-size detector are not decided yet, it is obvious 
that AMANDA-B can serve as a nucleus for the ICECUBE P detector. 

The very stable cold temperatures in the deep ice result in low noise-rates of PMTs and 
a high reliability. No photomultiplier, which was successfully deployed in 1994 or later, was 
lost since then. No daily or annual changes of the environment are observed and the PMT 
sensitivity is not degraded by effects such as sedimentation. 

The South-Pole site is largely used by various astrophysical facilities. One interesting 
aspect are coincidences between AMANDA and the surface air-shower detectors SPASE and 
GASP. Shower cores, which point from SPASE to AMANDA, give an important source of 
muons of known direction. These events are widely used for calibration and verification of 
reconstruction algorithms. Several hundred coincident events per year originate from cosmic 
rays of energies above the knee of the primary cosmic ray spectrum. SPASE aims to use the 
information e.g. on the muon number, measured with AMANDA, to determine the cosmic 
ray composition at high energies. 

Results from Amanda-B 

As a major result from the analysis of laser light propagating through ice it was found, 
that the optical absorption is extremely large, > 200 m at 800 m depths and about 100 m 
at 1800 m, considerably exceeding the values for clearest natural water. Residual bubbles 
at shallow depths of AMANDA-A lead to a small scattering length. Scattering of photons 
improves the calorimetric properties for energy reconstruction but lead to a loss of time- 
information and thus makes the reconstruction of track directions difficult. At the larger 
depths of AMANDA-B the scattering length is significantly improved (see figure ^). The 
challenge of track reconstruction for the remaining level of scattering was solved. The energy 
resolution seems to be superior to water experiments, e.g. BAIKAL. 

The data which was recorded in 1997 has been transferred from the Pole and is currently 
being analysed. However, using the data taken in 1996 with only 4 installed strings, first 
results have been obtained. 

As an example for neutrino search in AMANDA-B (4 strings), figure ^ shows the re- 
constructed zenith angle distribution of recorded muon events. From top to bottom quality 
criteria for filtering of up-going events are subsequently tightened. Without quality criteria 
(top) a large fraction of events are falsely reconstructed to originate from the lower hemi- 
sphere. These events can be rejected in good agreement with the MC expectation. The right 
picture shows a neutrino candidate event filtered out of the 4-string data. The event has a 
vertical length of about 300 m and is reconstructed as up-going. 

The performance for the completed AMANDA-B detector (10 strings) is expected to be 
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Photon arrival times (ns) Photon arrival times (ns) 

Figure 2: Ice properties measured at 1000 and 2000 m depth. The figures show the distribu- 
tions of photon arrival times relative to a laser source, measured with PMTs at 20 m (left) 
and 40 m (right) distance to the source. Scattering due to residual air bubbles leads to a sig- 
nificant increase of the photon propagation time (and a broadening of the spectrum) . This is 
a strong effect for AMANDA-A depths. It is much smaller at depths of 2 km (AMANDA-B). 
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Entries: 2637 
Mean: 5.3 
Median: 2.4 
MPV: 1.7 
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Figure 4: Performance expected for the 10 string AMANDA-B detector. The left figure 
shows the zenith angle distribution of up-going muon events passing all quality cuts. In this 
calculation minimum ionizing muons have been initially generated isotropically from the lower 
hemisphere. The difference in reconstructed and generated muon directions are shown on the 
right plot. A median accuracy of 2.4° is achieved. 

significantly better compared to the 4-string installation. The data is currently being analysed. 
Figure |^ shows results from a full Monte Carlo calculation. With background suppression 
(S/N > 10^) an effective area of the order of 10 • 10^ m^ and a median angular resolution of 
2.4° is achieved. Several atmospheric neutrinos per day should be detectable 

The full AMANDA-II detector is aimed to be installed until the year 2000. It should reach 
effective areas from 50 • 10^ m^ to 100 • 10^ m^. It is thus almost two orders of magnitude 
more sensitive than current detectors. The angular resolution is expected to be of the order 
of 1° g. 

Towards ICECUBE 

First experiences with data of a km scale telescope have already been gathered in 1996 with 
coincidences between the shallow AMANDA- A and the deep AMANDA-B detector [||, ||, ||] . 
The purpose of the current AMANDA-II installation campaign is to monitor the depth profile 
of optical properties from 1200 m down to 2500 m. The tests of a large variety of technical 
issues will help to define a robust technology, which is suitable for the construction of the large 
km^ scale telescope. Besides the tests of different PMTs, optical fibres and connectors also 
wavelength-shifters (to improve the UV sensitivity) and other strategies for signal processing 
(Digitising optical modules) are tested 

A straw-man ICECUBE detector might consists out of 60 additional 1 km long strings, 
each with 70 OMs. Based on the experience from AMANDA-B and AMANDA-II the prices 
per channel may be estimated to 6 k$. This includes the optical module, cables, calibrations 
sources, surface electronics (DAQ) and workshop labour. Adding 5 M$ for off-line computing 
and additional detector components the bare detector price would be of the order of 30 M$ Q. 
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However, in order to operate the detector and perform the instahation within a reasonable 
time (5 years) additional logistic costs have to be taken into account. These are mainly 
the costs of modernisation the drilling equipment, a new building, fuel, cargo and technical 
personal. These support costs can be estimated to about 10 M$ @. 

Adding the costs of AMANDA-B and AMANDA-II, the ICECUBE detector would cost 
of the order of 50 M$ §. 

It is planned to start the installation in the beginning of year 2001 and finish it by 2005. 
A workshop on the ICECUBE project will be held in March 1998 p. 
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Neutrino astrophysics has already generated considerable interest in both the particle physics 
and astrophysics communities, through the detection of solar neutrinos, atmospheric neu- 
trinos and neutrinos from supernova 1987A. ANTARES aims to extend this scientific 
reach to higher energies (10 GeV-10 PeV) where potential neutrino sources include Active 
Galactic Nuclei, Gamma Ray Bursters, pulsars. X-ray binaries, young supernovae remnants 
and neutralino annihilation in the center of the Sun or Earth. These neutrino observations 
should provide complementary information that might be difficult or even impossible to ob- 
tain through gamma or cosmic ray observations. The observable energy range for neutrinos 
extends well above 1-10 TeV where extra-galactic gamma rays are severely attenuated over 
cosmic distances by their interactions with infra red and cosmic microwave background pho- 
tons. Moreover, unlike charged cosmic rays, neutrinos are not deflected by magnetic flelds 
and thus point back to their source. 

There are several potential astrophysical sites for accelerating cosmic rays up to the ener- 
gies of lO^'^ eV that have been observed. In general, these sites will also produce significant 
fluxes of neutrinos and gamma rays through the decays of charged and neutral pions, which 
are produced in approximately equal numbers through soft interactions of the high energy 
protons with other protons or photons present in the acceleration site. The emitted fluxes of 
neutrinos and gamma rays from this mechanism would thus be roughly equal, however the 
gamma rays may be severely attenuated at the source by interactions with matter or photons, 
leaving the neutrinos as the best probe of the physics. 

Several authors (see and references therein) have constructed models which account for 
the observed TeV gamma ray emission of Active Galactic Nuclei (AGNs) in terms of proton 
acceleration in the jets which emerge from the central black hole of 10^-10^^ solar masses. 
These authors predict between 100 and 1000 events per year from all AGNs in a neutrino 
detector with 1 km^ effective area. Another possible cosmic acceleration site is the expanding 
fireball produced by coalescing neutron stars immediately after merging. These events are 
sufficiently energetic and numerous to explain the cosmic origin of the observed Gamma Ray 
Bursts (GRBs). Waxman and Bahcall Q predict between 30-100 neutrino events per km^ per 
year. Pulsed neutrino sources such as GRBs are easier to detect since the background from 
atmospheric neutrinos is greatly reduced by looking for time and direction coincidences with 
gamma-ray observations of GRBs with satelite experiments. They also offer the possibility of 
kinematic neutrino mass tests through the relative times of the neutrino and gamma signals at 
the Earth. Potential galactic acceleration sites of cosmic rays include pulsars. X-ray binaries 
and young supernova remnants. Calculations show that for these Galactic sources, total 
proton luminosities of a few percent of the Eddington limit could be detected by a km-^ 
neutrino detector. 

A very interesting signal for ANTARES arises if the non-baryonic dark matter is primarily 
in the form of weakly interacting massive particles, the favoured particle physics candidate 
being the lightest supersymmetric particle (LSP), normally assumed to be the neutralino 
[x)- A relic neutralino undergoing an elastic scatter while passing through a star or planet 
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could be gravitationally captured, the resulting accumulation of neutralinos continuing until 
balanced by XX annihilation. Final states from xx annihilation include many particles which 
subsequently decay to final states including energetic neutrinos which can be detected by 
ANTARES. If we assume that the dark matter in the halo of our galaxy consists of neutralinos, 
then a large fraction of the currently allowed SUSY parameter space predicts neutrino fluxes 
from XX annihilation in the centres of the Sun and of the Earth which would be detectable Q 
with exposure times of around 10^ m^ yr. 

High energy muon-neutrinos can be detected by observing long-range muons produced by 
charged current neutrino- nucleon interactions in the matter surrounding the detector. An 
array of photomultiplier tubes, placed in an optically transparent medium such as deep ocean 
water or antarctic ice, can be used to detect the Cerenkov photons produced along the muon 
tracks. Since Cerenkov radiation is emitted in a cone of fixed angle (42 degrees in water) the 
direction of the muon, which preserves the neutrino direction, can be reconstructed from the 
location and time of the PMT 'hits'. Simulation studies of prototype ANTARES detectors 
show that the angular resolution for reconstructed muons could be better than 0.2°. Moreover, 
for muons with energies above 1-10 TeV, the average angle between the muon and the parent 
neutrino is lower than 0.1°. Above 1 TeV, the average number of Cerenkov photons produced 
per meter of muon track is proportional to the muon's energy so that the number of detected 
photons can be used to measure the muon energy, albeit with poor resolution due to the 
stochastic nature of muon energy loss above 1 TeV. 

Muons from cosmic ray interactions in the Earth's atmosphere produce a large background 
of down-going muons which is many orders of magnitude greater than any potential neutrino 
signal. Thus the aperture of the telescope is restricted to the 27r solid angle of up-going muon 
events where this background is negligible, provided that the detector missreconstructs less 
than 1 in 10^ (for a detector at 2.5 km depth) down-going events as up-going events. If 
this is achieved then the dominant background source of up-going muons will be from up- 
going atmospheric neutrinos which can interact, just like the cosmic neutrinos, in the rock 
or water below the detector. This background is isotropic and falls steeply with energy: 
angular resolution and energy resolution will therefore play a critical role in determining 
the detector's ability to find point sources of neutrinos above this background. Electron 
and tau flavor neutrinos can also be detected via the Cerenkov photons produced from the 
electromagnetic and hadronic showers from neutrino interactions inside or near the detector. 
The pointing accuracy and event rates for these events will be significantly less favourable 
than for muon events but the energy resolution should be better. 

ANTARES is currently in an R&D phase where the technology for constructing a large- 
area neutrino detector in the mediterranean sea is being developed through deployments at 
a test site which is 30 km from the French coast near Toulon and at a depth of 2300 m. 
In 1998 a first prototype string will be deployed and data from 8 PMTs transmitted to 
shore by a 40 km long electro-optical cable which is ready to be laid in place by France 
Telecom. By the end of 1999 one or two fully-equipped strings and up to 100 PMTs will be 
deployed in an array that will be directly scalable to the large-area detector. Great emphasis 
is being placed on developing a reliable and cost-effective deployment strategy by using well 
established commercial technology from companies such as France Telecom and by involving 
two Oceanographies institutes (IFREMER and Centre d'Oceanologie de Marseille) directly 
in the collaboration. These two institutes have many years experience of deploying scientific 
instruments at depths relevant to ANTARES. 

In parallel with these activities, ANTARES has already started deploying autonomous 
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test strings at the Toulon site and at a site near Corsica. These strings wih be used to 
select the optimal site for the large-area detector by measuring the relevant environmental 
parameters: background noise on PMTs due to ^"K decays and bioluminescence, the rate at 
which glass surfaces lose transparency from sedimentation and biological growth and, most 
importantly, the transparency of the water in terms of absorption and scattering lengths as 
a function of wavelength. Twelve deployments of these strings have been carried out so far, 
showing encouraging results that can be accessed from the ANTARES web site 

After 1999, when the R&D phase is completed and the site for the large array chosen, 
ANTARES aims to deploy an array with an effective area for up-going muons of order 0.1 
km^ over a timescale of 3 years. Further development of the array up to 1.0 km^ area is 
anticipated but would depend on the results obtained to that point. 
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Introduction 

The first attempts to detect very high energy (VHE) gamma rays from cosmic point sources 
with ground based detectors go back to the early 1960s. After a long period in which there 
was only evidence from one point source of VHE gamma-rays (the Crab nebula) , in the last 
few years technological advances in background suppression have led to the discovery of six 
additional sources, three of them extragalactic. The study of these sources gives clues to the 
still very incompletely understood acceleration mechanisms and sites of cosmic particles at 
nonthermal energies. Moreover, the raised sensitivities for the first time allow to critically 
test predictions for theories of cosmic-ray origin. A very intense outburst of TeV gamma-rays 
observed in 1997 from the active galaxy Mkn 501 is an example for objects that emit a major 
part of their total energy output in the high energy range. VHE gamma-ray astronomy is thus 
becoming a part of main-stream astrophysics, but there are also possibilities for investigations 
which are unique to this energy range, particularly in cosmology. 

Status of Instrumentation 

In the energy range from 200 GeV to about 20 TeV the Cerenkov telescope technique to 
detect gamma-rays dominates presently. The primary gamma ray develops an electromagnetic 
cascade which emits Cerenkov light. This light is detected in an optical telescope with typical 
mirror areas from about 5 to 70 m^. The Cerenkov light is emitted over an area of about 20000 
m^ and each telescope situated whithin this pool can detect the shower. The resulting very 
large detection area leads to very high count rates. For example, the Whipple collaboration 
detected in one observing season (Spring 96) about 5000 photons from the active galaxy 
Mkn 421 above 300 GeV This is about the same number as a satellite experiment, 
EGRET at the Compton gamma-ray observatory, detects in one year above 100 MeV, where 
the flux is typically a few thousand times higher, from all active galaxies (about 50 active 
galaxies were detected above 100 MeV)! 

Until recently the sensitivities were seriously degraded by the fact that the hadronic cosmic 
radiation produces a background which even for the most intense sources is about a factor 
100 higher than the signal from photons. The imaging technique, pioneered by the Whipple 
collaboration allows to reduce this background. The image of the developing shower is 
registered in an array of fast photomulipliers which detect the few nanosecond long Cerenkov 
pulses. The showers produced by hadrons are much more irregular in shape than the ones 
from photons, and this allows to reject the background showers with an efficiency reaching 
e.g. 98 % for the HEGRA telescope 1, a single telescope of modest size with a camera of 
moderate angular resolution (127 pixels with a field of view of 3°) another telescope of 
the same experiment is shown in Fig. |^). 
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Figure 1: Typical Cerenkov telescope (HEGRA telescope 2 at La Palma). The main mirror 
consists of 30 submirrors with a total area of 8.5 m^. The energy threshold of this telescope 
for 7-rays is about 1.1 TeV. In the background quadratic boxes from the HEGRA array which 
operates at energies above 14 TeV are visible. 

This technique is still in a state of evolution, the most advanced camera presently in 
operation by the French CAT collaboration Q has 546 pixels. Its high angular resolution 
(0.12°) allows not only a further reduction in background but also a low energy threshold 
and a deduction of the core position of the shower which in turns allow a very good energy 
resolution. Another crucial technological advance for Cerenkov telescopes is the exploitation 
of stereo imaging. The same shower is viewed by several telescopes at the same time and its 
full three dimensional structure can be reconstructed. Besides advantages similar to a higher 
angular resolution the measurement of shower parameters is overdetermined and thus allows 
important consistency checks of the technique |^]. 

Array techniques where the lateral distribution of Cerenkov light and/or particles at 
ground level is sampled with high spatial resolution are an important technological alter- 
native to imaging telescopes. At higher energies (above 10 TeV in existing detectors) they 
offer the distinct advantage of a field of view on the order of a steradian, comparable to 
satellite detectors. The array technique will also be exploited by using heliostat mirrors in 
existing solar-power plants to reach low energy thresholds . 
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Figure 2: Sensitivities in units of photons cm~^ sec~^ of various existing and planned detectors 
for high and very-high energy gamma rays for various energies (fuU hues) or at the energy 
threshold (big dot). The numbers are for a total measuring time of 50 hours for devices with a 
small field of view and 1 year for devices with a field of view of 1 steradian or larger [satellites: 
EGRET (existing) and GLAST (planned) and the array MILAGRO (planned)]. WHIPPLE 
and HEGRA are existing telescopes, VERITAS and MAGIC proposed large future devices. 
STACEE, CELESTE and GRAAL are detectors which are currently being set up using the 
large mirror areas at solar power plants. The dotted line gives the extrapolation of a source 
spectrum from GeV to VHE energies (active galaxy Mkn 421). The arrows indicate the 
intensity of this source during various viewing periods ( "VP" ) of the EGRET satellite. 



Recent results 

This summer (1997) the Australian- Japanese CANGAROO collaboration announced detec- 
tion of radiation above 2 TeV from the remnant of the supernova (SNR) of A.D. 1006 0. The 
detected flux was close to previous theoretical predictions on the basis of the hypothesis 
that cosmic-ray electrons are accelerated by the shock front of SNRs. The accelerated elec- 
trons suffer inverse Compton scattering on the cosmological 3 K background radiation, which 
are then detected in the TeV range. This is strong experimental support for the idea that the 
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Galactic cosmic-ray electrons (which account for about 1 % of the total cosmic-ray flux at 
earth) are due to shock wave acceleration in SNRs. Perhaps this brings us near to a solution 
of the historic question of the origin of the main hadronic part of cosmic rays: it had long been 
speculated that all cosmic rays with energies below about 100 TeV are accelerated in SNRs. 

The VHE upper limits on the SNR G78. 2+2.1 shown as an example in Fig. |3| |^ present 
a problem for this scenario, though. This SNR lies near a dense molecular cloud which 
should act as a "target" for the accelerated cosmic-rays, leading to the prediction of a copious 
production of VHE gamma-rays via the reaction p + p — > ttq ^ 2 7. It is seen that the 
upper limits from various experiments shown in Fig. ^ are below theoretical predictions in 
the mentioned scenario This has recently led to a revival of speculations on an 

extragalactical origin of the main part of hadronic cosmic rays |11]. 

The active galaxy Mkn 501 at a distance of about 100 Mpc from earth showed a spectacular 
outburst of TeV radiation 1997 which was observed in detail by 5 different collaborations |12|. 
At the highest activity level in April 1997 the flux was about 40 times higher than during the 
discovery of the source in 1995. Figure § shows the light curve determined by the HEGRA 
telescope 1. This light curve is the most complete of all measured ones because data were also 
taken during moonshine. The outburst was also observed in the X-ray and optical region and 
these multiwavelength studies are especially valuable to discriminate between various models 
for the origin of the TeV radiation and to find the reasons for its fluctuation. The measured 
VHE spectra already allow interesting conclusions about the intensity of cosmological infrared 



background radiation field [13|. Interactions of VHE photons with cosmological background 



fields will allow a variety of cosmological studies unique to this wavelength range 1 14 1 . 




Figure 3: Upper hmits on VHE radiation from the SNR G78.2+2.1. Whipple and HEGRA 
are Gerenkov telescopes, AIROBICC, CASA-MIA and CYGNUS are arrays. It now seems 
likely that the detection of 7-rays from a satellite above 100 MeV (EGRET) is due to a pulsar. 
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Figure 4: Gamma-ray flux above 1.5 TeV from the direction of the active galaxy Mkn 501, 
determined with the HEGRA telescope 1. The various symbols are for different modes of 
operation with and without ambient moon light. The time axis is given in mean Julian days 
and extends from March 8 to October 13, 1997 (courtesy D. Kranich). 



Projects for the Future 

VHE astrophysics seems to be in a similar state as X-ray astronomy in the 1960s. The next 
generation of telescopes which will combine higher angular resolution, larger mirror area and 
stereoscopic imaging . It will probably allow to detect on the order of hundreds of sources 
and lower the energy threshold thus closing the present "gap" in detectable primary gamma 
energy between satellites (up to 20 GeV presently) and ground based devices. 
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Helium Absorption and Cosmic Reionization 

Craig J. Hogan 

Departments of Physics and Astronomy, University of Washington 
Box 351580, Seattle, WA 98195, USA 

In my talk at Ringberg I emphasized the current concordance of hght element abundances 
with the Big Bang predictions and with the observed density of baryons. The substance of 
most of these remarks together with further references can be found in recent papers |^. 
Here I want to highhght one aspect of this puzzle which has not been commented on very 
much, the interesting current state of observations on singly ionized helium at high redshift. 

The gas in the emptier regions of the universe at the highest redshift observed (almost 5) 
is aheady almost entirely ionized; the fraction of hydrogen which is neutral is less than 10~^. 
It is small enough that over most of redshift space, the absorption from hydrogen has small 
optical depth. Hell's higher ionization potential (54.4 eV as opposed to 13.6 eV for hydrogen) 
means that its ionization is delayed relative to hydrogen as such hard photons are relatively 
rare. Hell is therefore more abundant in absolute terms than any other species which makes 
it the best tool for studying matter in the emptier regions of space — the voids between the 
density concentrations. 

Currently there are two quasars where published high resolution HST/GHRS quasar ab- 
sorption line spectra permit the separation of the more diffuse component of the Hell res- 
onance absorption from the component in the concentrations of matter which appear as HI 
Lyman-a forest clouds ||^, ^. Both datasets imply an upper bound on the diffuse matter 
density at z ~ 3 of less than O.Ol/i^^'^, derived based on the maximum permitted mean Hell 
Lyman-a optical depth allowed after subtracting the minimal contribution from the detected 
HI Lyman-a forest clouds, while adopting the hardest ionizing spectrum allowed by the data. 
This upper bound is important for constraining ideas about the evolution of the baryons and 
especially the possibility of a large repository of baryons in the voids. Both datasets also agree 
that there is Hell absorption at nearly all redshifts — that is, there is at least some matter 
everywhere, even in the emptiest voids. 

However the current published results do conflict in one important respect. Hogan et 
al. IQ] find in Q0302-003 that although the Hell optical depth is rather high (greater than 1.3 
everywhere), there is also detectable quasar flux at all redshifts in between the identified HI 
clouds, whereas Reimers et al. Q find in another quasar (HE2347-4342) large regions where 
the flux is consistent with zero. The first result suggests that there is rather little Hell and 
that the helium is mostly doubly ionized Helll already (unless the voids are swept implausibly 
clean of gas), whereas the second result suggests that large volumes of space are still mostly 
Hell and that helium ionization is not yet complete. 

Since the ionization is expected to occur in patchy domains around the strongest sources 
of ionizing photons (indeed near 0302-003 itself the sphere of influence is observed in the 
spectrum to about 4000 km/sec in radius), it is possible that both interpretations are correct 
and that the different directions are just different. This would be interesting because the 
scale of the inhomogeneity in ionization history would be much larger than expected, with 
likely consequences for large scale structure in the galaxy distribution. On the other hand 
it could also be that one or both of the published spectra have an incorrect zero level — a 
natural suspicion since they are taken with the one-dimensional GHRS for which background 
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subtraction requires considerable modeling. If this is the case our data are more suspect than 
those of Reimers et al. since Q0302-003 is a fainter quasar. 

The reason I highlight this situation now is that it is likely to be resolved soon as the 
new STIS detector is now being trained on both of these targets and should with its better 
background subtraction resolve this issue. By the time this note is published there should be 
a new result which will give us a direct insight into the final stages of cosmic reionization. It 
should also allow a much tighter measurement of the range of densities allowed; at present the 
lower limit on the density is about two orders of magnitude below the upper limit, because of 
the uncertainty in modeling the ionization. In addition, new higher resolution observations 
will allow a better defined separation of diffuse and cloud components, and observations of 
a third quasar (PKS 1935-692) should help to differentiate between the accidental properties 
of particular sightlines and universal properties of the intergalactic medium. 
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Non-Standard Big Bang Nucleosynthesis Scenarios 

K. Jedamzik 

Max-Planck-Institut fiir Astrophysik, Karl-Schwarzschild-Str.l, D-85748 Garching 
Abstract 

An overview of non-standard big bang nucleosynthesis (BBN) scenarios is presented. Trends 
and results of the hght-element nucleosynthesis in BBN scenarios with smah-scale or large- 
scale inhomogeneity, the presence of antimatter domains, stable or unstable massive neutrinos, 
neutrino oscillations, neutrino degeneracy, or massive decaying particles are summarized. 

Introduction 

Light-element nucleosynthesis during the BBN epoch below cosmic temperatures T ~ 1 MeV 
contributes significantly to "^He, ^He and ''Li abundances and likely, to all of the abun- 
dance observed throughout the universe. BBN is a freeze-out process from nuclear statistical 
equilibrium such that light-element abundance yields are sensitively dependent on the cosmic 
conditions during the BBN era as well as the properties of neutrinos governing the freeze-out 
process from weak equilibrium. Calculations of abundance yields in a standard BBN scenario 
are performed under the assumptions of a universe homogeneous in the baryon-to-photon 
ratio, with massless neutrinos and vanishing neutrino chemical potentials, and in the absence 
of massive decaying particles or other degrees of freedom. For reviews on the physics of 
standard BBN see [||]. I presently summarize results and trends in theoretical calculations of 
non-standard BBN scenarios where one of the above assumptions is relaxed. This summary 
is in no way intended to be complete in discussing all possible modifications to a standard 
BBN scenario or in referencing the thousands of articles on non-standard BBN. I apologize for 
including only key references due to the limited scope of these proceedings and for presenting 
my personal view of the field of non-standard BBN. For two excellent reviews on non-standard 
BBN the reader is referred to The determination of observationally inferred primordial 
abundance constraints, which represents possibly the most important branch of the field of 
big bang nucleosynthesis at present, will not be discussed here. In what follows, abundance 
yields in non-standard BBN will either be given in relation to abundance yields in standard 
BBN or in absolute values and shall be understood as indicative of approximate trends. 

Non-Standard BBN 

In the following a list of non-standard BBN scenarios, their respective modifications to stan- 
dard BBN (hereafter; SBBN), as well as trends and results in these scenarios are given. 

a) Inhomogeneity 

The baryon-to-photon ratio r] is the one free parameter in SBBN. Any inhomogeneity in this 
quantity results in modified nucleosynthesis yields which depend on the typical amplitude 
and spatial seperation scale of inhomogeneities. Substantial changes in the abundance yields 
only result when 5r]/r]^ 1. 
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i) Inhomogeneity in the Baryon-to-Photon Ratio on Small Mass Scales: 

Fluctuations in rj which may arise on sub-horizon scales at earlier cosmic epochs as, for exam- 
ple, possibly during a first-order QCD transition or electroweak transition, result in a highly 
nonstandard BBN scenario. The nucleosynthesis in an environment with ry-fluctuations is 
characterized by coupled nuclear reactions and hydrodymamic processes, such as baryon dif- 
fusion and late-time expansion of high-density regions . Fluctuations in r] persist down to 
the onset of BBN provided the mass of an individual high-density region exceeds 10~'^^ Mq. 
One of the main features of such scenarios is the differential diffusion of neutrons and protons 
leading to the existence of neutron- and proton-rich environments. The trend in inhomoge- 
neous BBN is the overabundant production of ^He when compared to SBBN at the same 
average rj, nevertheless, there exists parameter space where less ^He than in SBBN is synthe- 
sized. Scenarios which don't overproduce ^He typically have high (^H/H) ~ 1—2 x 10~^ and 
high (''Li/H) ~ 10~^ — 10"*^. Such BBN may agree with observational abundance constraints 
for fractional contributions of the baryon density to the critical density Of, about 2-3 times 
larger than in SBBN, but only in the seemingly unlikely advent of efficient ^Li depletion in 
population II stars. 

ii) Inhomogeneity in the Baryon-to-Photon Ratio on Large Mass Scales: 

When the baryonic mass within a typical fiuctuation exceeds (M^ 10~^^Mq), baryon diffu- 
sion and hydrodymanic processes during the BBN era are of no significance such that BBN 
abundance yields may be given by an average over the SBBN abundance yields of seperate re- 
gions at different ij. For non-linear fiuctuations exceeding the post-recombination Jeans mass 
M> lO^M©, which may exist in primordial isocurvature baryon (PIB) models for structure 
formation, early collapse of high-density [r]) regions is anticipated [^. The nucleosynthesis 
yields of collapsing regions may be excluded from the primordial abundance determination if 
either dark objects form or significant early star formation in such high-density regions occurs. 
If only low-density regions contribute to the observable primordial abundances, characteristic 
average abundance yields for scenarios designed to possibly agree with observationally in- 
ferred primordial abundances are: (^H/H) 1—3 x 10~^, (^Li/H) 5 x 10^^°— 2 x 10^^, and 
^He mass fraction Yp ~ 0.22—0.25 at a total $75^0.2 (i.e. including possible dark objects), 
larger than inferred from a SBBN scenario [||, ^. One feature of such models is the predic- 
tion of fairly large intrinsic spatial variations in the primordial abundances, which may be 
observationally tested by (^H/H) determinations in Lyman-limit systems ||5|. These models 
may only agree with observationally inferred abundance limits when there are no fiuctuations 
below M ^ 1O^M0 and collapse efficiencies of high-density regions are large. 

Hi) Matter/ Antimatter Domains: 

A distribution of small-scale matter /antimatter domains in baryon- asymmetric universes 
(i.e. where net r/ 7^ 0) may result from electroweak baryogenesis scenarios. If the baryon 
(antibaryon) mass in individual domains is ^ 10~^^ Mq the BBN process in such scenarios 
is characterized by differential diffusion of neutrons (antineutrons) and protons (antiprotons) 
which causes a preferential annihilation of antimatter on neutrons 0. When annihilation 
of antimatter occurs before significant ^He synthesis (T ^ 80keV) but after weak freeze-out 
(T ^ 1 MeV) a modest to substantial reduction of Yp results. When annihilations occur mainly 
after ^He synthesis the dominant effect is significant production of ^He and |^, with •^He/^H 
ratios likely to be in confiict with observational constraints. 
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b) Non-standard Neutrino Properties 

The BBN process may be approximated as the incorporation of all available neutrons into 
^He nuclei at a temperature T ~ 80keV. The neutron abundance at this temperature, and 
hence the final ^He mass fraction Yp, may be increased with respect to a SBBN scenario due 
to an increased expansion rate of the universe during the BBN era. An increased expansion 
rate raises the neutron-to-proton ratio (hereafter; n/p) at weak freeze-out and reduces the 
time for neutron decay to decrease the n/p ratio between weak freeze-out and significant ^He 
synthesis. We will refer to this effect as the "expansion rate effect." In addition, the n/p 
ratio at T ~ 80 keV may be either decreased or increased by introducing additional electron- 
and/or anti-electron neutrinos into the plasma. In SBBN it is assumed that the left-handed 
neutrino and right-handed antineutrino seas of three massless, stable neutrino flavors i^e, 
and Vr are populated and that neutrino chemical potentials do vanish. Modifications of these 
assumptions usually result in either the expansion rate effect or additional (anti) electron 
neutrinos, or both. The principal effect of such modifications is to change the ^He abundance, 
and to a less observationally significant degree the abundances of other light-element isotopes. 

i) Massive, Long-Lived t -Neutrinos: 

Neutrinos are considered massless and long-lived in the context of BBN for neutrino masses 
rriu ^ 100 keV and lifetimes Ti, ^ 10^ s [see ii) below, however, for possible photodisintegration] . 
A massive, long-lived r-neutrino leads to the expansion rate effect since the contribution to 
the total energy density from the rest mass of r- neutrinos continually increases as the universe 
expands between weak freeze-out and ^He synthesis, possibly even resulting in matter domi- 
nation during the BBN era. BBN with massive, long-lived r neutrinos and for experimentally 
allowed z^T-masses therefore results in increased ^He and useful limits on the allowed mass of 
a long-lived r-neutrino have been derived Q. 

ii) Massive, Unstable r-Neutrinos: 

The effects of decaying r neutrinos on the light-element nucleosythesis [^] sensitively depend 
on the decay products. One distinguishes between (i) decay into sterile particles, in particular, 
particles which interact neither weakly with nucleons interchanging neutrons and protons nor 
electromagnetically with the ambient plasma (e.g. Vj- ^ v^+cj), where cj) is a, weakly interacting 
scalar), (ii) decay into sterile and electromagnetically interacting particles (e.g. Vr ^ t^^i + j), 
and (iii) decay into (anti) electron neutrinos and sterile particles (e.g. v-r ^ i^e + ) 0' llll- 
For decay channel (i) and t^-^Is increased ^He mass fraction results due to the expansion 
rate effect which, nevertheless, is weaker than for long-lived r-neutrinos since the energy 
of the (massless) decay products redshifts with the expansion of the universe. In contrast, 
for life times t^^'Z^s and m^^ ^ 10 MeV it is possible to reduce the Yp since effectively the 
distributions of only two neutrino flavors are populated [13]. Decay channel (ii) would have 



interesting effects on BBN but is excluded by observations of supernova 1987A for T,y^ ^ 10^ s. 
For decay via channel (iii) additional (anti) electron neutrinos are injected into the plasma 
and their effect depends strongly on the time of injection [0]. When injected early (tj^^ ~ 1 s), 
the net result is a reduction of Yp [|l^] since weak freeze-out occurs at lower temperatures. 
In contrast, when injected late (r^^ ~ 10^ — lO^s) the resulting non-thermal electron neutrinos 
affect a conversion of protons into neutrons, yielding higher Yp and/or higher depending 
on injection time. It had been suggested that a scenario with late-decaying I'j- via channel (iii) 



may result in the relaxation of BBN bounds on by a factor up to ten |11|, nevertheless, this 
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possibility seems to be ruled out now by the current upper laboratory limit on the z/^-mass. 
For life times t^^ ^ 10^ s, modifications of the light-element abundances after the BBN era by 
photodisintegration of nuclei may result (cf. radiative decays). 

in) Neutrino Oscillations: 

Neutrino oscillations may occur when at least one neutrino species has non-vanishing mass 
and the weak neutrino interaction eigenstates are not mass eigenstates of the Hamiltonian. 
One distinguishes between (i) flavor-changing neutrino oscillations (e.g. v^. i/^) and (ii) 
active-sterile neutrino oscillations (e.g. ^ t^s)- Here Vs may be either the right-handed 
component of a [v^, Vj-) Dirac neutrino, or a fourth family of neutrinos beyond the standard 
model of electroweak interactions. In the absence of sterile neutrinos and neutrino degeneracy, 
neutrino oscillations have negligible effect on BBN due to the almost equal number densities 
of neutrino flavors. When sterile neutrinos exist, neutrino oscillations may result into the 
population of the sterile neutrino distribution, increasing the energy density, and leading to 
the expansion rate effect. The increased Yp has been used to infer limits in the plane of 
the neutrino squared-mass difference and mixing angle In the presence of large initial 
lepton number asymmetries (see neutrino degeneracy) and with sterile neutrinos, it may be 
possible to reduce Yp somewhat independently from the detailed initial conditions, through 



the dynamic generation of electron (as well as /i and r) neutrino chemical potentials [17| 



c) Neutrino Degeneracy 

It is possible that the universe has net lepton number. Positive net cosmic lepton number 
manifests itself at low temperatures through an excess of neutrinos over antineutrinos. If 
net lepton number in either of the three families in the standard model is about ten orders 
of magnitude larger than the net cosmic baryon number BBN abundance yields are notably 
affected. Asymmetries between the (i^r) and (v'r) number densities result in the expan- 
sion rate effect only, whereas asymmetries between the and number densities induce a 
change in the weak freeze-out (n/p) ratio as well. Since the expansion rate effect leads to 
increased ^He production (i^^) degeneracy may be constrained. However, one may find 
combinations of f^, Vr, o,nd chemical potentials which are consistent with observational 



abundance constraints for fit much larger than that inferred from SBBN |18]. Nevertheless, 
such solutions not only require large chemical potentials but also an asymmetry between the 
individual chemical potentials of fg and {vt)- Asymmetries between the different flavor 
degeneracies may be erased in the presence of neutrino oscillations. 



d) Massive Decaying Particles 

The out-of-equilibrium decay or annihilation of long-lived particles (r^^O.ls), such as light 
super symmetric particles, as well as non-thermal particle production by, for example, evapo- 
rating primordial black holes or collapsing cosmic string loops, during, or after, the BBN era 
may significantly alter the BBN nucleosynthesis yields |]l9|. The decays may be classified ac- 
cording to if they are (i) radiative or (ii) hadronic, in particular, whether the electromagnetic 
or strong interactions of the injected non-thermal particles are most relevant. 

i) Radiative decays: 

Electromagnetically interacting particles (i.e. 7, e^) thermalize quickly in the ambient photon- 
pair plasma at high temperatures, such that they usually have little effect on BBN other than 
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some heating of the plasma. In contrast, if radiative decays occur at lower temperatures after a 
conventional BBN epoch, they result in a rapidly developing cascade which only subsides 
when individual 7-rays of the cascade do not have enough energy to further pair-produce 
on the ambient cosmic background radiation. The net result of this cascade is a less rapidly 
developing 7-ray background whose properties only depend on the total amount of energy 
released in electromagnetically interacting particles and the epoch of decay. At temperatures 
T ^ 5 keV the most energetic 7-rays in this background may photodissociate deuterium, at 
temperatures T^0.5keV the photodisintegration of ^He becomes possible. Radiative decays 
at T « 5 keV may result in a BBN scenario with low and low ^He if an ordinary SBBN 
scenario at low rj is followed by an epoch of deuterium photodisintegration |2^. Radiative 
decays at lower temperatures may produce signficant amounts of and ^He. Nevertheless, 
this process is unlikely to be the sole producer of ■^H due to resulting observationally disfavored 
^He/^H ratio. The possible overproduction of ^He and by the photodisintegration of ^He 
has been used to place meaningful limits on the amount of non-thermal, electromagnetically 
interacting energy released into the cosmic background. These limits may actually be more 
stringent than comparable limits from the distortion of the cosmic microwave background 
radiation Q. 

ii) Hadronic decays: 

The injection of hadrons into the plasma, such as vr^'s, vr'^'s and nucleons, may affect light- 
element nucleosythesis during, or after, the BBN era and by the destruction and production 
of nuclei. In general, possible scenarios and reactions are numerous. If charged hadrons are 
produced with high energies below cosmic temperature T ^ keV they may cause a cascade 
leading to the possibility of photodisintegration of nuclei (see radiative decays). Through 
charge exchange reactions the release of about ten pions per nucleon at T ~ 1 MeV results in 
a significant perturbation of weak freeze-out and increased Yp [|2^. A fraction of only ~ 10~^ 
antinucleons per nucleon may cause overproduction of ^He and through antinucleon-^He 
annihilations. It is also conceivable that the decaying particle carries baryon number such 
that the cosmic baryon number is created during the BBN era. A well-studied BBN scenario is 
the injection of high-energy (~ 1 GeV) nucleons created in hadronic jets which are produced 
by the decay of the parent particle ||2^. If released after a conventional BBN era these 
high-energy nucleons may spall pre-existing ^He, thereby producing high-energy ^H, ^H, and 
^He as well as neutrons. Such energetic light nuclei may initiate an epoch of non-thermal 
nucleosythesis. It was found that the nucleosynthesis yields from a BBN scenario with hadro- 
destruction/production and photodisintegration may result in abundance yields independent 
of r] for a range in total energy injection, and half-life of the decaying particles. Nevertheless, 
such scenarios seem to produce ^Li in conflict with observational constraints. 

e) Other Modifications to BBN 

There are many other variants to a standard big bang nucleosythesis scenario which have 
not been mentioned here. These include anisotropic expansion, variations of fundamental 
constants, theories other than general relativity, magnetic fields during BBN, superconducting 
cosmic strings during BBN, among others. The influence of many, but not all, of those 
scenarions on BBN is due to the expansion rate effect. Studies of such variants is of importance 
mainly due to the constraints they allow one to derive on the evolution of the early universe. 
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Big Bang Nucleosynthesis With Small-Scale 
Matter- Antimatter Domains 

J.B. Rehm and K. Jedamzik 

Max-Planck-Institut fiir Astrophysik, Karl-Schwarzschild-Str.l, D-85748 Garching 

Big Bang Nucleosynthesis (BBN) is one of the furthest back-reaching cosmological probes 
available. By means of comparing the predicted and observationahy inferred light element 
abundances the physical conditions as early as a few seconds after the Big Bang can be 
scrutinized [p. Nevertheless, presence of matter-antimatter domains during BBN on length 
scales comparable in size to the neutron diffusion length at temperatures between 10 and 
10~^ MeV (see e.g. Q) has so far not been investigated. To our knowledge, the influence 
of antimatter on BBN has only been studied |^] with respect to a homogeneous injection of 
antiprotons into the primordial plasma after the end of BBN (T < 10 keV). The conclusion 
of such studies is that only a small fraction of about 10~^ antiprotons per proton is allowed 
to be injected after the BBN epoch for the abundances of deuterium and ^He produced in 
p^He reactions not to exceed the observationally inferred values. 

In this work Q we discuss the BBN in a baryo-asymmetric universe with a inhomogeneous 
distribution of matter-antimatter domains. For possible scenarios which may yield matter- 
antimatter domains during an electroweak baryogenesis scenario see e.g. Refs. ^. As long as 
the diffusion length of neutrons and antineutrons is shorter than the size of a typical domain, 
matter and antimatter remain segregated. Note that proton diffusion is suppressed due to 
Coulomb scattering. Therefore, annihilations may likely only proceed when (anti)neutrons 
diffuse out of their domains. We identify three main scenarios: 1. Annihilations before weak 
freeze-out: No effect on the light-element abundances because the n/p-ratio, which governs the 
abundances, is reset to the equilibrium value by the fast weak interactions. 2. Annihilations 
between weak freeze-out (« 1 MeV) and the onset of ^He synthesis (T ~ 0.08 MeV): The 
n/p-ratio is strongly affected. Annihilation of antimatter mostly occurs on neutrons, since 
neutrons can both be annihilated by antineutrons diffusing into the matter region and diffuse 
themselves into the antimatter region to annihilate there on antiprotons and antineutrons. 
Protons on the other hand are confined to the matter region and can only be annihilated by 
diffusing antineutrons. By providing this effective sink for neutrons, the ^He abundance may 
be substantially lower than in a standard BBN scenario at the same net entropy-per-baryon, 
simply because there are less neutrons left to build up ^He. In extreme cases, when all neutrons 
have annihilated on antimatter before the onset of ^He synthesis, zero net ^He production 
results. This may well represent the only known BBN scenario in a baryo-asymmetric universe 
which precludes the production of any light elements. 3. Annihilations well after BBN: This 
scenario was already mentioned above. Annihilations of only small amounts of antimatter on 
^He may easily overproduce the cosmic deuterium and ^He abundances. 

In summary, the figure illustrates that matter-antimatter domains present at the BBN 
epoch may have dramatic effects on the net synthesized ^He abundance. By comparison, the 
abundances of other light isotopes are generally less affected. A more detailed analysis of 
BBN with matter-antimatter domains will be presented in a forthcoming work |^]. 
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Figure 1: Predicted ^He mass fraction as a function of the typical antimatter domain size. 
Length scales are given in proper units at a temperature of T = 100 MeV. The different lines 
give results for different amounts of antimatter, i.e. the ratio of R = A^antimatter /-^matter takes 
values of 0, 1, 10, and 20%, respectively. All models have the same net baryon density. 
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Neutrinos and Structure Formation in the Universe 

Matthias Bartelmann 

MPI fiir Astrophysik, P.O. Box 1523, D-85740 Garching, Germany 
Abstract 

I review the standard theory of structure formation in the Universe, starting from the assump- 
tions that (i) the Universe is globally described by a Priedmann-Lemaitre model, (ii) structure 
forms via gravitational instability from primordial density fluctuations, and (iii) these primor- 
dial density fluctuations were Gaussian in nature. I describe how density fluctuations grow 
in time, what the density perturbation power spectrum looks like, and why hot dark matter 
like neutrinos has a considerable impact on the process of structure formation. 

Introduction 

The simplest and most widely accepted theory of structure formation in the Universe starts 
from the following three assumptions: 

1. Globally, the Universe is well described by the Priedmann-Lemaitre-Robertson-Walker 
(FLRW) model. 

2. Structure formed via gravitational instability from primordial density fluctuations. 

3. These primordial density fluctuations were Gaussian in nature. 

By the first assumption, the global dynamics of the Universe is described in terms of four 
parameters, viz. the density of ordinary and relativistic matter, the Hubble constant, and the 
cosmological constant. The second assumption, together with the dynamics implied by the 
first, specifies how structures grow in time. The third assumption asserts that the statistics of 
the density fluctuation field is completely specified by its mean and variance. By definition, 
the mean density contrast is zero, and the variance is specified by the power spectrum, 

P{k) = {\5\k)\) , (1) 

where k is the wave vector of the density perturbation. The power spectrum itself does not 
depend on the direction of k because of the (assumed) isotropy of the density fluctuation field. 

These three basic assumptions provide the structure of this article. I will first describe the 
global dynamics of FLRW models, then discuss how density fluctuations grow in an expanding 
universe, and finally describe the gross features of the density fluctuation power spectrum, 
before I come to the conclusions. 

Global Dynamics 

The scale factor a of a FLRW universe is governed by Friedmann's equation. 
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where I have ignored the cosmological constant for simphcity. Neglecting the curvature term 
(as one can do at early times, see below), and approximating a ~ a/t, eq. (§) implies the 
familiar result that the expansion time scale t oc The scale factor a grows monotonically 

with the cosmic time t and can therefore be used as a time variable. Boundary conditions 
for (^) are chosen such that a(to) = 1 at the present cosmic time to) and t = for a = 0. 
The present value of the Hubble parameter, H = a/a, is called the Hubble constant Hq. K is 
the curvature of spatial hypersurfaces of space-time, and p is the total matter density. The 
quantity 

p,,^|^«2xl0-29/.2gcm-3 (3) 

is called the critical density, for reasons that will become clear below. The density parameter 
r^o is the current total matter density /o(to) = Po in units of 

f^o-^. (4) 

Per 

The Hubble constant is commonly written as 

Hq = 100/ikms~^Mpc~^ « 3.2 x 10"^*^ /is^^ , (5) 

where 0.5 < /i < 1 expresses our ignorance of Hq. Note that the Hubble parameter has the 
dimension (time)~^, so that H^^ provides the natural time scale for the expansion of the 
Universe. Evaluating Friedmann's equation (^) at time it follows 



K 



In other words, a universe which contains the critical matter density is spatially flat. 
The first law of thermodynamics, dU -\- pdV = 0, can be written as 

For ordinary matter, pcP S> p ~ 0, hence p oc a~^. For relativistic matter, p = pc^/S, hence 
p oc a~^. Therefore, the matter density changes with a as 

p{a) = po a"" = Per fio a"" , (8) 

where n = 3 for ordinary matter ( "dust" ) , and n = 4 for relativistic matter ( "radiation" ) . 
Summarising, we can cast Friedmann's equation into the form, 



H\a) = Hi 



17o a"" - (^^0 - 1) o~ • (9) 



Since the universe expands, a < 1 for t < to, and so the expansion rate H{a) was larger in 
the past. At very early times, a ^ 1, the first term in eq. (P) dominates because n > 3, and 
we can write 

H{a) = Ho 0^/2 a-"/2 _ (10) 



This is called the Einstein-de Sitter limit of Friedmann's equation. 
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Photon and Neutrino Backgrounds 

The Universe is fiUed with a sea of radiation. It has the most perfect black body spectrum 
ever measured, with a temperature of T^^ = 2.73 K. Apart from local perturbations and 
kinematic effects, this Cosmic Microwave Background (CMB) is isotropic to about one part 
in 10^. It is a rehc of the photons that were produced in thermal equihbrium in the early 
Universe. The energy density in the CMB is 

and the equivalent matter density is p^^ = u^fl/c^. We have seen before that the matter 
density in radiation changes with the scale factor like a~^, hence 

r^(a) = a"^T^,o • (12) 
Converting the matter density in the CMB to the density parameter as in eq. (^, we find 

o = ^ = 2.4 X 10-5 /i-^ (13) 

The energy density in radiation today is therefore about five orders of magnitude less than 
that in ordinary matter. 

Like photons, neutrinos were produced in thermal equihbrium in the early universe. Be- 
cause of their weak interaction, they decoupled from the plasma at very early times when 
the temperature of the universe was kT ~ 1 MeV. Electrons and positrons remained in equi- 
librium until the temperature dropped below their rest-mass energy, kT ~ 0.5 MeV. The 
consequent annihilation of electron-positron pairs heated the photons, but not the neutrinos 
since they had already decoupled. Therefore, the entropy of the electron-positron pairs Se 
was completely dumped into the entropy of the photon sea S.y. Hence, 

('S'o ~l~ 'S'-j- ) before ~ 

(5^) 

after ) 

(14) 

where 'before' and 'after' refer to the electron-positron annihilation time. Ignoring constant 
factors, the entropy per particle species is S" oc gT^, where g is the statistical weight of the 
species. For bosons, g = i, and for fermions, g = 7/8 per spin state. Before the annihilation, 
we thus have the total statistical weight (/before = 4-7/8+2 = 11/2, while after the annihilation, 
Rafter = 2 because only photons remain. From eq. (|l4[). 



-^after \ ffbefore H /i r\ 

(15) 



before / 9after ^ 

After the electron-positron annihilation, the neutrino temperature is therefore lower than the 
photon temperature by the factor (4/11)^/'^. In particular, the neutrino temperature today is 

/ 4 \ i/3 

7;,o=fYT) ^7,0 = 1-95 K. (16) 



Although the neutrinos are long out of thermal equilibrium, their distribution function re- 
mained unchanged since they decoupled, except that their temperature gradually dropped. 
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We can therefore calculate their energy density as for the photons, and convert this to the 
equivalent cosmic density parameter. The result is 

J7,,o = ^ = 2.8x10-6/1-2 (17) 

per neutrino species. Likewise, their number density is given by the Fermi distribution with 
temperature T^Oi 



per neutrino species. Suppose now there is one neutrino species which has rest mass rrii, > 0, 
while the other species are massless. For the dark matter to be dominated by neutrinos, rrii, 
needs to satisfy 

rriu riufi = VLq pcr , (19) 

hence, 

m^^c^ ~ 95eV(J^o/i^) • (20) 

Adding the energy density of the two massless neutrino species to that of the photons, the 
total equivalent mass density in relativistic particles is 

1-5 1,-2 



17r_o = 3.0 X 10"^ h'^ . (21) 



The Horizon 



The size of causally connected regions of the Universe is called the horizon size. It is given 
by the distance a photon can travel in the time since the Big Bang. Since the appropriate 
time scale is provided by the inverse Hubble parameter, the horizon size is = cH{a)~^, 
and the comoving horizon size is 

dnia) = ^ = ^n-'/'a-/'-\ (22) 
an [a) Hq 

where we have inserted the Einstein-de Sitter limit ( [lO| ) of Friedmann's equation. cHq^ = 
3h~^ Gpc is called the Hubble radius. 

The previous calculations show that the matter density today is completely dominated 
by ordinary rather than relativistic matter. But since the relativistic matter density grows 
faster with decreasing scale factor a than the ordinary matter density, there had to be a time 
Oeq <C 1 before which relativistic matter dominated. The condition 

n^^oa-^ = noa-^ (23) 

yields 

Oeq = 3.0 X 10"^ i^oh^r^ . (24) 

We shall see later that the horizon size at the time a^q plays a very important role for 
structure formation. Under the simplifying assumption that matter dominated completely 
for all a > Ceq (i.e. ignoring the contribution from radiation to the matter density), eq. ( p2|) 
yields 

(^H(aeq) = ^ n-'^^ al{^ ^ 13 {^lo h^r' Mpc . (25) 
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The (photon) temperature at Oeq is 

fcr^(aeq) = a^^ kT^fl ~ 8 eV . (26) 

This is much lower than the rest-mass energy of a cosmologicaUy relevant neutrino. Hence, 
if there is a neutrino species with mass given by eq. (pO[), it becomes non-relativistic much 
earlier than Ogq. 

Growth of Density Fluctuations 
a) Linear Growth 

Imagine a FLRW background model into which a spherical, homogeneous region is embedded 
which has a slightly different density than the surroundings. Because of the spherical symme- 
try, this sphere will evolve like a universe of its own; in particular, it will obey Friedmann's 
equation. Let indices and 1 denote quantities within the surrounding universe and the 
perturbed region, respectively. Then, Friedmann's equation reads 

^ . Kc^ SttG 

Hi = ^Po (27) 

for the perturbation and for the surrounding universe, respectively. We have ignored the 
curvature term in the second equation because the background has K = at early times. 

We now require that the perturbed region and the surrounding universe expand at the 
same rate, hence 

Ho = Hi. (28) 
Assuming that the density inside the perturbed region only weakly deviates from that of the 



background, we can put ai ^ ao, and eqs. (|27D and (|2^ ) yield 



^ _ pi- PO 3Kc 1 2x-l /on\ 

6 = = — -— 2 /'oa ) . (29) 

po SttG poa^ 

We have seen before that po oc a~", with n = 4 before Oeq and n = 3 thereafter. It thus 
follows that 

6o.h" before ae„ 
1 a alter Ogq. 

This heuristic result is confirmed by accurate, relativistic and non-relativistic perturbation 
calculations. 

b) Suppression of Growth 

A perturbation of wavelength A is said to "enter the horizon" when A = dn- If A < (iH(ffleq)) 
the perturbation enters the horizon while radiation still dominates. Until Oeq, the expansion 
time scale, iexpj is determined by the radiation density /jr, and therefore it is shorter than 
the collapse time scale of the dark matter, Idm- 

texp ~ (Gpr)-1/2 < (Gpdm)"'/' ~ iDM • (31) 
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In other words, the fast radiation-driven expansion prevents dark-matter density perturba- 
tions from cohapsing. Light can only cross regions that are smaller than the horizon size. 
The suppression of growth due to radiation is therefore restricted to scales smaller than the 
horizon, and larger-scale perturbations remain unaffected. This suggests that the horizon size 
at Oeq, d}i{aeq), scts an important scale for structure growth. 
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Figure 1: Sketch illustrating the suppression of structure growth during the radiation- 
dominated phase. The perturbation grows oc before Ceq, and oc a thereafter. If the 
perturbation is smaller than the horizon at Ogq, it enters the horizon at Oenter < Oeq while 
radiation is still dominating. The rapid radiation-driven expansion prevents the perturba- 
tion from growing further. Hence it stalls until Ogq. By then, its amplitude is smaller by 
/sup = (oenter/oeq)^ than it would be without suppression. 

Figure || illustrates the growth of a perturbation with A < (iH(ffleq)) that is small enough 
to enter the horizon at agnter < Oeq- It can be read off from the figure that such perturbations 
are suppressed by the factor 

f^^^=(^] . (32) 

\ "cq / 

According to eq. (22), the comoving horizon size scales with a as dn oc a"/^~^. When a 
perturbation enters the horizon, A = dnloenter), which yields Oenter c>c A^/^"""^-* = k~'^/'^^~'^\ 
or Oenter OC A = during the radiation-dominated epoch. Thus we obtain from (^) 



where /cq is the wave number corresponding to di{{a, 



eqy 



The Perturbation Spectrum 

Consider now the primordial perturbation spectrum at very early times, Pi{k) = \5f{k)\. Since 
the density contrast grows as (5 oc a"~^, the spectrum grows as P{k) oc a^*^""^). At Oenter, the 
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spectrum has therefore changed to 

/Inter (A:) = 0,^"'^ n{k) OC P,{k) . (34) 

It is commonly assumed that the total power of the density fluctuations at Center should 
be scale-invariant. This implies A;^Penter(^) = const., or Penter(fc) . Accordingly, the 

primordial spectrum has to scale with k as -Pi (A;) oc k. This scale-invariant spectrum is 
called the Harrison-Zel'dovich-Peebles spectrum. Combining that with the suppression of 
small-scale modes, we arrive at 

forA:>A:o. ^^^^ 



Damping by Free Streaming 

If the particles of the dominant dark- matter component are fast enough ("hot"), perturbations 
can be washed out by free streaming. To see this, consider the so-called free- streaming length 
Aps, that is the length that particles can cover within the available time t. It is given by 

AFs(t)=a(t) r^dt'. (36) 
Jo a{t'j 

While particles are relativistic, v = c, and later v oc because the momentum is red-shifted. 
We thus introduce a new scale, the scale factor a^^ where the dark-matter particles become 
non-relativistic. It is defined by the condition kT = a^^kT^^ ~ myC? . Inserting T^^o = 2.73 K 



and the neutrino mass rriy from eq. (20), we find 

a^, = 4.7 X 10"^ {^oh^y^ . (37) 
Evaluating eq. ([3^), the comoving free-streaming length turns out to be 



2ct, 



anr 



for a < a 



Till 



1 + In y — j for anr < a < Oeq, ^3g-j 

5 . -, f fleq \ 3 / flgq \ ^ 



2+'" e -5 if ) 



where Aps is plotted in Fig. |2| as a function of a. This figure shows that Aps quickly reaches 
an asymptotic value A^. For dark matter dominated by neutrinos, this asymptotic free- 
streaming length is 

Ai?s = 11-3 Mpc {Vt^h^)-^ . (39) 

We therefore conclude that if neutrinos were to dominate the dark matter, fluctuations 
with A < Apg would be wiped out by free streaming. Such perturbations cannot keep the fast 
neutrinos bound. Consequently, the power spectrum is exponentially cut off at wave numbers 
k > feps = (-^Fs)"^- Since /cps is only slightly larger than the wave number where the power 
spectrum reaches its peak, ko, the neutrino-dominated power spectrum is cut off right above 
the peak. This is illustrated in Fig. |3| which contrasts cold with hot dark matter spectra. For 
cold dark matter (CDM), free streaming is unimportant, so the spectrum behaves as described 
by eq. (|35|). For hot dark matter (HDM), the spectrum is cut off exponentially at large 
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Figure 2: The comoving free-streaming length AFs(a) as a function of the scale factor a. Aps 
grows rapidly until am and quickly approaches an asymptotic value Ap^ thereafter. 




k [h/Mpc] k [h/Mpc] 

Figure 3: Left panel: Power spectrum P{k) for cold and hot dark matter (dashed and solid 
lines, respectively). For cold dark matter, the spectrum falls off oc k"^ at large k, while free 
streaming of hot dark matter like neutrinos exponentially damps away perturbations smaller 
than the asymptotic free-streaming length Xpg. Right panel: Total power k^P{k) for cold 
and hot dark matter. While the power continues to increase monotonically iox k ^ oo for 
cold dark matter (solid curve), it reaches a peak at A; = fco for hot dark matter and decreases 
rapidly on smaller scales. 

k. Consequently, the total power of the density fluctuations, k'^P{k), grows monotonically 
with k for CDM, while it peaks around ko for HDM. Structure begins to form on scales 
where the power is largest. For CDM, the smallest structures are formed first, and structure 
formation proceeds towards larger scales. For HDM, structures on scales of A ~ A:q"^ form first. 
They correspond to objects somewhat more massive than galaxy clusters. Smaller structures 
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Figure 4: Growth of structure in CDM (top) and HDM (bottom) from the same initial random 
density perturbations; from left to right: z = 0, z = 1, and z = 2. High-density regions are 
black. 

form later by fragmentation. Thus, free streaming has the important consequence that the 
bottom-up scenario of CDM is changed to the top-down scenario of HDM. The formation of 
structures smaller than galaxy clusters, e.g. galaxies, is therefore strongly delayed in HDM 
models. Figure ^ illustrates this. It shows two cosmological simulations starting from the 
same initial random density perturbations, one with CDM and one with HDM. Three epochs 
of the evolution at z = 2, 1, are shown. While significant structure on small scales is already 
seen at z = 2 in CDM, structure then just begins to grow on much larger scales in HDM. 

Summary 

The main results of the preceding discussion can be summarised as follows: 

• The matter density of the Universe today is dominated by ordinary rather than rela- 
tivistic matter. Going back in time, the radiation density grows faster than the matter 
density. The epoch when they equal each other plays an important role for structure 
formation. It is characterised by the scale factor 

= 3.0 X 10"^ (1^0^^)"^ • 

• Linear density perturbations grow oc before Oeq, and oc a thereafter. 
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• Rapid radiation-driven expansion earlier than Ogq suppresses the growth of modes which 
are smaher than the size of the horizon at Oeq, 

dH(aeq)~13 (Oo/t')-^ MpC . 

Density perturbations with wave number k > d^^{aeq) = ko are suppressed by the 
factor 



sup 



k 



Together with the assumption of a scale-free primordial power spectrum, this yields the 
cold dark matter spectrum 



P{k) oc 



k {k < ko), 
k~^ (fc>fco). 



• Free streaming of hot dark-matter particles like neutrinos erases structures smaller than 
the asymptotic, comoving free streaming length 

Afg ~ ll(^7o/i^)"^Mpc. 

The hot dark matter spectrum therefore decreases exponentially for k > (A^)~^. 

• This implies that structure formation in a neutrino-dominated universe proceeds ac- 
cording to the top-down scenario: small structures like galaxies form later than large 
structures like galaxy clusters, in conflict with observational results. 

Note: There is a wealth of literature on the subject briefly reviewed here. For further reading, 
and for references on details, I can particularly recommend the following text-books: 
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Neutrino Experiments with Cryogenic Detectors 
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In general an ideal cryogenic thermal detector consists of an absorber for radiation or particles 
where the energy released has to be rapidly converted into thermal energy, changing the 
detector temperature. The absorber is coupled to a thermometer which responds to a change 
of temperature by a change of resistance. In the case of current biasing of the thermistor, 
the temperature change is measured by observing the voltage drop across it. In order to have 
an increase of the detector temperature which is not negligible (in the range of /xK) for a 
very small amount of energy release (100 eV), the detector heat capacity has to be extremely 
small. Since heat capacities decrease at low temperature, the sensitivity is enhanced by 
running the detector at very low temperatures (^ 1 K). In the case of cryogenic thermal 
detectors characterized by a small mass and a heat capacity of the order of 10~^^ J/K at 
0.1 K, it is theoretically possible to obtain an energy threshold and an energy resolution 
of a few eV. The intrinsic slowness of the developed low temperature detectors (total pulse 
duration is about 100 ms) limits their application only to low activity measurements. 

A possible application of this kind of detector is the measurement of the (5 calorimetric 
spectrum with a high energy resolution. The investigation of the (3 decay spectrum of par- 
ticular isotopes, such as ^^'^Re and ^^'''Ho can provide more stringent kinematics limits on 
anti-neutrino and neutrino mass respectively. The isotope ^^^Re is a j3~ emitter; it has the 
lowest value of the end-point Q among the known (5~ isotopes. The study of the /3~ spectrum 
of ^^''Re can be an alternative to the tritium experiment to set a limit on the anti-neutrino 
mass. The effect of a finite anti- neutrino mass should manifest itself as a count deficit, in 
first approximation proportional to m^, near the spectrum end-point. The isotope ^^"^Re is 
naturally occurring at the 62% level and the element Re is a superconductor at the detector 
working temperature of 0.08 K. Therefore it is possible to use a crystal of natural Rhenium as 
a detector absorber to perform the ^^''Re spectrum calorimetric measurement. As a thermistor 
an heavily doped germanium thermistor has been used. 

In the preliminary measurement ^ of the ^*^Re beta spectrum performed in Genoa, an 
energy resolution of 30 eV FWHM has been obtained. The end-point energy, obtained from 
the linear fit of the Kurie Plot in the interval 112-2360 eV, is £^end-point = (2482 ± 12) eV. 
The detector energy calibration has been made using a multi-lines fiuorescence x-ray source; 
the Kq and x-ray lines of CI, Ca, Va, Mn have been detected, in the energy range 100- 
7000 eV (see Fig. 1). Using this external source it has been possible to investigate the detector 
linearity and the energy dependence of the detector energy resolution. A fit of the measured 
energy lines versus the expected values with the function y = ax^ + bx + c has given the 
following result: a = (-0.7 ± 0.7)^-^eV~\ b = (1.000 ± 0.005), c = (2.2 ± 1.5) eV. The 
rms value of the energy resolution has been plotted against the energy lines, and a fit of 
this plot with the function y = Ax + B has given the following result: A = (0.9 it 2.5) E~^, 
B = (15.3 lb 1.4) eV. It should be possible to set a kinematics limit to the anti-neutrino 
mass by means of this calorimetric method lower than 20 eV/c^ in 50 days of data taking. A 
high-statistics measurement of the ^^"^Re spectrum is under development. 




Figure 1: Calorimetric [i spectrum of Re plus a multi-lines X-ray fluorescence spectrum 
due to an external calibration source (see text). 

The isotope ^^^Ho decays by electron capture (EC) with the subsequent emission of an 
electron neutrino. The EC decay is followed by the emission of Auger electrons and X-rays 
due to the electronic cloud excitation of the daughter atom. For each occupied energy level 
from which the electron capture is energetically allowed, there is a line in the calorimetric 
energy spectrum. The effect of a finite neutrino mass should manifest itself as a variation of 
the ratios of the relative EC probabilities. Measuring lines of the calorimetric spectrum, 
it is possible to set N — 1 simultaneous constraints on the neutrino mass and the end-point 
value Q from the relative ratios of the integral value of each spectral line. 

Enclosing the ^^^Ho source in a suitable absorber of a micro-calorimeter it has been 
possible to measure [§, ^ four lines of the calorimetric spectrum (see Fig. 2. a), but the 
statistics of this measurement is not high enough to set a significant limit on the neutrino 
mass. Considering a zero value of the neutrino mass, the end-point energy has been calculated 
by a fit procedure and it is turned out to be £^end-point = (2800 ± 50) eV. Taking into account 
this end-point value it should be possible to set a kinematical limit to the neutrino mass by 
means of this calorimetric method lower than 170 eV/c^ in 50 days of data taking. Running 
more than one detector at the same time it will be possible to further improve the achievable 
mass limit. Although the limit on neutrino mass could not be competitive with the achievable 
limit on electron anti-neutrino mass, we believe it is meaningful to set independent limits on 
the masses of the electron neutrino and anti-neutrino. 

Another important application of this kind of detector is the measurement of the solar 
neutrino flux in radiochemical solar neutrino experiments. In these experiments the inte- 
grated value of the solar neutrino flux is measured counting the extracted isotope which is 
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produced in the target by neutrino capture. The use of a detector characterized by a low 
energy threshold and a good energy resolution can improved the results obtained so far with 
proportional counters. 

A prototype of a new radiochemical lithium detector is now under development at the 
Institute for Nuclear Research at Moscow The interest in this new experiment is related 
to the fact that it could provide new data about the solar neutrino flux at intermediate 
energy. Even taking into account the hypothesis of a strong attenuation of the neutrino flux 
of intermediate energies, this experiment can measure the flux with an accuracy of 12% in 
only one year of measurement time. 

The full scale experiment will have a target mass of only 10 tons, which is the absolutely 
minimal mass among all detectors of solar neutrinos. This experimental scale will be possible 
only if ''Be can be detected with an efficiency close to 100%. The counting of the "^Be isotope 
can be performed measuring the electron capture decay of this isotope to "^Li, but in the 
90% of the cases ''Be decays into ^Li ground state with the subsequent emission of a very 
low detectable energy, of the order of 100 eV. The only possible way to solve this problem 
appears to be the use of cryogenic detectors for "^Be counting. The detection of the ''Be EC 
decay by means of cryogenic detector with an efficiency close to 100% has been successfully 
realised g, |. 




Energy (eV) Energy (eV) 



Figure 2: (a) Calorimetric EC spectrum of a ^^^Ho source enclosed in a superconductive 
tin absorber, (b) Calorimetric ''Be EC spectrum. The source has been produced by proton 
irradiation of a BeO absorber. 

As an absorber a beryllium oxide crystal of mass 3 fig has been used, as this is one of 
the possible forms in which Be can be chemically extracted form the Li target. The energy 
resolution has been found to be 24 eV FWHM for the 112 eV ^Be line (see Fig. 2.b), which 
would enable one to discriminate the background of the counting system in the lithium Solar 
Neutrino Experiment. 
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OMNIS — A Galactic Supernova Observatory 

P.F. Smith 

Rutherford Appleton Laboratory, Chilton, Oxfordshire, 0X11 OQX, UK 

A type II (or lb) supernova explosion releases most of its energy as neutrinos and antincutrinos 
with a time constant of a few seconds, with all three neutrino types and their antiparticles 
produced in comparable numbers. The emission time profile from a typical model is shown 
in Fig. 1 [1]. This provides a unique opportunity to discover neutrino properties which are 
difficult or impossible to determine using terrestrial neutrino sources. In particular a non-zero 
neutrino mass will alter the time profile arriving at the earth (Fig. 2), allowing direct time- 
of-flight measurement of the mass of at least one neutrino type. In particular the most likely 
distance range for Galactic supernovae (~ 2-20 kpc) is ideal for time-of-flight measurement 
of a "cosmologically significant" neutrino mass — i.e. a mass in the range 10-100 eV, for which 
it would form a major component of the mass of the universe and could be a candidate for 
the Galactic dark matter. 
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Figure 1: Time profiles of neutrinos emitted by supernova explosion. The mu and tau neu- 
trinos have identical profiles when emitted. 




Figure 2: Effect of non-zero mass on mu/tau neutrino arrival time profile from distance 8 kpc. 
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The frequency of Galactic supernovae is uncertain. There is a substantial discrepancy 
between the expectation of 1 per 30-100 years, from various astrophysical estimates, and the 
historical record which indicates a higher frequency [2, 3]. Figure 3 shows the locations of 
known supernovae in our Galaxy during the past 1000 years. Only those within 4 kpc of 
the sun are visible optically with high efficiency. A prediction of 1 in 50 years for the whole 
Galaxy would imply an expectation of only 1 within 4 kpc of the sun, whereas in fact 4 or 5 
type Il+Ib have been observed within this radius in the past millennium. This observation 
appears to exclude, with 90% confidence, an interval as large as 30 years or more, and is 
consistent with a total Galactic rate of 1 in 15 it 5 years. 



R = 4 kpc 

(5% sample of Galaxy) 




Figure 3: Galactic map of recorded supernova for past 1000 years, showing that only the 
closest 5% are visible optically. The much larger number from the whole Galaxy would be 
observable as neutrino bursts. 

With a time interval of this order, it becomes reasonable to set up an underground de- 
tector array for such an event, run on a part-time basis in parallel with other astrophysics 
and particle physics projects. For a statistical precision sufficient to determine a neutrino 
mass in the 10-100 eV range, we require about 2000 events (1000 each of mu and tau neutri- 
nos). A number of existing world detectors are sensitive to supernova neutrinos, in particular 
Super-Kamiokande, LVD, MACRO, and SNO. However, these detect principally the electron 
antincTitrino component throTigh charged current interactions, and have relatively low sensi- 
tivity to neutral current events [4]. Numerical estimates in Table 1 show that from an 8 kpc 
supernova Super-Kamiokande would register ~ 8000 charged current events, but only 1% 
of that number from neutral currents. We describe here a low cost detection scheme which 
would be sensitive principally to the mu/tau neutrino component. This would provide the 
first direct observation of a cosmologically significant neutrino mass, and would complement 
the electron antineutrino signal from other world detectors to provide data on mixing between 
all three neutrino generations. Information on the profile shape for all three neutrino types 
would be also of considerable astrophysical interest. 
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The proposed detector array is referred to as OMNIS — an observatory for multiflavour 
interactions from Galactic supernovae. It has evolved from the Supernova Neutrino Burst 
Observatory (SNBO) proposed by D.B. Cline et al., based on neutrino detection by neutral 
current nuclear excitation in natural underground rock [5]. The excited nucleus would emit 
neutrons, which could be captured by counters embedded in the rock. The E2 dependence of 
the excitation cross section makes the detector sensitive principally to the higher temperature 
mu/tau neutrinos, producing a dominant neutral current signal. The evolution of OMNIS 
from SNBO is described in a more detailed paper [6]. There are two major improvements: 
(a) neutron collection efficiencies higher by an order of magnitude axe obtained by siting 
optimised detectors in open caverns, and (b) additional target materials are used, in particular 
Fe and Pb, which give higher neutron production rates and with Pb also having a charged 
current excitation route [7, 8, 9]. Thus differences between lead, iron and rock signals can be 
used to infer mixing between mu/tau and electron neutrinos. 

A possible arrangement of detector /target modules for OMNIS is shown in Fig. 4. Stan- 
dard Gd or Li-loadcd scintillator could be used for the neutron detection, in modules of about 
0.5-1 ton. About 200-300 modules would be needed to observe 2000 events from an 8 kpc 
supernova. This demonstrates the uniquely high efficiency of the detector array, giving it 
a unit cost much lower than other world neutrino detectors. An average of over 10 events 
per ton scintillator could be achieved using indirect neutron detection compared with the 0.3 
events per ton scintillator achieved with direct charged current detection in MACRO. The 
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Figure 4: Possible arrangement of target materials and neutron detectors for OMNIS. 
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Table 1: Comparison of proposed OMNIS multi-target observatory with world detectors based 
on direct detection with water and scintillator targets. Approximate event numbers for each 
neutrino type shown for a supernova at 8 kpc. 
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estimated neutron production rates in salt rock, iron, and lead are 0.06, 0.15, 1 per ton for a 
supernova at 8 kpc. In addition, higher collection efficiency is possible from Fe and Pb due 
to lower absorption [6], giving them an effectiveness in the ratio 1:4: 30. Thus lead is by far 
the most effective target to maximise event numbers, but the design must include a sufficient 
quantity of the other targets to provide information on mixing between mu/tau and electron 
neutrinos, based on the fact that this would enhance the charged current excitation in Pb 
but not in Fc or rock. 

Tabic 1 shows the event numbers of different neutrino types from principal world de- 
tectors, compared with those from the proposed OMNIS detector. Note that the neutron 
background from U/Th in the rock would not be significant compared with the signal burst, 
and for sites with depths > 600 m, the neutron background from muons would also be an 
order of magnitude below the signal. Thus the project is technically straightforward and 
based on conventional low-cost neutron detection technology. Neutron background is dis- 
cussed in more detail in [6] where some discussion is also included on the enormously greater 
background problem for extra-Galactic supernova detection. For discrimination of supernova 
events from local neutron showers from cosmic ray muons, it is desirable to split the array 
between several underground caverns, and also between several world sites, also permitting 
estimates of supernova direction by arrival time differences. 

In conclusion the relatively low unit cost of the OMNIS scheme makes it possible to propose 
a dedicated Galactic supernova neutrino observatory, capable of observing the time profile 
of the mu and tau neutrino component. This would provide the first direct measurement of 
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a cosmologically-significant neutrino mass, and the use of several different target elements, 
in conjunction with the large electron antineutrino signal from SuperKamiokande would also 
provide data on the presence or absence of neutrino mixing on the Galactic distance scale. 
It could be designed to run largely unattended, with minimal maintenance, so would be a 
part-time activity for the participating groups, alongside other underground experiments. 
The current collaboration on this project involves groups from RAL, Manchester, Sheffield, 
Imperial College (for the UK Boulby Mine site) and UCLA, UCSD, Ohio, LLNL, LANL, 
UT-Dallas, PNL, for the US Carlsbad site, both sites in low activity salt rock and already 
containing suitable infra-structure installed for other purposes. First prototype modules could 
now be set up in these sites and funding is being sought for this. New collaborators would be 
welcome to strengthen the case for construction of the full array. 
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Theory of Neutrino Mass 

Neutrinos are the only apparently massless electrically neutral fermions in the Standard Model 
(SM). Nobody knows why they are so special when compared with the other fermions. The 
masslessness of neutrinos is imposed in an ad hoc fashion, not dictated by an underlying 
principle such as gauge invariance. Moreover, massless neutrinos may be in conflict with 
present data on solar and atmospheric neutrino observations, as well as cosmological data on 
the amplitude of primordial density fluctuations |l| . The latter suggest the need for hot dark 
matter in the Universe. On the other hand, if massive, neutrinos would pose another puzzle, 
namely, why are their masses so much smaller than those of the charged fermions? The key 
to the answer may lie in the fact that neutrinos could be Majorana fermions. Such fermions 
are the most fundamental ones. If neutrinos are majorana particles the suppression of their 
mass could be related to the feebleness of lepton number violation. This, in turn, would point 
towards physics beyond the SM. 

Lepton number, or B — L can be part of the gauge symmetry Q or, alternatively, can be 
a spontaneously broken global symmetry. In the latter case there is a physical pseudoscalar 
Goldstone boson generically called majoron j^. One can construct an enormous class of such 
models Q which may have important implications not only in astrophysics and cosmology 
but also in particle physics 

Many mechanisms exist to explain the small masses of neutrinos as a result of the violation 
of lepton number. The first is the seesaw mechanism, which is based on the existence of some 
relatively large mass scale. However, neutrinos could acquire their mass radiatively. In this 
case the smallness of the mass holds even if all the new particles required to generate the 
neutrino masses are light and therefore accessible to present experiments . The seesaw and 
the radiative mechanisms of neutrino mass generation may be combined. Supersymmetry 
with broken R-parity provides a very elegant mechanism for the origin of neutrino mass 
1^] in which the tau neutrino i^r acquires a mass due to the mixing between neutrinos and 
neutralinos. This happens in a way similar to the seesaw mechanism, in which the large mass 
scale is now replaced by a supersymmetry breaking scale characterizing the neutralino sector. 
In supergravity models with universal soft breaking scalar masses, this effective neutralino- 
neutrino mixing is induced only radiatively Q]. As a result the Majorana mass is naturally 
suppressed, even though there is no large mass scale present. From this point of view, the 
mechanism is a hybrid between the see-saw idea and the radiative mechanism. Despite the 
small neutrino mass, many of the corresponding R-parity violating effects can be sizeable. An 
obvious example is the fact that the lightest neutralino decay will typically decay inside the 
detector, unlike the case of the minimal supersymmetric model. 

Other than the seesaw scheme, none of the above models requires a large mass scale. In all 
of them one can implement the spontaneous violation of the lepton number symmetry leading 
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to neutrino masses that vanish as the lepton-number scale goes to zero, in contrast to the see- 
saw model. Such low-scale models are very attractive and lead to a richer phenomenology, 
as the extra particles required have masses at scales that could be accessible to present 
experiments. One remarkable example is the possibility of invisibly decaying Higgs bosons 

The large diversity of neutrino mass schemes and the lack of a theory for the Yukawa cou- 
plings imply that present theory is not capable of predicting the scale of neutrino masses any 
better than it can fix the masses of the other fermions, such as the muon: one should turn to 
observations as means for constraining their properties. Direct laboratory observations leave 
a lot of room for massive neutrinos, except beta decay experiments which are more restrictive. 
A complementary approach to the problem of neutrino mass tries to get information from as- 
trophysics and cosmology. It constitutes a promising interdisciplinary field, now often-called 
astroparticle physics. In this talk I will briefly mention recent work devoted to the constraints 
that one can place on neutrino properties from cosmological and astrophysical observations. 



Limits from Cosmology 

The first cosmological bound on neutrino masses follows from avoiding the overabundance of 
relic neutrinos 

^m,, < 921^,/i2eV, (1) 

where Q^h'^ < 1 and the sum runs over all stable species of isodoublet neutrinos with mass 
less than 0(1 MeV). Here = Pu/Pc, where pi, is the neutrino contribution to the total 
density and pc is the critical density. The factor /i^ measures the uncertainty in the present 
value of the Hubble parameter, 0.4 < h < I, and is smaller than 1. For the and u-j- 

this bound is much more stringent than the laboratory limits. 

Apart from the experimental interest, a heavy tau neutrino (say in the MeV range) could 
also be interesting from the point of view of structure formation Q. According to eq. (|l|) 
such neutrino can not exist if it has only the interactions prescribed by the SM. However 
a heavy tau neutrino is theoretically viable if in models with spontaneous violation of total 
lepton number since these contain new interactions of neutrinos majorons which may cause 
neutrinos to decay into a lighter neutrino plus a majoron, for example |^, 

Ur^iy^ + J. (2) 

or have sizeable annihilations to these majorons, 

zy^ + i/^ — > J + J . (3) 

The possible existence of fast decay and/or annihilation channels could eliminate relic neutri- 
nos and therefore allow them to have higher masses, as long as the lifetime is short enough to 
allow for an adequate red-shift of the heavy neutrino decay products. These 2-body decays 
can be much faster than the visible modes, such as radiative decays of the type u' ^ v + 'j. 
Moreover, the majoron decays are almost unconstrained by astrophysics and cosmology (for 
a detailed discussion see ref. @). 

A general method to determine the majoron emission decay rates of neutrinos was first 
given in ref. [^]. The resulting decay rates are rather model-dependent and will not be 



discussed here. Explicit neutrino decay lifetime estimates are given, e.g. in ref. [Q, 11]. The 
conclusion is that there are many ways to make neutrinos sufficiently short-lived and that all 
mass values consistent with laboratory experiments are cosmologically acceptable. 
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Figure 1: A heavy annihilating to majorons can lower the equivalent massless-neutrino 
number in nucleosynthesis. 



Stronger limits on neutrino masses, lifetimes and/or annihilation cross sections arise from 
cosmological Big-Bang Nucleosynthesis. Recent data on the primordial deuterium abundance 



[12 1 have stimulated a lot of work on the subject |13]. If a massive Uj- is stable on the 
nucleosynthesis time scale, {vr lifetime longer than ~ 100 sec), it can lead to an excessive 
amount of primordial helium due to their large contribution to the total energy density. This 
bound can be expressed through an effective number of massless neutrino species (N^)- If 



Ni, < 3.4—3.6, one can rule out i/^ masses above 0.5 MeV [14|. If we take N^, < 4 the 
nii^^ limit loosens accordingly. However it has recently been argued that non-equilibrium 
effects from the light neutrinos arising from the annihilations of the heavy i/j- 's make the 
constraint a bit stronger in the large m^^ region [^]. In practice, all Ut masses on the few 
MeV range are ruled out. One can show, however that in the presence of Vr annihilations 
the nucleosynthesis mj^^ bound is substantially weakened or eliminated [|^]. In Fig. |l| we 
give the effective number of massless neutrinos equivalent to the contribution of a massive 
I'r majoron model with different values of the coupling g between i^t-'s and J's, expressed in 
units of 10^^. For comparison, the dashed line corresponds to the SM g = case. One sees 
that for a fixed N^^^, a wide range of tau neutrino masses is allowed for large enough values 
of g. No Vr masses below the LEP limit can be ruled out, as long as g exceeds a few times 
10"^. One can express the above results in the mj^^-^-plane, as shown in Fig. |2|. Moreover the 
required values of gim,^^) are reasonable in many majoron models |l^. Similar depletion 
in massive z/,- relic abundance also happens if the i'r is unstable on the nucleosynthesis time 



scale 1 18] as will happen in many majoron models. 



Limits from Astrophysics 

There are a variety of limits on neutrino parameters that follow from astrophysics, e.g. from 
the SN1987A observations, as well as from supernova theory, including supernova dynamics 
|jl9| and from nucleosynthesis in supernovae [^0| . Here I briefly discuss three recent examples 
of how supernova physics constrains neutrino parameters. 

It has been noted a long time ago that, in some circumstances, massless neutrinos may 



be mixed in the leptonic charged current |21|. Conventional neutrino oscillation searches in 
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Figure 3: SN1987A bounds on massless neutrino mixing. 



vacuo are insensitive to this mixing. However, such neutrinos may resonantly convert in the 
dense medium of a supernova |22|. The observation of the energy spectrum of the SN1987A 
Pe's [23 1 may be used to provide very stringent constraints on massless neutrino mixing angles, 
as seen in Fig. |^ The regions to the right of the solid curves are forbidden, those to the left 
are allowed. Massless neutrino mixing may also have important implications for r-process 
nucleosynthesis in the supernova For details see ref. [ |22| . 

Another illustration of how supernova restricts neutrino properties has been recently con- 
sidered in ref. |24]. In this paper flavour changing neutral current (FCNC) neutrino interac- 
tions were considered. These could arise in supersymmetric models with R parity violation. 
These interactions may induce resonant neutrino conversions in a dense supernova medium, 
both in the massless and massive case. The restrictions that follow from the observed en- 
ergy spectra from SN 1987A and the supernova r-process nucleosynthesis provide constraints 
on R parity violation, which are much more stringent than those obtained from the labora- 
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Figure 4: SN 1987A bounds on FCNC neutrino interactions. 
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Figure 5: Supernovae and FCNC neutrino interactions. R-process. 



tory. In Fig. ^ we display the SN 1987A constraints on explicit i?-parity- violating FCNCs in 
the presence of non-zero neutrino masses in the hot dark matter eV range. In Fig. ^ we dis- 
play the sensitivity of the r-process argument to the presence of FCNC neutrino interactions. 
Taken altogether, Fig. ^ and Fig. |5| disfavour a leptoquark interpretation of the recent HERA 
anomaly. 

As a final example of how astrophysics can constrain neutrino properties we consider 
the case of resonant fe Vs and Ve — > i>s conversions in supernovae, where Vg is a sterile 



neutrino |25], which we assume to be in the hot dark matter mass range. The implications 
of such a scenario for the supernova shock re-heating, the detected i>f> signal from SN 1987A 
and for the r-process nucleosynthesis hypothesis have been recently analysed |25]. In Fig. 
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Figure 6: Supernovae shock reheating argument and — Vg neutrino conversions. 

taken from ref. |^5|, we illustrate the sensitivity of the supernova shock re-heating argument 
to active-sterile neutrino mixing and mass difference. Notice that for the case of r-process 
nucleosynthesis there is an allowed region for which the r-process can be enhanced. 



Acknowledgements 

Supported by DGICYT grant PB95-1077 and in part by EEC under the TMR contract 
ERBFMRX-CT96-0090. It is a pleasure to thank the organizers for a very friendly atmo- 
sphere. 



References 

[1] For reviews on the status of neutrino mass see A. Yu. Smirnov, Neutrino Masses and 
Oscillations, Plenary talk given at 28th International Conference on High energy physics, 



July 1996, Warsaw, Poland, hep-ph 9611465; and Neutrinos Properties Beyond the 



Standard Model, J. W. F. Valle, Plenary talk, WIN97, Capri, Italy, June 1997, ^epj 



ph/9709365 



[2] M Gell-Mann, P Ramond, R. Slansky, in Supergravity, ed. D. Freedman et al. (1979); 
T. Yanagida, in KEK lectures, ed. O. Sawada et al. (1979); R.N. Mohapatra and G. Sen- 
janovic, Phys. Rev. D23 (1981) 165 and Phys. Rev. Lett. 44 (1980) 912 and references 
therein. 

[3] Y. Chikashige, R. Mohapatra, R. Peccei, Phys. Rev. Lett. 45 (1980) 1926 



[4] For reviews see J. W. F. Valle, Gauge Theories and the Physics of Neutrino Mass, Prog. 
Part. Nucl. Phys. 26 (1991) 91-171 



176 



[5] A. Joshipura and J. W. F. Valle, Nucl. Phys. B397 (1993) 105 and references therein 
[6] A. Zee, Phys. Lett. B93 (1980) 389; K. S. Babu, Phys. Lett. B203 (1988) 132 



[7] Marco A. Di'az, Jorge C. Romao and Jose W. F. Valle, [hep-ph/97063T5| ; M. A. Diaz, A 



Joshipura, J. W. F. Valle, in preparation; For additional papers on neutrino masses in 
broken R-parity models see: F. Vissani and A. Yu. Smirnov, Nucl. Phys. B460, 37-56 



(1996) . R. Hempfling, Nucl. Phys. B478, 3 (1996), and |hep-ph/9702412| ; H.R Nilles and 
N. Polonsky, Nucl. Phys. B484, 33 (1997); B. de Carlos, P.L. White, Phys.Rev. D55 
4222-4239 (1997); E. Nardi, Phys. Rev. D55 (1997) 5772; S. Roy and B. Mukhopadhyaya, 
Phys. Rev. D55, 7020 (1997). 

[8] E. Kolb, M. Turner, The Early Universe, Addison- Wesley, 1990, and references therein 

[9] J. Bardeen, J. Bond and G. Efstathiou, Astrophys. J. 321 (1987) 28; J. Bond and G. 
Efstathiou, Phys. Lett. B265 (1991) 245; M. Davis et al.. Nature 356 (1992) 489; S. 
Dodelson, G. Gyuk and M. Turner, Phys. Rev. Lett. 72 (1994) 3754; H. Kikuchi and E. 
Ma, Phys. Rev. D51 (1995) 296; H. B. Kim and J. E. Kim, Nucl. Phys. B433 (1995) 
421; M. White, G. Gelmini and J. Silk, Phys. Rev. D51 (1995) 2669 A. S. Joshipura 
and J. W. F. Valle, Nucl. Phys. B440 (1995) 647. 

[10] J. Schechter, J. W. F. Valle, Phys. Rev. D25 (1982) 774 

[11] J. W. F. Valle, Phys. Lett. B131 (1983) 87; G. Gelmini, J. W. F. Valle, Phys. Lett. 
B142 (1984) 181; J. W. F. Valle, Phys. Lett. B159 (1985) 49; M. C. Gonzalez-Garcia, 
J. W. F. Valle, Phys. Lett. B216 (1989) 360. A. Joshipura, S. Rindani, Phys. Rev. D46 
(1992) 3000; J. C. Romao, J. W. F. Valle, Nucl. Phys. B381 (1992) 87-108 

[12] R.F. Carswell, M. Ranch, R.J. Weynman et al., MNRAS 268 (1994) LI; A. Songalia, 
L.L. Cowie, C. Hogan and M. Rugers, Nature 368 (1994) 599; D. Tytler and X.M. Fan, 
Bull. Am. Astr. Soc. 26 (1994) 1424; D. Tytler, talk at the Texas Symposium, December 
1996. 

[13] N. Hata et al., Phys. Rev. Lett. 75 (1995) 3977; C.J. Copi, D.N. Schramm and M.S. 
Turner, Science 267 (1995) 192 and Phys. Rev. Lett. 75 (1995) 3981; K. A. Olive and 
G. Steigman, Phys. Lett. B354 (1995) 357-362; S. Sarkar, Rep. on Prog, in Phys. 59 

(1997) 1493; P. J. Kernan and S. Sarkar, Phys. Rev. D 54 (1996) R3681 

[14] E. W. Kolb, M. S. Turner, A. Chakravorty and D. N. Schramm, Phys. Rev. Lett. 67 
(1991) 533; A.D. Dolgov and I.Z. Rothstein, Phys. Rev. Lett. 71 (1993) 476. 

[15] A.D. Dolgov , S. Pastor, and J.W.F. Valle, Phys. Lett. B383 (1996) 193-198; §ep] 
I ph/9602233| . For related papers see Hannestad and Madsen Phys. Rev. Lett. 77 (1996) 

5148, Phys. Rev. D54 (1996) 7894; J. Rehm, G. Raffelt and A. Weiss Astron. Astron. 

& Astrophys. 327 (1997) 443; B. Fields, K. Kainulainen and K. Olive, Astropart. Phys. 

6 (1997) 169 

[16] A.D. Dolgov, S. Pastor, J.C. RomSo and J.W.F. Valle, Nucl. Phys. B496 (1997) 24-40 



| hep-ph/9610507 | 



[17] A Masiero, J. W. F. Valle, Phys. Lett. B251 (1990) 273; J. C. Romao, C. A. Santos, 
and J. W. F. Valle, Phys. Lett. B288 (1992) 311 



177 



[18] For a recent ref. see S. Hannestad, hep-ph/971124£ ; for a review see G. Steigman; in 
Cosmological Dark Matter, p. 55, World Scientific, 1994, ISBN 981-02-1879-6 

[19] G. Raffelt, Stars as laboratory for fundamental physics, Univ. of Chicago Press, Chicago- 
London (1996). 

[20] Y.-Z. Qian, G. M. Fuller, G. J. Mathwes, R. W. Mayle, J. R. Wilson and S. E. Woosley, 
Phys. Rev. Lett. 71 (191965) 93. 

[21] J.W.F. Valle, Phys. Lett. B199 (1987) 432 

[22] H. Nunokawa, Y.Z. Qian, A. Rossi, J.W.F. Valle, Phys. Rev. D54 (1996) 4356-4363, [ 



hep-ph/9605301 | 



[23] A. Yu. Smirnov, D. N. Spergel and J. N. Bahcah, Phys. Rev. D49 (1994) 1389. 

[24] H. Nunokawa, A. Rossi and J. W. F. Valle, Nucl. Phys. B482 (1996) 481. 

[25] H. Nu nokawa, J.T. Pe ltoniemi, A. Rossi, J.W.F. Valle, Phys. Rev. D56 (1997) 1704- 
1713; |iep-ph/9702372 l 



178 



Some Neutrino Events of the 21st Century 

Leo Stodolsky 

Max- Planck- Institut fiir Physik, Fdhringer Ring 6, 80805 Miinchen, Germany 

For this short concluding talk about the future of neutrino physics, I would like to give the 
only talk of the meeting which has a really good chance of being 100 percent wrong. Let us 
imagine ourselves — or our grandchildren — here a hundred years from now. We're reviewing 
neutrino occurences of the last century. Sticking our necks way out, here are some guesses as 
to what happened. 

Present— 2015: Masses and Mixings 

In this period, we may suppose the various masses and mixings of the neutrinos will be 
elucidated — or at least the larger ones where the oscillation length is not greater than the 
earth-sun distance and the mixing angle is not microscopic. I recall that a 30 eV neutrinos 
once was the natural candidate for dark matter; so this point, that is if there are neutrino 
masses in the eV range, is important in telling us if neutrinos play any role in cosmology 
following the big bang. 

2000-2025: Point Sources 

In this period the existence of high energy point sources of neutrinos should be established, 
by very large detectors in the water or in the ice. If the recent spectacular observations Q of 
point sources in high energy photons have a hadronic mechanism as their origin, then there 
should be an equal number of similarly energetic neutrinos. That is, if we have a hadronic 
accelerator making pions then the photons come from vr*^ decay and the neutrinos from vr^'s, 
giving v/j ~ 1. On the other hand the photons could be originating from something non- 
hadronic, such has high energy electrons. Thus the neutrino observations will tell us about 
the nature of the production mechanisms in these high energy sources. 

2010: Coherent Scattering in the Lab 

The coherent scattering of neutrinos is an old favorite topics. Due to the (neutron number)^ 
factor in the cross section, such high rates result ^ that one can imagine smallish neutrino 
detectors for many interesting applications. This process is totally uncontroversial in the 
Standard Model. But it still would be nice to see it in the laboratory, which should be 
possible with cryogenic detectors early in the century. 

2023: Next Nearby Supernova 

All neutrino folk want to know when the next Milky Way supernova is going to be. Let's put 
our money on 2023. By then several big detectors should be operating smoothly, allowing a 
full exploitation, of the thousands of events they will observe, both in neutral and charged 
current modes. 
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2030: Neutrino Thermometer 

Now we're getting to some harder stuff, real challanges for the next generation of experimen- 
tahsts. As you know, the ''Be neutrinos from the sun are essentially monochromatic, coming 
from an electron capture process. However the Be nuclei in the center of the sun are moving 
around due to the high temperature, giving a spread and shift of the line. A calculation Q| 
gives the shift as 1.29 keV. A measurement, performed perhaps by comparing with a terres- 
trial source, would give the temperature in the center of the sun, an amusing use of neutrinos 
as a thermometer. 

2040: Neutrino Geology 

Around this time we see the active development of neutrino geology. The various radioactive 
processes in the earth give rise to neutrinos (actually anti-neutrinos, except for the neutrino 
"line" from electron capture in ^^K). Measurement of these terrestrial neutrinos will give 
a direct snapshot of the energy production in the earth and allow us to answer many in- 
teresting questions of geophysics and the thermal history of the earth. Since the coherent 
superconducting detector provides a light and portable instrument, local and "tomographic" 
studies will be possible, as well as the investigation of the planets and their moons. I set this 
relatively late since the background from solar neutrinos and probably also nuclear reactors 
must be well understood to see this signal (see Fig. 11 of Ref. [2]). 

2050: Extra-Galactic Neutrino Burst Observatory 

Now we're getting to the big stuff: the extra-galactic neutrino bursts. If we had a detector that 
could see the neutrinos from stellar collapse in nearby clusters of galaxies we wouldn't have 
to wait decades for the next event. With a thousand galaxies in the Virgo cluster at 10 Mpc, 
we will be having them every few weeks and supernova neutrino observations will become a 
systematic affair. This will permit the study of many interesting points, such as flavor "echos" 
due to mixing, tests of CP for neutrinos and so forth Q. The enormous distance involved 
means various time-of-flight effects due to neutrino masses are greatly magnified (in fact there 
is a danger this could become too much of a good thing since too much spreading of the pulse 
could make it disappear into the background |^]). Such a detector, or perhaps we should say 
observatory, is very ambitious, but not inconceivable. With the coherent scattering process 
and cryogenic detection to see the small recoils, about a megaton of cold material could 
suffice [^. An alternative idea is the "magic mountain" or OMNIS Q where neutral-current 
induced neutrons from natural Ca are the signal. The joker in these proposals, as usual, is 
background. 

2050: Neutrino Technology 

It took the laser about 50 years to get from the lab to the checkout counter at the supermarket. 
Since neutrinos are more difficult than photons, let's give them a hundred years from their 
discovery to get into the economy. About mid-century, then, the light and portable detector 
will allow us to monitor nuclear power stations from the outside and to make geological 
investigations for minerals and petroleum. Here again an understanding of the background 
is essential. I'm not so sure about the sometimes mentioned neutrino telecommunication 
channel, because I don't know what the transmitter is supposed to be. 
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2099: Relic Neutrinos 

And last but not least, we have the observation of the big bang neutrinos. This at the end of 
the century since I don't know how this is going to be done. These relic neutrinos are a simple 
consequence of big bang nucleosynthesis and should exist just as the relic microwave photons 
do, with about the same density: 100/cm'^ for each species. While their existence is not at 
all controversial, it would be very satisfying to observe them and study their properties. The 
great difficulty is their low energy, (actually momentum, which is the correct thing to redshift 
if they are massive Q) which is only millivolts. Previous schemes for low energy neutrino 
detection involving some kind of energy transfer would seem to be unworkable, leading one 
to consider coherent effects like forces and torque 0. Two schemes have survived various 
criticisms. In one, a net helicity in space because of an excess of one kind of neutrino (not 
a feature of big bang cosmology) would exercise a torque on an object with spin, in another 
a sponge-like object with pores on about the wavelength of the neutrinos would get random 
kicks from scattering, leading to a hopefully detectable force. A variant of this would be to 
observe myriads of little grains of this size and wait til one of them hops Or, on another 
tack, there is the possibility that extremely high energy neutrinos, near the cosmic ray cutoff, 
hitting the background neutrinos and making a Z boson, lead to an observable characteristic 
cascade None of these ideas seem particularly practical, but then 2099 is a long way off 
and maybe our grandchildren won't be so dumb . . . 

Comments 

Some people have objected that I didn't include any results from accelerator physics like 
observation of right-handed currents or new heavy neutrinos. True, but I don't know what 
these may be. I stress that I have interpreted the word "neutrino" narrowly, as the familiar 
known, almost massless, object. 

Others have suggested I am much too conservative and are willing to bet that it's all 
going to happen more and faster. That would be great, let's hope so. Indeed, concerning the 
terrestrial neutrinos the Borexino group ||l^ says they are seriously considering looking for 
an anti-neutrino signal, which could be detectable according to estimates of Raghavan. 
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